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Introduction 


The report about the fact that the Kovalevskaya case of the dynamics of a rigid body can be 
generalized to a gyrostat was made by Yehia at Demin’s seminar and Kozlov’s seminar at the Moscow 
State University in 1985. At the same time, Yehia presented an article in Journal de Mechanique 
Theoretique et Applique which, for unknown reasons, was not published. In this regard, the official 
opening date of this integrable cases is 1986 when the paper [69] was published. The Russian version of 
this paper [54] was presented much later, in April 1986. In fact, in [69], the Kovalevskaya integral was 
generalized in two ways: to a gyrostat and to a double field, which simulates the action of the gravity 
and a constant magnetic field. Earlier, an analogue of the Kovalevskaya integral for a double field was 
found by Bogoyavlensky [5]; however, Yehia’s generalization turned out to be crucial: the introduction 
of the gyrostatic moment violated the classical structure of the integral (the sum of squares) and its 
homogeneity: the new summand, proportional to the gyrostatic moment, has the third degree in the 
angular velocity, similarly to the Goryachev-Chaplygin integral. In 1987, the papers [59, 64] appear; 
they contain a generalization of the the Kovalevskaya integral to a gyrostat in a uniform field that 
simulates the gravity force, but these publications cannot be considered original. 

Since, in the general case, the introduction of a double field breaks the symmetry of the problem 
(the area integral becomes invalid), Yehia noted the following two cases of complete integrability: a 
gyrostat of the Kovalevskaya type in the gravity field and a gyrostat in the double field of a specific 
structure that admits a singular symmetry. In the present review, we we consider the first problem. 
Until now, it has not been reduced to quadratures. However, it surprisingly turned out that all motions 
on critical manifolds of the integral mapping that play a key role in the topological analysis of the 
problem had been detected and integrated long before the paper [69] and the emergence of the modern 
approach to the study of integrable Hamiltonian systems. In 1991, scientists began a systematic study 
of the phase topology of the Kovalevskaya-Yehia case. The first step of the study is the examination of 
singularities of the integral mapping (in modern literature, it is usually called the moment mapping). 
These singularities are trajectories on the phase space on which the first integrals are dependent in 
the sense of the linear dependence of differentials. The corresponding solutions of the Euler-Poisson 
equations are either fixed points (from the physical standpoint, they correspond to uniform rotations 
of the body) or periodic trajectories of a special nature, on which the rank of the moment mapping 
decreases. This subset of the phase space is a critical set of the moment mapping. In [25], a system 
of equations that determines the critical set was obtained; considering sections of this set and taking 
into account obvious symmetries, one can obtain the parametric equations of the surfaces that contain 
critical values of the moment mapping (the bifurcation diagram of the problem). We emphasize 
the connection of the bifurcation diagram with Appelrot classes of the classical Kovalevskaya top 
(see [4]). It turned out that surfaces in the space of constants of first integrals corresponding to the 
four Appelrot classes in the Kovalevskaya case are transformed into two surfaces in the Yehia case at 
nonzero gyrostatic moment; one of these surfaces consists of two connected components. At the same 
time, in the paper [11], all solutions of the Euler-Poisson equations that lie in the preimage of the 
first surface have been reduced to quadratures. 

It was stated in [12] that solutions that correspond to singular periodic solutions (i.e., consist of 
critical points of the moment mapping, see [25]) have been obtained and studied earlier in [36] as a 
part of the class of motions of a heavy gyrostat (this class is a generalization of the Bobylev-Steklov 
solution). Further, in [13, 26] the exhaustive study of the structure of the bifurcation set was proposed; 
It was noted in [13] that the equations of critical points in the preimage of the second surface coincide 
with the invariant relations that describe a particular solution detected in [37]. A more accurate 
analysis shows that to describe all critical motions in the preimage of the second bifurcation surface, 
it is necessary to consider the set of solutions found in [37, 51]. Thus, the whole set of critical motions 
of the Kovalevskaya-Yehia top is an open set described by quadratures in [36, 37, 51] without analogs 
of the Kovalevskaya integral that were unknown at that time. 
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Note that the bifurcation diagram of the moment mapping depends on a significant physical pa¬ 
rameter (namely, the gyrostatic moment), and, therefore, the reduced systems with two degrees of 
freedom, obtained by the reduction of the order, depend on two parameters (the gyrostatic moment 
and the area constant). In [27, 28], the classification problem for bifurcation diagrams for reduced 
systems was solved and the separating set in the plane of these parameters was constructed. Possible 
variants of classification of topological invariants (Fomenko graphs), still without complete analytical 
classification, were considered in [14, 46]. 

In the present paper, we make an attempt of systematic treatment (from the modern point of view) 
of existing analytical and qualitative results. The idea of the work belongs to M. P. Kharlamov who 
devised it in 2011 as a review article for the 25th anniversary of the discovery of this integrable case. 
However, it was found that numerous papers treating these topics contain a lot of discrepancies, errors, 
and gaps in proofs. Therefore, we continued our work as an analytical research with the aim of giving 
an exhaustive description of the phase topology with a complete argumentation. This leads us to 
another series of publications supported by the Russian Foundation for Basic Research in 2011-2015. 
All these results are presented below. A. Yu. Savushkin and E. G. Shvedov participated in the research 
and development of computer visualization programs of the objects considered; they created a software 
package in the Mathematica system, which can be called a constructor of topological invariants. 

The research has been continued for several years under the support of the Russian Foundation for 
Basic Research (project Nos. 10-01-00043, 10-01-97001, 13-01-97025, 14-01-00119, and 15-41-02049). 

We do not include a detailed presentation of results regarding the two problems that can be con¬ 
sidered as degenerate cases of the general case. First, this is the classic Kovalevskaya problem (with 
zero gyrostatic moment) and, second, the Kovalevskaya-Yehia gyrostat in the case where the constant 
cyclic integral (the are constant) vanishes. The rough topology of the first problem is completely 
described in [38-40]. In [23], the sequence of bifurcations described in [39] was combined in Fomenko 
graphs on isoenergetic levels and a complete list of these graphs was found (i.e., a classification of 
graphs in the energy-moment urn plane was proposed). Marks on graphs (fine topological invariants or 
Fomenko-Zischang invariants) are calculated in the paper [7], which also contains a consistent presen¬ 
tation of the method of loop molecules in subtle topological analysis. Topological analysis of a family 
of systems that occur for the zero area constant (the free parameter is the value of the gyrostatic 
moment) is described in a series of papers [25-29]. Marks on graphs (fine topological invariants or 
Fomenko-Zischang invariants) have also been calculated in [20, 21], where a number of general asser¬ 
tions important for applications of the method of loop molecules were proved. We appeal to these 
problems only when necessary, if we want to generalize any facts. 

1. Analytic Results 

1.1. Equations and integrals. In the general case, the Euler-Poisson equations of motion of a 
gyrostat in the gravitational field have the form 

M = Mxu;-|-cxq!, a = a x u>, (1.1) 

where uo is the angular velocity, a is the unit vector of the direction of the gravitation force (it is 
called the vertical), and c is the vector connecting the immovable point O of the body and its center 
of mass whose magnitude is equal to the product of the weight of the system “body+rotor” by the 
distance between O and the center of mass. All objects are assigned to the movable axes. The vector 
of kinetic moment M is related with the angular velocity by the formula 

M = Ilj + A, (1.2) 

where I and A are the tensor of inertia at the point O and the vector of gyrostatic moment, respectively; 
they are constant in the movable coordinate system. For brevity, objects that are constant in the 
movable coordinate system are said to be “constant in the body.” The system (1.1) is obtained from 
the mechanical system with the configuration space SO( 3) (a Hamiltonian system with three degrees 
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of freedom) by factorization with respect to the action of the group of rotations about a fixed (in 
the inertial space) vertical direction; the phase space of this system is the five-dimensional manifold 
P 5 =R 3 {u)xS 2 (a). 

Here and below, we denote by M n (rti,..., u n ) the space 
functions u\,... ,u n . We denote by S 2 (a) the unit sphere in 
The definitions above imply the relation 


with a fixed system of coordinate 
ex) = M 3 (ai, a 2 , a 3 ). 


called the geometric integral, where 


r = 1 


r — Gy T g2 T G3. 


(1.3) 


Other general integrals are the Hamiltonian (the total energy) 


H = —Iuj 
2 


Ul) — c • a 


and the cyclic integral (the area integral) 


L = -M • a. 


(1.4) 


(1.5) 


The coefficient 1/2 is introduced following by S. V. Kovalevskaya. The Poisson bracket of this me¬ 
chanical system on the space {(M,cx)} can be represented by its values on pairs of the coordinate 
functions (see, e.g., [5]): 

{ -1 A ■ 1 / / } — -y’/A'-1 h.: ■ { 1 h . (\j } — ^ijk^-ki {Gy , (Vj } — 0. (1-6) 

The functions T and L are Casimir functions of the bracket (1.6). Therefore, on each level of the area 
integral fixed in P 5 , 

Pf = { L = l}cP\ (1.7) 

this bracket induces a symplectic structure, whereas the restriction of the dynamical system of the 
Euler-Poisson equations is a Hamiltonian system with two degrees of freedom. This restriction corre¬ 
sponds to the reduction of order in the Routh sense in the system with three degrees of freedom on the 
phase space TSO[2>). For the complete integrability, we need one more general integral independent 
with H and L almost everywhere on P~‘. (Here the independence of functions at a point means the 
linear independence of their differentials.) The independence of two functions almost everywhere is 
sometimes called the functional independence [ 8 ]. Cases where one can indicate submanifolds in P 5 
of positive codimension that are distinct from Pf and are invariant for the system (1.1), on which this 
system is integrable in quadratures, are called particular cases of integrability. 

The following assumptions are called the Kovalevskaya-Yehia conditions. As a movable coordinate 
system, we take the orthonormal trihedron Oeie 2 e 3 of the principal axes of the tensor of inertia. 
Assume that the corresponding principal moments of inertia are in ratio 2 : 2 : 1 , the center of mass 
lies in the equatorial plane c • e 3 = 0 , and the gyrostatic moment is directed along the axis of 
dynamical symmetry: A = Ae 3 . Then, using the remaining freedom of the choice of movable axes and 
appropriately taking measurement units, we can achieve that 

c = {1,0,0}, I = diag{ 2 , 2 , 1 }, A > 0; (1.8) 

then Eqs. (1.1) take the form 

2 chi = ^ 2(^3 — A), 2u>2 = —Wl (<^3 — A) — G 3 , UJ 3 = G 2 , 

. . . 

d\ — CK 2 CC /3 — 0136^25 Oi 2 — C^3^1 — <^1^3} — C^ 1^2 — C^ 2 ^ 1 * 

This system possesses an additional Kovalevskaya-Yehia integral, which we denote by K. Thus, there 
exists a complete involutive set of integrals 


H — U)\ + + —UJ 3 — Gl, 


L = uqai + u; 2 a 2 + ^(^3 + A)a 3 , 


I\ — (u)\ — u>2 + G1) 2 + { 2 (jj\UJ 2 + G2) 2 + 2A [((U3 — A)(a; 2 + 0J2) + 2CU1G3] . 


( 1 . 10 ) 
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We note another important aspect of terminology, which sometimes leads to a false sense of the 
presence of more general problems than the problem described by the last formulas. Since the center 
of mass lies in the equatorial plane of the ellipsoid of inertia, each orthonormal basis in this plane is a 
basis of principal axes of inertia. Taking into account the arbitrariness of the choice of measurement 
units, we can write conditions (1.8) in the form 

c = {ci,C2,0}, I = diag{2C, 2C, C}, Ael. 

In the literature, systems determined by these conditions are called “Kovalevskaya-Yehia-type sys¬ 
tems.” We emphasize that each of these systems is completely equivalent to the system with condi¬ 
tions (1.8), i.e., can be obtained from it by a global nondegenerate linear change of variables. Therefore, 
in our opinion, the use if this term is incorrect. 


1.2. Exact solutions. Before the paper [69], the following three general cases of integrability for 
A / 0 were known: the generalization of the Euler case proposed by N. E. Zhukovsky (c = 0), the 
generalization of the Lagrange case (A x c = 0), and the generalization of the Goryachev-Chaplygin 
proposed by L. N. Sretensky. A series of particular cases of integrability of a heavy gyrostat (c / 0) 
was indicated in the papers of P. V. Kharlamov and E. I. Kharlamova. As was noted above, two of 
them are fundamentally important for the study of the system (1.9). Let us dwell on them in detail. 

1.2.1. Solution of P. V. Kharlamov. The first exact solution 1 was published in [36]. Assume that in 
system (1.1), (1.2) 

u> = ojq + r(t)e, r / const, (1-11) 

where u>o and e are constant in the body. Without loss of generality, we set |e| = 1 and cjo • e = 0. 
The representation (1.11) is a most general for solutions that admit two particular integrals that are 
linear with respect to the angular velocity. It is well known that all solutions of the system (1.1) are 
defined and bounded for all t since all levels of the energy integral are compact. This implies that if 
a function r(t) satisfies an equation of the form 

k 

br = ^ air 1 (1.12) 

i=0 

with constant coefficients b and a*, then all these coefficients are equal to zero. 


Theorem 1 (P. V. Kharlamov [36]). For each motion of the form (1.11), the following conditions 
hold: 

(1) the axis Oe is a principal axis of inertia at the immovable point ; 

(2) the center of mass lies in the plane passing through the immovable point perpendicularly Oe, i.e., 

c • e = 0; 

(3) the total kinetic moment of the system is coplanar to the vectors c and e. 


Proof. We multiply the first vector equation (1.1) (the Euler equation) scalarly by a nonzero constant 
vector c. Then we obtain an equation of the form (1.12) in which k = 2 and b = c • Ie, 02 = c • (Ie x e). 
Thus, 

c • Ie = 0, c • (Ie x e) = 0. (1.13) 

We write the following geometric identity: 


1 


a = — [(c • a)c + (c x a) x c]. 


(1.14) 


We denote by h the constant of the integral (1.4); then from (1.1) and (1.4) we obtain the expressions 

c X a = K - (la; + A) x w, c ■ a = h — -la; • a ;. 


1 An exact solution is a general or particular case of integrability explicitly reduced to quadratures. 
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Substitute them on the right-hand side of (1.14); then the relation L = t and the area integral (1.5) 
yield an equation of the form (1.12) in which k = 3 and 

b = (Iu> + A) • (Ie x c), a 3 = 2(c • e)(Ie • Ie) — (c • Ie)(Ie x e). 

Therefore, taking into account (1.13), we have 

c • e = 0, M • (Ie x c) = 0. (1.15) 

For nonzero c, the first of these relations together with (1.13) imply Ie x e = 0; therefore, the axis Oe 
is a principal axis of the tensor I. Then the relation (1.15) express the remaining assertions of the 
theorem. □ 

Note that we do not impose any requirements for the vector uj o, which is constant in the body. In 
particular, for ujq = 0 we see that the theorem is valid for all pendulum motions of a gyrostat in the 
gravitational field. 

In accordance with the above, we introduce a movable coordinate system so that e 3 = e and the 
following notation for the components of the constant vectors and the tensor of inertia: 

u 0 =pei+qe 2 , A = Aiei + A 2 e 2 + A 3 e 3 , I = diag {A,B,C}. 

Let Mo be the projection of the vector M on the plane Oeie 2 . By Theorem 1, 

Mo = me, m = const. (1.16) 

The projection of the Euler equation on the plane Oeie 2 takes the form 

r (me — Cuo) + cc 3 c — A 3 a^o = 0 . 

Since, by the assumption, c / 0 and r / const, we see that 


u 0 = £c, (1.17) 

(m — eC)r + cc 3 — A 3 e = 0. (1.18) 

From (1.16) and (1.17) we obtain Mo x cjo = 0; therefore, the projection of the Euler equation on the 
axis Oe 3 has the form 

Cf + c 2 ai — cia 2 = 0. (1.19) 

At the same time, differentiating (1.18) and taking into account (1.17), we obtain 

(to — eC)f + e(c 2 ai — ci<a 2 ) = 0 . ( 1 . 20 ) 

The compatibility condition for (1.19) and (1.20) yields 


to = 2 Ce. ( 1 . 21 ) 

Finally, taking arbitrarily the constants A, B, C , ci, c 2 , A 3 , e, and h, we obtain a solution in the form 

Wl=£Cl, w 2 = ec 2 , u 3 = r, 


1 

~2 i 2 

4 + 4 L 


^2 2 i 2 

cf+4 L 


t;CV 2 + K) ci - C\J R(r)c2 


Cr 2 + /i* ) c 2 + CifW)ca 


( 1 . 22 ) 


a 3 = e(A 3 - Cr), 

where the following notation is used: 

h* = £ — (Aci + Bel) - h, R(r) = ^ j (c 2 + cf) [l - £ 2 (A 3 - Cr) 2 } - Q <> 2 + h^j j . 

The two components of the gyrostatic moment, by (1.16) and (1.21), are defined by the relations 

Ar = (2(7 - A)eci, A 2 = (2 C - B)ec 2 , (1.23) 
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and the evolution of the variable r, which has been chosen as the principal variable, is described by 
the equation 

r = yj R(r). (1-24) 

Due to the arbitrariness of the choice of the measurement units, only three physical parameters among 
A, B, C, ci, C 2 , and A 3 are substantial. The parameters e and h related to the choice of an initial 
point are arbitrary. Therefore, formulas (1.22) and (1.24), for fixed physical parameters, describe a 
three-dimensional invariant submanifold of the phase space stratified into periodic trajectories (with 
possible bifurcations). This solution generalizes the classical Bobylev-Steklov case. 

Now we consider a gyrostat under the Kovalevskaya conditions (i.e., A = B = 2 C and the center of 
mass lies in the equatorial plane) and write the solutions (1.22), where the axis of dynamical symmetry 
is the axis Oe (which is a principal axis by Theorem 1). From (1.23) we see that Ai = A 2 = 0. 
An arbitrary axis lying in the equatorial plane is principal; therefore, without loss of generality, we 
assume that C 2 = 0 and hence 0 J 2 = 0. We choose measurement units so that C = 1 and ci = 1; as an 
independent parameter, we take, instead of e, the first component of the angular velocity p. Omitting 
the index of the sole nonzero component A3, from ( 1 . 22 ) we obtain the parametric equations of the 
three-dimensional manifold, which will be denoted by Adi: 


Mi : 


LOl = p, UJ 2 = 0, 
ai = ^r 2 +p 2 - h, 


W 3 = r, 

a 2 = yj R( r ), 


a 3 = -p{r - A). 


(1.25) 


Here 

R = —-j-r 4 ~ (2p 2 — h)r 2 + 2A p 2 r + 1 — (p 2 — h ) 2 — p 2 A 2 . (1-26) 

The dynamics on M\ is governed by Eq. (1.24), which determines periodic solutions and their possible 
bifurcations if immovable points appear. The exact solution (1.24), (1.25), (1.26) generalizes the family 
of remarkable motions of the fourth Appelrot class of the classical Kovalevskaya top to the case of a 
gyrostat. Note that generalizations of the fourth Appelrot class to the cases of a top and a gyrostat in 
double potential fields (including fields with gyroscopic terms, see [33]) were obtained in [41, 61, 62], 


1.2.2. Solution of P. V. Kharlamov and E. I. Kharlamova. Another exact solution for a gyrostat in 
the gravitational field was constructed in [37, 51]. The initial assumptions are as follows: first, the 
center of mass and the gyrostatic moment belong to one of the principal planes of the tensor of inertia 
and second, there exist particular algebraic integrals that guarantee the possibility of expression of all 
phase variables through one additional variable. Under the conditions (1.8), the most general form 
of the corresponding solutions was presented in [37], where they was united into a single family. We 
write it in an appropriate notation, as was proposed in [48]. We fix the constant t of the area integral; 
let s be a nonzero constant. We set 


x* = e + aV, p 2 = 1 — 


2 M 


g 2 = - 


K 


+ 


i A 2 

pY + — ) -1 

sx J 


(X,Y) = 


{ (cos cr, sin a), p 2 > 0 , 
(cosher, i sinher), p 2 < 0 . 


(1.27) 


Here i is the imaginary unit and a is an auxiliary variable. The invariant manifolds .M 2 and M 3 are 
defined, respectively, for s < 0 and s > 0 by the following system of parametric equations: 


M 2,3 : { 


oj\ = -xpK, 


W2 = - py/sQ , 


A sX + tpY n 2^2 0 v- m 

a 1 =- Zk Y , 0 L 2 = —ZxY y/sy, 03 = 


CU3 — A + 2 A, 
tX - A spY 


(1.28) 


X 


X 
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The dynamics induced by the system (1.9) is governed by the equation 

& 2 = sgn(p 2 ) s Q 2 . (1-29) 

Since £ and s are free parameters, for given physical parameters of the problem we have a three- 
dimensional invariant submanifold in the phase space stratified into periodic trajectories (as in the 
first solutions, bifurcations of periodic trajectories can occur on this manifold). 

We can immediately verify that for A = 0, on the family of trajectories (1.28), (1.29), there exists 
the following relation between the integral constants: 

(h - 2£ 2 ) 2 - k = 0 , 

which corresponds to remarkable motions of the 2 nd (s < 0) and 3rd (s > 0) Appelrot classes and 
the constant s is a multiple root of the polynomial in the Kovalevskaya differential equations. There 
exist corresponding generalizations for double fields (see [41, 61, 62]), but in the case A / 0 they are 
not reduced to quadratures. 


1.2.3. Solution of I. N. Gashenenko. Attempts at an explicit separation of variables in the reduced 
system on P 4 for the Kovalevskaya-Yehia gyrostat were unsuccessful even in the case £ = 0. In [9, 52, 
65], changes of variables on Pq that lead to equations of AbeKJacobi type with complex variables were 
proposed. Attempts at the reification of these equations falied, so they are not suitable for the study 
of the system (in particular, for the search for conditions of the existence of real solutions; C. Jacobi 
supposed that this is the main goal of the transition to equations with radicals; see [53]). 

I. N. Gashenenko [11, 12] proposed the reduction to quadratures of solutions of the system (1.9) 
on a four-dimensional submanifold in P 5 distinct from the level of the area integral. The family of 
solutions obtained is completely classified by the type of motions (periodic, asymptotic to periodic, or 
two-periodic). We recall the main results of the papers [11, 12]. Consider the relations between the 
constants h, £, and k of the first integrals (1-10) generated by the family of solutions (1.25): 

£ = — ^p{ 2 h — A 2 ) +p 3 , k = 1 — p 2 ( 2 h — A 2 ) + 3p 4 . (1.30) 

The complete preimage of the set of such values ( p and h are arbitrary) in P 5 is the stratified manifold 
G 4 = IV 4 U.Mi, where dim N 4 = 4 and dN 4 = M\. We pass on N 4 from (u, a) to the new coordinates 
{x,y,z,a,/3, 7 ) by setting (see [11]) 

LOl = p — xM~ 4 , U 2 = —yM ~ 4 , 

w 3 = 2 z + A-MxyM -1 , M = x 2 +y 2 ^ 0 , 

a 1 = —2 a + 4(.x 2 — y 2 )^ 2 M~ 2 + 2px(l — 2 px)M~ l + 4A xyM -1 , (1-31) 

012 = —2(3 + 8 xy 7 2 M ~ 2 + 2py(l — 2 px)M ~ 4 + 4A yyM -1 , 

0:3 = 2(1 — 2px)'yM~ 1 — 2 pz. 


By (1.30), from the set of equations of motion and first integrals, one can extract the following closed 
subsystem: 


z = Vf( Z ), +\(p 2 ~ h )r P = Vf(z), 




+ W~ h) 


(1.32) 


its solution z(i), aft), (3{t) can be expressed through elliptic functions of time. Introducing the notation 

(1.33) 


Ll = ±-(2h-2p 2 -\ 2 ), L 2 = -r{2h — 6 p 2 — A 2 ), L 3 = 4(p 2 + X 2 )L, - \ 
lb 4 4 
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for independent parameters, the function <J>, and the auxiliary variable y by the formulas 

$(*) = (z ~ A ) 2 - L 2 , y = y^~ 1 / 2 ( z ), 

we obtain an elliptic quadrature for y and finite expressions for other variables: 


(1.34) 


/ 


dr / 


Li — L 2 L 3 y 2 


+ 


dz 




= const, 


y = y^ 1 / 2 (z), x = - QpL 2 1 + y/3$ 1/2 (z) + y$ 1/2 (z)(z - X)L 2 ^ , 
7 = ^PzL^ 1 + \y/3<f>~ 12 (z) + y<S> 1/2 (z)(a + p 2 )^ 1 . 


(1.35) 


the solution 

obtained 

is 

performed with 

following cases (: 

see 

[11]): 




(I) 

L\ 

> 

0, 

L 2 

> 

0, 

l 3 

> 

0; 

(II) 

Li 

> 

0, 

L 2 

< 

0, 

l 3 

< 

0; 

(III) 

L] 

> 

0, 

L 2 

< 

0, 

l 3 

> 

0; 

(IV) 

L] 

> 

0, 

L 2 

> 

0, 

l 3 

< 

0; 

(V) 

Li 

< 

0; 







(VI) 

L\ 

= 

0; 







(VII) 

L 2 

= 

0; 







(VIII) 

l 3 

= 

0. 








(1.36) 


As was noted in [11], in the cases I and II, two-periodic motion are possible; in the cases III, IV, 
VII, and VIII, motions distinct from (1.25) asymptotically tend to trajectories on Adi; in the cases V 
and VI, motions distinct from (1.25) are impossible. Below, we discuss these classes and find the 
relationship between them and types of singular points of the moment mapping. 

2. Critical Set of Moment Mapping 

2.1. Lax representation. Although in [69] the Kovalevskaya integral was generalized to the case 
of a double field, this problem was not treated as an integrable problem since double fields prevents the 
existence of the area integral. However, later in [24], Reiman and Senrenov-Tyan-Shansky obtained 
a Lax representation with a spectral parameter for a gyrostat of the Kovalevskaya type in a double 
Held, which implied the existence of another first integral, which turns into the square of the area 
integral if the second field disappears. We present the corresponding Lax representation and resulting 
consequences under the Kovalevskay-Yehia conditions. 

Following S. V. Kovalevskaya, we introduce the complex variables (i 2 = —1): 

X\ = Qq + i CC 2 , X2 = cr\ — i Ct 2 , 

W\ = coi + iw 2 , u >2 = Ufi - iw 2 , (2.1) 

z = a 3 , w 3 =u 3 

(the last row is accepted for convenience and uniformity). We denote the differentiation d/d(it) by 
the prime and write the system (1.9) as follows: 

2w[ = — (w\W3 + z), 2w' 2 = W2W 3 + z, 2iu' 3 = x 1 —x 2 , 

x\ = — X\W 3 + ZWl, x ' 2 = X2W 3 — ZW 2, 2 z' = X\W 2 — X 2 W\. 

The Lax representation for this system, as a particular case of results obtained in [24], has the form 

B' = BA- AB, (2.3) 


( 2 . 2 ) 
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where 


B = 


2A 

X2 

— 2w\ 

z 


w 3 

0 

UJ2 

0 


K 

K 


2 


2 


Xl 

—2A 

z 

2w\ 


0 

w 3 

0 

W\ 

K 


K 

A — 


2 


2 

— 2w± 

z 

-2w 3 

— — — 4k 

) ^ — 

W\ 

0 

W 2 

K 


K 

K 


2 


2 


z 

2u> 2 

X 2 

- h 4 k 

2 W 3 


0 

W 2 

— K 

W 3 

K 

K 



2 


2 


(2.4) 


k is the spectral parameter, and the derivative in (2.3) is taken with respect to (2.2). The equation 
for eigenvalues fj, of the matrix B defines an algebraic curve associated with this representation [55]. 
We set s = 2k 2 . The equation of the algebraic curve becomes 

T 


/-4 


- (2h + A ) + 2s 


/i 2 + 4 


\ + -(M 2 - 2 h- A 2 ) + (4 k + 8 A 2 h) - 8 A 2 s 
■s z s 


= 0. 


(2.5) 


We introduce the integral mapping (the moment mapping) generated by the integrals (1.10): 

J = LxHxK : P 5 -> M 3 . (2.6) 


The set X of its critical values is called the bifurcation diagram; it is a classifying set in the study of the 
phase topology of the system. Research shows, in the presence of a Lax representation, the bifurcation 
diagram of the mapping J is contained in the set of values (£, h, k) for which the curve (2.5) ceases to 
be nonsingular, i.e., either becomes reducible (the left-hand side of Eq. (2.5) splits into the product 
of rational expression) or has a singular point in the standard sense. Thus, we can predict the result 
of the following assertion. However, for a rigorous proof, one needs to perform direct calculations on 
the set of critical points of J. This set will be stated below. 


Theorem 2. The bifurcation diagram of the integral mapping LxHxK is contained in the union of 
the following ( intersecting ) surfaces in M 3 (£, h, k): 


IIi : 



1 ^ 2,3 : 



(h~~^ - s 2 , s € M; 

_ 2A 2 y-2s) -A 4 + ^, 

\{ h + Y)~ x 2 s 2 ~Ts' s€M \W- 


Here s < 0 for H 2 and s > 0 for n 3 . 


(2.7) 


( 2 . 8 ) 


Remark 1. In the form presented above, the equations of bifurcation diagram have been obtained 
from (2.5). As independent parameters on the surfaces, the parameters h and s were chosen. Since in 
this case h and l 2 are linearly dependent, we can express k and h through s and i. There is no problems 
on the surfaces n 2 and ^ 3 . On ni, due to the relation s = 0, an artificial singularity appears. From 
the geometric point of view, this corresponds to the self-intersection line of the surface Hi, which is not 
a singularity in the parametric representation (2.7). However, for the study of bifurcation diagrams of 
reduced systems on Pf, we need their expressions through s and £, which can be written in the form 

f f 2 A 2 f 4 1 

ni = |h = ^ + Y + S , = —+ i, 

Lljfc=l, ^ = 0 }u|a: = 1+ , £ = 0 j . (2.9) 
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2.2. Description of the critical set. In [25], the set of critical points of the moment mapping was 
first described and parametric equations of singular surfaces 11, were also obtained. A more detailed 
presentation is contained in [26]. In [13], similar results were independently obtained; it was also 
proved that the surfaces II, naturally appear as discriminant sets of some polynomials. The equations 
of the surface III follow from the results of [ 11 ], but in this paper, the corresponding conditions for 
the constants of integrals were not related to the notion of critical points. 

The set C of critical points of the moment mapping is stratified by the rank of this mapping. Since 
the integral L is regular everywhere and stratifies P r> into smooth symplectic, leaves (1.7), we naturally 
use the following terms. 

Definition 1. The rank of a point x € Pf C P 5 is the rank at this point of the restriction mapping 

J t = HxK |p 4 : Pf -> M 2 . (2.10) 

Hence, the rank of a point is one less than the rank at this point of the mapping (2.6). Therefore, 
we have 

C = C°UC\ C = {x € Pf | rankJi(x) = *}. ( 2 . 11 ) 

The bifurcation diagram defined above is the J-image of the set C. 

The presence of two singular surface (2.7) and (2.8) generates another partition of the critical set. 

Theorem 3 ( [13, 25, 26]). The set of critical points of the moment mapping has the form 

C = M\ U M -2 U M 3 , (2.12) 

where M\, M 2 , and M 3 are defined by Eqs. (1.25) and (1.28). In this case, each set Mj is the 
closure of the corresponding set J _ 1 (n,) n C 1 (j = 1,2,3). 

Equality ( 2 . 12 ) follows from the analytic description of the set C obtained in [25, 26] by examining 
the minors of the Jacobi matrix 

d(L,H,K,T) 
d (u>,cx) 

An elegant proof based on the reduction of integral equations to a symmetric system of six complex 
coordinates on M 6 (lu, a) can be found in [13]. The fact that the image of each of the manifolds Mj is 
contained in the corresponding surface n, can be immediately verified. For M 1 , we must take 

s = h~Y~P 2 ( 2 - 13 ) 

in the formulas (1.30). 

Definition 2. The manifolds Mi, M 2 , and M 3 with induced dynamics are called the first, second, 
and third critical subsystems in the Kovalevskaya-Yehia case, respectively. 

Remark 2. We use the term “manifold” for the set Mj meaning their invariant properties (“invariant 
manifolds”). In fact, it may occur that they are smooth manifolds only almost everywhere. 


The representation of C in the form of the set of solutions of the system of invariant relations can 
be obtained as a particular case of the corresponding result from [41]. 

Theorem 4. The set of critical points of the integral mapping J consists of the following subsets 
in P 5 : 

(1) the set determined by the system 

Ri = 0 , R 2 = 0 , (2.14) 

where 

R\ = UJ2, R -2 = (<^3 — A)uq + 03 ; 
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( 2 ) the set determined by the system 
where 


F\ =0, F 2 = 0, 


F\ = ^012^2 + [lJi + 6^2 + C^l) J |^(^1 + OL 2)^3 — ^0^3 ( CX\U)i + O' 2^2) 

+ ^ — 2CX3UJ1 (to\ F C0 2 ) — 0L20L^,U} 2 (3a3 + 2W1CJ3) 

+ Ql — SttgCJl + 202^2 ( W 1 + W2) + 03 ( — Wf + CO 2)^3 

+ 02 a 3 w 3 ~ w i ( w i + <^2 — co |) a 2 wf + a3a;3 + wi (w 2 + w i) 

■P ^ — CK3 — 2 a|wiw 3 T (cr 2 — ^2)^1 T 2aia2W2 ^3 “I - crier 3 (ca 2 T to 2 F cu 2 ) ^ A^ 
F aia 3 co 3 X 3 , 


F 2 = 


F < a 2 


a 2 coi F <o 2 (cof F <o 2 — ai) (a ; 2 + a 2 ^co 3 — 2a^(^aiioi F oc 2 to 2 ) 


co\( — a\ F 2ai ( to\ F co 2 )) + a 3 (cj 2 — wl ) w 3 + co 2 co 2 (co'{ + co\ F W 3 ) 


— W2 


2a|(w 2 + w 2 ) -aices(a 3 - 2 wiu; 3 ) + a\(coj F to 2 — wf) 


A 


+ |a|w 2 W 3 - (af + 2a|)a;2W3 + «2 2aiwiw 3 + a 3 (w 2 F to\ F a|) 

+ a 2 Q!3W3A 3 . 


A 2 


(2.15) 


We can prove that the systems (2.14) and (2.15) are invariant by differentiation with respect to the 
system (1.9). We immediately verify that the points (1.25) satisfy the system (2.14) and the points 
(1.28) satisfy the system (2.15); taking into account the dimensions of the corresponding subsets, we 
obtain the relation ( 2 . 12 ). 

From the results of [41] we obtain, as a particular case, that the systems (2.14) and (2.15) describe 
the set of critical points of the restriction to P 5 of the function with undefined Lagrange multipliers 

2L 2 F{t -l)H F sK. (2.16) 


The differential of this function is preserved by phase flows; therefore, s and r expressed through the 
phase variables become particular integrals of the critical subsystems A 4j: 


Mi : 


M2,3 : 


co\ 

co 3 - 2 la 3 
2(ca 2 F CO 2 F Ac^ 3 )A 


r = 1 + 2 lto\, 
t = 2A 2 s. 


(2.17) 

(2.18) 


These values of s determine the corresponding points of singular surfaces IF, in (2.7) and (2.8). We 
note that, as (2.17) shows, on the set M\ the function (2.16) can be “cancelled” by L without loss 
of critical points. On the set .M 2,3 this cannot be done since there exists a set of points of global 
dependence of the functions K and H on P 5 : on this set L vanishes and hence the undefined multiplier 
of dL in the nontrivial zero combination of differentials is proportional to L. 

We also note that the results of [41] imply that the Poisson brackets of pairs of functions from (2.14) 
and (2.15) have the form 


{R\,R 2 } 



(2.19) 
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{Fi,F 2 } = ^|(l - 8A 2 s 3 ) y 2s 2 “ 2 (/» + y) « + 1, (2-20) 

where 

-Fg = (wj + ca 2 + Xuj^j (aicui + Q2 w 2 + Aa^) + (a.2UJ\ — (I 1 W 2 )" • (2.21) 

Vanishing of the brackets (2.19) and (2.20) corresponds to the cases of degeneration of the symplectic 
structure induced on the two-dimensional manifolds which are the phase spaces of Hamilton¬ 

ian systems with one degree of freedom. Below, we clarify the relation of this phenomenon with types 
of critical points. Looking ahead, without rigorous proofs, we say that, since nondegenerate critical 
points are organized into symplectic submanifolds (see the corresponding assertion in [8]), points on 
the critical set at which the induced symplectic structure is degenerate must be degenerate as critical 
points of the moment mapping. 


3. Relative Equilibria Are Critical Points of Rank 0 

3.1. Dependence of the energy and area integrals. Immovable points of the Euler Poisson 
equations are the projections on P 5 of motions of the body in which the trajectory in 50(3) coincides 
with an orbit of the symmetry group, i.e., rotations about the vertical with a constant angular velocity. 
In Routh’s theory of reduction of order such immovable points are called relative equilibria. In Smale’s 
theory (see [68]), these points are critical points of the energy-moment mapping. In the dynamics of 
a rigid body, this means that the energy integral and the are integral are dependent. Bifurcation 
diagrams of the mapping 

LxH : P 5 -> M 2 (3.1) 

are called Srnale diagrams. 

In its turn, on a symplectic leaf (1.7), critical points of the mapping (3.1) are critical points of the 
“reduced” Hamiltonian 

H( = H\ P i : Pf -> M. (3.2) 

If we want to emphasize the dependence on A, we write Pf( A) and The results of [30] imply that 
almost all such points are nondegenerate in the sense of Morse (below we consider a classification of 
these points in more detail). However, a nondegenerate critical point of the Hamiltonian is critical 
for all other first integrals. Therefore, the set considered coincides with C°. Thus, the problem of 
the study of critical points of the mapping (3.1) coincides with the problem of the study of the set of 
rank-0 critical points of the mapping (2.6) in the sense of Definition 1. 

In the sequel, we will use the following assertion for the search for critical points of functions on P 5 
without introducing an additional undefined multiplier corresponding to the constraint (1.3). 


Lemma 1. Let f be a smooth function on M b (u;, a). The set of critical points of the restriction of f 
on the submanifold P 5 is determined by the system of equations 

df _ df _ 

du ‘ ^ da ' (3 " 3) 

In the complex variables (2.1) this system has the form 

d^ = d^ = d^ = Q 

dwi dw 2 dw 3 ’ , . 


9 df df 

2Z T - 

0X2 dz 


x 17- X 2W— 

OX 2 OX 1 


The proof is obvious. Note that the last three equations in both systems (3.3) and (3.4) are linearly 
dependent. The choice of an independent pair in the case of general position is determined by reasons 
of convenience. 
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We write Eq. (3.3) for the function fn = H — 2 aL. We have 

u\ = aai, U2 = cr CH2, cu 3 = ( 7 a 3 , 

«3 - [«i(w3 + A) — 2a S 0Ji]a = 0, a 2 - 2(oqu; 2 - o 2 wi)cr = 0, (3.5) 

[2a 3 cu 2 - a 2 (w 3 + A)]cr = 0. 

This rank-5 system together with Eq. (1.3) allows one to express the phase variables and the undefined 
multiplier a through one variable chosen as an independent parameter. In [29], the parametrization 
of the set C° by the variable a\ was used. 


Proposition 1 (P. E. Ryabov [29]). The set of rank-0 critical points is described by the following 
system of equations: 


Wl = 2 Tf x - <l2) ’ 

oq = x, 


U 2 = 0 , 

0L2 = 0 , 


w 3 — - d 2 ), 

a.3 = d\, 


where 


d 2 = ±\j A 2 - -(1 - x 2 ), 


d\ = ±v1 — x 

the signs of d\ and d 2 are arbitrary, and the parameter x satisfies the conditions 

x Jl-1.0)U[CA,1]. 4 > 0. 

\(-l,0)U[c A ,l), d 2 < 0, 8 V ) 

In this case, the values of the first integrals and the undefined multiplier a are as follows: 

A(1 — x 2 ) — d 2 ( 1 + 3x 2 ) 


1 

d~i L 


(3 — ar)A — (1 + ar)d 2 


h = 


2 x(X + d 2 ) 


,. = ^(A-d 2 ). 


(3.6) 


(3.7) 


Note that the system (3.5) has no solutions in P 5 with aq = 0 and hence the singularity x = 0 is 
nonremovable. 

Equation (3.7) is a parametric equation of the Srnale diagram. In the study of bifurcation diagrams 
of the mapping (2.10), the first equation in (3.7) allows one to find all rank-0 critical points containing 
in a given level Pf. 

Another parametrization of the set C° was proposed by Gashenenko (see [13]); he noted that such 
points belong to the intersection of with the union of the set M 2 an d M 3 . Therefore, it is 
convenient to take cj 3 = r as an independent parameter in the relations (1.25). In this case, the 
polynomial (1.26) must have a multiple root, which allows one to express all unknowns through r. 
Surely, the same expressions follow from the system (3.5). Taking into account the results of [13, 16], 
we arrive at the following assertion. 


Proposition 2 (I. N. Gashenenko [13]). The set of rank-0 critical points is described by the following 
system of equations: 


oj 1 = A 


r — A 


oq = 


-r -I- -d 

r — A 


-r H- -d 

r — A 


iw 2 = 0, u 3 = r, 


«2 = 0, a 3 = =F(r- A) 1 


—r 


r — A 


where 

d = ±y/r 2 {r- A) 2 + 4, 

the parameter r runs over the set 

r € (—00, 0] U [0, A) U (A, +00), 


(3.8) 


(3.9) 

(3.10) 
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and the sign of d coincides with the sign of r{r — A) for r/ 0 and is arbitrary for r = 0. The values 
of the first integrals are as follows: 


1 


= - \) + d\J- 


-r + 


1 


r — A 


, 1 / n 2r — A 


(3.11) 


and the undefined multiplier is equal to 

cr = =F 


Hr- A) + d[ 

T/ie signs of u i, 03 , and a are matched (all upper or all lower). 

Comparing with Eqs. (1.25), we see that rank-0 points appear in the subsystem M\ if we set 

1 


P =2 


-r + 


r — A 


(3.12) 


Remark 3. We verify that at all point found above the differential d(K\ P 4.) also vanishes. For this, 
we write Eq. (3.3) for the function 


fr< = K - 4 y L, 7 = =F 


r — A 


-r + 


r — A 


For Smale diagrams, the value of the integral K is irrelevant, but we need this value at the points (3.8) 
below. By [13], we can represent it in the form 

A 


k = 


r(r — A) — d r(r — A)(4r — 3A) — A d 


(3.13) 


4(r — A ) 2 

We also calculate the values of the particular integrals s in the systems Aij at the points (3.8), i.e., 
the values of the parameter s on the surfaces II j at rank-0 point. Using (2.17) and (2.18), we find 


Mi : s = ^[A(r - A) + d], 
M 2 ,3 - s = —^-\ r {r ~ A) - d]. 


(3.14) 

(3.15) 


From (3.8) and (3.10) we immediately obtain that the set C° has exactly four connected components 
that are homeomorphic to M. In correspondence with the domain of r, we introduce the following 
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notation for subsets in C° determined by the formulas (3.8): 


'0 . 
1 • 

r £ [0, A), 

d < 0 , 

lim d = — 2 

r— »+0 

'0 . 

'2 • 

r £ (-oo, 0 ], 

d > 0 , 

lim d = 2 , 
r^—0 

'0 . 

3 • 

r £ (A, + 00 ), 

d > 0 . 



(3.16) 


The first two sets are connected, whereas the last set consists of two components that differ by the 
sign of uq. In C® and C!J, each value of r / 0 corresponds to exactly two points, and in C® the value 
r = 0 corresponds to the point u> = 0 , a = { 1 , 0 , 0 } of the lower equilibrium position of the body (for 
further convenience, we denote it, in correspondence with the fact that r —>• +0, by c+). In the set C 
the zero value of r leas us to the point u> = 0 , a = {— 1 , 0 , 0 } of the upper equilibrium position (we 
denote it by c_). On the sets (3.16), there exists an obvious symmetry 


symm : (uq, a 3 ) (-uq, -a 3 ), (3.17) 

which changes the sign of the area constant £, preseves the connected sets C® and Cf, and in the set C 3 
swaps the connected components. The structure of the family of sets C°(A) is shown in Fig. 1. 


3.2. Classification of rank 0 critical points. Here we present results related to an analytic 
classification of types of rank-0 critical points in correspondence with [44]. 

We recall some facts and definitions related to the notion of the type of a critical point of an 
integrable system (see [ 8 ]). 

Let M be a symplectic manifold. To each smooth function F : M -£ M, we put in correspondence 
a Hamiltonian vector held sgrad-F on M. We denote the Poisson bracket generated by the symplectic 
structure by {•, •}; then the differential equations of the system sgradF have the form 

x = {F, x}. 

We note that usually the arguments of the bracket in the right-hand side are written in the reverse 
order. It can be easily verified that just the expression {F, x} allows one to obtain Eqs. (1.1) from 
the definitions (1.6) and (1.4). There also exist other approaches to this inconsistency; for example, 
in [10], opposite signs in the definition of the bracket (1.6) are used, whereas in [23] the components 
of the kinetic moment with opposite sign are taken as momenta. 

Let dim M = 2 n and let a Hamiltonian system sgradF have n functionally independent, involutive 
first integrals T),..., F n , i.e., {Fj, Fj} = 0. A point C € M is called a critical point of rank k < n if 
the rank of the system of vectors sgradFj at the point C is equal to k. We introduce the notion of the 
type of a critical point following [ 8 ]. 

Consider a critical point C of rank n — m (m > 0). We assume that below in this section the 
subscript i runs over the set 1,..., to. By a linear change of variables with constant coefficients in the 
system of functions F\,..., F n . we can achieve that the point C be critical for each of the functions F\ 
and regular for all other points. Then sgradFj(C) = 0 and the linearization of this field at the point C 
is a symplectic operator a. t : T^M —> 7} M. The linear span 2l(£) of such operators is a subalgebra in 
the algebra of all symplectic operators on T^M. 

Definition 3 (see [ 8 ]). A point £ is called a nondegenerate critical point of rank n — m if 2l(£) is a 
Cartan subalgebra; this is equivalent to the following conditions: 

(i) symplectic operators Oj are linearly independent (dim 2 l(£) = to) ; 

(ii) there exists an operator a £ 21 (C) whose eigenvalues are distinct. 

We recall that eigenvalues of an symplectic operator can be split into groups: pairs of purely 
imaginary ±ia, pairs of real ± 6 , and quadruples of complex ± 6 ±ia (ab / 0). For an operator a,;, the 
projections of the field sgradF) on root subspaces of such groups have respectively a center, a saddle, 
or a focus. An operator a £ 21(C) with distinct eigenvalues is called a regular element. Choosing a 
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regular element at a nondegenerate point, we denote by mi, m 2 , and m 3 , the numbers of centers, 
saddles, and focuses, respectively (m = mi + 777,2 + 2777 . 3 ). These nonnegative numbers are independent 
of the choice of the regular element. 

Definition 4 (see [ 8 ]). A quadruple (n — m, mi, m 2 , m 3 ) i called the type of a nondegenerate critical 
point £. 

For systems with two degrees of freedom, the following more illustrative terms are usually used. 
For rank-0 points, the type (0,2,0,0) is called “center-center,” the type (0,1,1,0) is called “center- 
saddle,” and the type (0, 0,2,0) is called “saddle-saddle.” For rank-1 points, the type (1,1, 0,0) is 
called “center” and the type (1,0,1,0) is called “saddle.” In the problem considered, other types (in 
particular, focuses) do not occur. 

The type of a nondegenerate critical point largerly (but not completely) determines the Liouville 
fibration (i.e., the fibration of common levels of n involutive first integrals) in a neighborhood of this 
point. To construct a rough topological description of the Liouville fibration (for example, as a loop 
molecules without marks in a system with two degrees of freedom), it suffices additionally to know 
the number of connected components of regular and critical integral surfaces. For critical levels of low 
complexity (i.e., when each connected component contains few critical orbits), for example, for rank -0 
critical points with one point on a fiber, these data are sufficient for the search for a fine topological 
invariant (a marked loop molecule) since one has descriptions of all possibilities (see [ 8 ]). Below, we 
present such a description for the problem considered. 

Definition 5. A critical point £ is said to be strongly degenerate if the condition (i) of Definition 3 
is violated and weakly degenerate in the opposite case. 


We note that these terms are not conventional; we use them here for convenience. 

For the Euler-Poisson equations, the 2-form induced on P 5 by the symplectic structure of the 
manifold TSO( 3) is degenerate. The notions introduced above can be considered from the point 
of view of systems on Pf. The explicit transition to these systems makes calculations boundless. 
However, these calculations are not necessary. To the function F, the degenerate Poisson bracket 
( 1 . 6 ) assings the field sgradF determined by the equations 


• ,, OF dF 

M = Mx 8 M + “ x &C 


OL = OL X 


(3.18) 


Remark 4. If F is a Casimir function, i.e., it is identically dependent with L and T = a 2 , then the 
right-hand sides of Eqs. (3.18) identically vanish. In this connection, in the calculation of linearizations 
of the fields sgrad-F of arbitrary functions F and eigenvalues of the corresponding operators in M 6 , 
there is no need to consider undefined Lagrange multipliers for the functions L and T; we must only 
take into account Definition 1 and discard two zero eigenvalues, which necessarily exist. 


The physical model of a gyrostat is a system with four degrees of freedom (a body and a rotor) for 
which A is the constant of a cyclic integral. The Euler-Poisson equations are obtained by reduction 
of order in this system. Therefore, all reasonings on the Constance of any properties in the space of 
integral or other parameters must be naturally considered in the extended space of these parameters 
including the axis K = M(A). In this connection, we will use the following notation. Let A be a set 
and B( A) be a family of its subsets depending on a parameter A. We introduce the notation 

A (A) = A x R, A (B) = [J B( A) x {A} C A(A). (3.19) 

A 

Hence we consider the extended set A(C°) C A(P 5 ) = P 5 x 1. Recall that by (1.8) we assume that 
A > 0; the case A = 0 is considered only as a limit possibility where explicitly stated. By Proposition 2, 
the set A(C°) continuously twice covers the domain 

V° = {(r, A) £« 2 : A > 0, r ± A}. 
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Remark 5. In the sequel, we denote the images of the sets A(C?) in various spaces of parameters (con¬ 
stant general integrals and particular integrals, physical parameter A) by eh, i = 1,2,3. In particular, 
as subsets in T>°, they have the form 

c)i : {(r, A) : 0 < r < A, A > 0}, 

d 2 : {(r,A) : r < 0, A > 0}, (3.20) 

d 3 : {(r, A) : r > A, A > 0}. 

Subsets of these sets obtained as results of further detailed elaboration are marked by double indices. 

We recall the existence of a symmetry (3.17) on C°. 

Definition 6. We say that points ^ 1,^2 £ A(C°) belong to the same class if there exists a continuous 
path in A(C°) connecting these points (or connecting the point £1 with the point synnn^)) along 
which the type of critical points is preserved. 

Let (r, A) € V°. We denote by £±(r, A) the point (3.8) for r / 0 and, due to Proposition 2, for the 
chosen sign of (3.9), sgnd = sgn[r(r — A)]. By the notation adopted, 

lim. £±(r, A) = c + eCf, lim £±(r, A) = c_ € C 2 . 

1 —^+0 r—>•—0 

Definition 7. A point (r, A) € T>° is said to be separating if any its neighborhood contains images of 
points from A(C°) belonging to different classes. 

The ray of prohibited points r = A, A > 0 is a separating ray; we denote in by W. Since £1 € C® and 
£2 £ C 2 cannot belong to the same class, points of the form (0, A) are always separating. We denote 
the semiaxis r = 0, A > 0 by ttq. For r / 0, the types of critical points £±(r, A) always coincide, hence 
the point (r, A) is separating if and only if both critical points £±(r, A) are degenerate. 

Let an and ax be symplectic operators of the linearizations of the fields sgradR and sgradit at 
the point (3.8). 

Proposition 3. The set of strongly degenerate critical points £±(r, A) corresponds to the curve 

7 T 2 i : r = A — 0 < A < 1, (3.21) 

in the domain 5 2 . 


Proof. At strongly degenerate points, we consider the combination 

b = iq an + v 2 clk = 0. (3.22) 

Arranging the variables and, respectively, the elements of matrices in the order u\, u 2 , cu 3 , oq, a 2 , a 3 , 
we have 

bi 2 = ~(?* A) zq - 2 v 2 \(Q 2 + A + r) , (3.23) 

where we have introduced the notation 

Q = 

Obviously, Q / 0. From (3.22) and (3.23) we obtain 

zq = 2v 2 \(Q 2 T A + r), v 2 / 0. 

Substituting this in (3.22), we have 

b = 2 zz 2 (A — r)Q(Q 2 + A)<£ = 0, 
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where 



f 0 

0 

0 

0 

0 

0 \ 


0 

0 

0 

0 

0 

0 


0 


0 

0 

A + r 
Q(\-r) 

0 

£ = 

0 

—2 ry/r 

0 

0 

r( A + r) 
Q(\-r) 

0 


2 ry/r 

0 

-A Q 

r(A + r) 

0 

2 Xy/f 


Q(A-r) 

A — r 


0 

2rQ 

0 

0 

y/r(\ + r) 
A — r 

0 i 

a nonzero matrix. Therefore, Q 2 + A 

= 0, which is equivalent to the equation 


1 + A(r - A) 3 = 0 


(3.24) 


with the condition Q 2 < 0. But sgnQ 2 = sgn[(r — A )d\ coincides with sgnr (see Proposition 2). 
Therefore, from the solutions (3.24) we must take only points lying in 82 , which yields the curve (3.21). 


□ 


Lemma 2. Under the condition 


all critical points are nondegenerate. 


(2r — A)(r — A) + d = 0 


(3.25) 


Proof. Let the condition (3.25) hold. If we assume that d > 0, then r € (— 00 ,0] U (A,+ 00 ). But on 
this interval (2r — A)(r — A) > 0, which contradicts (3.25). Therefore, d < 0, r € [0, A), and 


1 , 


-Ml = -) (r-A) 

and finally r € [0, A/2]. By squaring, we obtain the equation 

(r — A) 3 (3r — A) — 4 = 0. (3.26) 

We represent its solutions on the required interval in the parametric form; for this, we introduce the 
variable x by the formula 


(see [48]). From (3.26) and (3.27) we find 


x = A — r 
x 4 + 4 


(3.27) 


d= — 


2 x 2 


and 


x 4 - 4 A 3.x 4 - 4 
r = - „ , A = 


2 x 3 ’ 


2 x 3 


The condition r € [0, A) is fulfilled if 


x > V2. 

The characteristic polynomial of the operator an takes the form 


(3.28) 

(3.29) 


Xh(u) = H + 


2 , \Jx 4 + 4 


2 x 


(see (3.28)); we see that it is useless in the problem on nondegeneracy. However, we calculate the 
characteristic polynomial of the operator ajc- 


Xk(u) = 


■ 2 (x 4 -4)2(4+ x 4 ) 

/i “T 


~14 


' 2 + ( x 4 + 4)(3x 8 - 7x 4 + 4) 2 


..14 
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Under the condition (3.29), all its roots are distinct; therefore, it is a required regular element of the 
algebra. □ 

Proposition 4. Assume that a point (r, A) € V () does not lie on the curve (3.21) and does not satisfy 
any of the equations 


r + A = 0, 

(3.30) 

(2r — A )(r — A) — d = 0, 

(3.31) 

(2 r — A )(r — A )r + Ad = 0. 

(3.32) 


Then the critical points f±(r. A) are nondegenerate. 

Proof. The characteristic polynomial of the operator an at the point £±(r, A) cancelled by p? in 
correspondence with Remark 4 has the form 


where 


1 


a = 


b = 


8 (r — A) L 

1 


Xh{t) = T 4 - 2a/U 2 + b, 


— (3 r — A)(2r — A)(r — A) + (r — 3A )d 


(3.33) 


(r — A) 3 (4r — A)r — 4A — (2r — A)(r — A ) 2 d 


8 (r — A) L 

The discriminant of the polynomial (3.33) 

« 2 " b = 6 l; + _ A i )2 [( 2r - A)(r - A) + d ] 2 

by Lemma 2 vanishes only if the condition (3.30) holds. The roots of xh with respect to p 2 can be 
explicitly expressed: 


2 1 
^ = "4 

= 


(2 r — A)(r — A) — d 

1 


2(r - A) 


(2r — A)(r — A)r + A d 


(3.34) 

(3.35) 


These values that determine the types of rank-0 points were first calculated in [30]. Thus, outside the 
set determined by the equations (3.30)-(3.32), all roots of xh are distinct and outside the curve (3.21) 
the algebra generated by the operators an and ax is two-dimensional. The proposition is proved. □ 


Proposition 5. On the curve 

vr 2 2 : r = -A, A > 0, 

i.e., under the condition (3.30), all critical points are degenerate. 


(3.36) 


Proof. We calculate the characteristic polynomial of the combination b = viaH+^ax at the points (3.8) 
under the condition (3.30). We have 


X(h)= h 2 + 


(Z 2 - 2)(mZ + u 2 ) 
2 Z 3 


, Z = \ 2 + VA 4 + 1. 


Therefore, the linear span of the operators an and ax does not contain regular elements. 


□ 


Note that if v\Z + zz 2 = 0, then all eigenvalues of the operator b are equal to zero, but the operator 
itself is nonzero, except for the value A = 1 /2 3 / 4 , which corresponds to the intersection point 7r 22 with 
the curve (3.21) of strong degeneration. 
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Proposition 6. The condition (3.31) is realized in the domain V° on the following curves: 


7 t 2 3 : r = 


^31 


r = 


x 4 — 4 
2x 3 
x 4 -4 
2 x 3 


A = 


A = 


3x 4 — 4 
2 x 3 

3x 4 — 4 
2x 3 


x € 


(y4/3,V2 


G f-v/473,0). 


All corresponding critical point are degenerate. 


(3.37) 

(3.38) 


Proof. From Eq. (3.31) we obtain Eq. (3.26), whose solutions can be represented in the form (3.28). 
If we assume that d < 0, then r € [0, A), and hence r € [0, A/3) and (2 r — A)(r — A) > 0, i.e., 
Eq. (3.31) does not have such solutions. If d > 0, then in the domains 62 and £3 we must substitute 
the corresponding range for x in (3.37) and (3.38). At such points, the characteristic polynomial of 
the combination v\an + 1^2 a k has the form 


. 4 (x 8 + 2x 4 - 8) [ - 2z^x 2 + V 2 (x 4 - 4)1 2 2 

X(M) = M 4 - 1 -^10- -T 2 


(3.39) 


and always has two zero roots. Therefore, the algebra of operators does not contain regular elements. 

□ 


Proposition 7. The condition (3.32) is satisfied in the domain V Q on the curve 

7T24 : r = - (a — \/A 2 + 4A 2 / 3 ^ , A > 0. (3.40) 

All corresponding critical points are degenerate. 


Proof. The equation 

r{r — A) = A 2 / 3 (3-41) 

is a consequence of Eq. (3.32) Obviously, it has exactly two solutions, one in each of the domains £3 
(r > A) and 82 (r < 0), where we must choose d > 0. But in the domain £3 we have (2r —A)(r —A)?’ > 0, 
so Eq. (3.32) does not hold. Therefore, only one root of Eq. (3.41) is a solution, namely, the lower root; 
this leads to the points of the curve (3.40). Setting Z = A 2 / 3 + y4 + A 4 / 3 , we write the characteristic 
polynomial of the combination v\an + ^2 &k a t these points in the form 


x{p) = p 4 + 


(Z 2 -8)(4 + Z 2 )|^2 Z 2 (Z 2 — 8) 2 — 64 -4iqZ 3 } 


512Z 7 


-T 


(3.42) 


It always has two zero roots. Therefore, the algebra of operators does not contain regular elements. 


□ 


The domain T>° with separating curves is shown in Fig. 2 (see also the notation for classes of 
nondegenerate rank-0 points. We see that the set A(C 3 ) contains a single class, which is denoted, 
as a subdomain in X>°, by A. It also includes the subset with the value r = 0; through this subset 
one can pass to symmetric points, so that this class consists of a single connected component. The 
set A(C3) contains two classes, which are denoted by £31 and 832 ', they do not contain points with 
r = 0 and hence each such class consists of two connected components. The set A(C 3 ) contains eight 
classes 821, ■■■ ,82s- Three of them, namely, 82 1, 82s, and 827, contain points with r = 0 and hence 
have one connected component; the other classes consist of two connected components. Figure 2 
also contains the curve £q, which is generated by the possibility specific for gyrostats, namely, by the 
presence of uniform rotations about the vertical axis (r 7^ 0), for which the area constant vanishes. 
This means that a pair of symmetric points lies on the same level of all first integrals (but not in 
the same connected component). If we cross the curve £ 0 , then the type of critical points does not 
change, but we show below that it causes a surgery of the Srnale diagram and changes the topology 
of the common level of first integrals and topological invariants. We see from (3.11) that the relation 
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I = 0 for r^O implies A(r — A) + d = 0. In correspondence with the definition of the sign of d, this is 
possible only on 82 ] therefore, the curve £q is determined by the conditions 

£ 0 : (r — A) 3 (r + A) + 4 = 0, r < 0. (3.43) 

We recall once more that, according to the agreement, A > 0. The curve £q divides the classes 827 and 
e>28 into subdomains, which we will mark with primes (see the figure). 

Since we have analyzed all cases where the values (3.34) and (3.35) vanish, the analysis of their 
signs in the generated subdomains of the (r, A)-plane leads to the following classification of rank-0 
points. 

Theorem 5 . In the extended phase plane A(P 5 ), the rank-0 critical points of the Kovalevskaya-Yehia 
case form four manifolds diffeomorphic to M 2 . Degenerate critical points correspond to the following 
five separating curves on the parameter plane: the curves tt 2 i~tt 2 A and 7r3i. 

According to the equivalence relation introduced above (the presence of a continuous path with a 
constant type of critical points that contains the first point with the second point or with the point 
symmetric to the second point with respect to the given inversion mapping), nondegenerate critical 
points can be divided into 11 classes. The four classes that contain one of equilibria are connected. 
The other classes that contain only the images of uniform rotations (relative equilibria) consists of two 
connected components that are mutually symmetric with respect to the inversion mapping. According 
to the notation of classes in Fig. 2, critical points have the following type (we list them in the order 
corresponding to the order of pairs of roots yf\ and 

(i) 821 , 822 , and 82 s are “saddle-saddle 

(ii) 823 and <531 are “saddle-center”-, 

(iii) 827 is “center-saddle”] 

(iv) 82 A, <$25; <$28; 8 \, and 832 are “center-center.” 
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Now we compare the classification of relative equilibria with respect to their types types and the 
classification of of relative equilibria in Gashenenko’s solution (1.31), (1.32), (1.35) with respect to the 
parameters (1.33). In Fig. 3, we show the signs of triples (Iq, L 2 , Lfi) and zero values of parameters 
on separating curves. It turns that the curve 7^4 does not affect the signs of triples, L 3 vanishes 
on the curve 7721 , but in the adjacent domains the signs of triples are the same. The parameters L 2 
and L 3 vanish on the curve 1:22 ■ Finally, we obtain the following correspondence between domains 
and classes (1.36): 


(I) 

^24) ^25 5 <^32 ; 

(V) 

r c■// e// 

^27>^28? 

(II) 

^27) <^28 i 

(VI) 

4>; 

(III) 

C>22, <^23, £ 31 ; 

(VII) 

7722, 7T23, 77 3 i 

(IV) 

<W 26 ; 

(VIII) 

7721,7722- 


We formulate as separate assertions an important property that was implicitly used by Ryabov and 
Gashenenko for classification of bifurcation diagrams and rough Fomenko invariants (we discuss these 
problems in more detail below). Proposition 8 shows that all rank-0 singularities have complexity 1 
(the notion of the complexity was introduced in [ 8 ]). Proposition 9 shows that, moreover, two such 
points cannot lie on the same level of first integrals if the area constant is nonzero. The possibility 
of coincidence of values of all first integrals at distinct rank -0 points is important in the study of 
bifurcation diagrams of various moment mappings that appear in this problem. Note that in the 
literature there is no explicit and clear proofs of these two assertions. 

Proposition 8. For all values of first integrals, the connected component of the integral manifold 
cannot contain more than one rank-0 critical point. 

Proposition 9. One of common levels of first integrals (different components of this level) contains 
two rank -0 points corresponding to the curve Iq on the plane (r, A). All other common levels of first 
integrals contain a rank-0 point, and such a point is unique. 
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Proof. All rank-0 point belong to the subsystem .Mi. We fix A / 0, By (1.25), a rank-0 point is 
uniquely determined by the values h, p, and r, where r is a multiple root of the polynomial (1.26). If 
we assume that there exist two multiple roots, i.e., that 4 R(r) = —(r — ?q) 2 (r — r 2 ) 2 , where r\ / r 2 , 
then we obtain the system 

V 2 = — ri, p = 0, r 2 = 2 h, r| — 4 h 2 + 4 = 0. 

Thus, for given A, h, and p a rank-0 critical point is unique (if exists). Now we fix A and h and assume 
that a (£, k ) is determined by two different p\ p 2 - From (1.30) we obtain two possibilities. The first 
possibility is p\ = — p 2 , and from the relation h — A 2 /2 = p\ we obtain £ = 0; this is the curve £q whose 
preimage contains two points with the same r but with opposite p. If P 1 +P 2 / 0, then we obtain the 
system 

2 2 , A 2 2 2 K A 2 . 

Pi +P1P2 +P2 = h - P1+P2 = y), 

which is obviously inconsistent. Thus, if R(r) has a multiple root, then this root is determined by 
h and p uniquely and, moreover, these values are uniquely determined by £ and k. Proposition 9 is 
proved. □ 


Note that for £ 0 this also implies Proposition 8. For £ = 0, the proof of Proposition 8 easily 

follows from Gashenenko’s solution presented above. 

Thus, in the Kovalevskaya-Yehia case, all rank-0 points have complexity 1. 

In [18], a weaker assertion was stated (without proof) in which cases that separate various types 
of bifurcation diagrams of the mapping (2.10) in the plane (£, A) are not considered. This separating 
set ®l was found by Ryabov [27-29]. In [18], nonseparating values of the pair (£, A) were called 
nonbifurcation. As the proof of Proposition 9 shows, the bifurcation character of the pair (£, A) is 
irrelevant to the complexity of a rank-0 point. 

The equations of the curves of the Ryabov separating set ®l can be found, e.g., in [46]. This set 
classifies bifurcation diagrams of the mappings (A). We calculate the set 0, i.e., the image in the 
octant {(£, A) : £ > 0, A > 0} of the curves 7 r^, which are separating curves in the classification of 
rank-0 points, together with their limit points for A = 0. Keeping the notation of the images of curves 
we obtain 

71-21 : 0<A<1; 


+23 


7T22 : £= -j= (yi + A 4 -A 2 J 


2 A 3 / 2 


A > 0; 


£ = 
A = 


(4 - x 4 ) 3 / 2 
4x 3 
3x 4 — 4 
2x 3 ’ 


x € 


y^,V2 


7T 2 4 : £ = 


| V4 + A 4 / 3 — 2A 2//3 | 

\/2 ^\/4 + A 4 / 3 — A 2 / 3 ^ 1/2 


A > 0; 


(3.45) 


+31 


(4 — .x 4 ) 3 / 2 


A = 


4x 3 

3x 4 — 4 
2x 3 ’ 


x € 


-^ 4 / 3,0 


Recall that the curve +) = {r = 0, A > 0} is not a separating curve in the class 82 in the sense of 
the equivalence relation introduced above. Therefore, its whole image {£ = 0, A > 0} is not contained 
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in the separating set: @l contains only endpoints of the curves 7T2 j : A = 0, 1,V2. Thus, we see that 
&L \ © consists of one curve £ = (4A) _1 whose sense will be clarified below. It was proved in [50] 
that this curve corresponds to the extremal value of the integral L in the family of rank-1 degenerate 
critical points. Therefore, the majority of surgeries of bifurcation diagrams of the mappings Jn{ A) 
occur in the cases where the reduced system possesses a degenerate rank-0 critical point, whereas the 
set &l of bifurcation pairs (£, A) is, except for a single curve, the image of degenerate rank-0 points. 
The set ©l is irrelevant to the possibility of the situation where several rank-0 critical points are 
contained in the same integral level. 

Note that Theorem 9 in [18] (and the same Theorem 8 in [32]), in some sense, plays a key role 
for the nondegeneracy of rank-0 points corresponding to nonbifurcation values (£, A). This fact is the 
basis of their further classification and description of marked loop molecules and, as a consequence, 
we can find some marks on Fomenko graphs. However, this part of the assertion, which is obvious in 
the analysis of Eqs. (3.45), has not been proved in [18, 32], It was indicated that for all nonbifurcation 
values (£, A), the verification of nondegeneracy conditions can be performed by a computer. However, 
the characteristic polynomials of the corresponding symplectic operators cannot be expressed through 
£ and A. The same holds true for the verification of the condition that the subalgebra generated by 
the operators an and <Xk is two-dimensional. It seems that this subalgebra is always two-dimensional, 
but this hypothesis is invalid. Above we have presented a set of values of parameter for which the 
subalgebra is one-dimensional. 

We also note that the paper [32] contains a reference to [3], concerning, in particular, information 
on the number of families of tori in chambers, which allows one to construct loop molecules for rank-0 
points (in [18], the number of families is invalid). The paper [3] also does not contain a proof since it 
is based on numerical simulations with a hypothetical algorithm that allows one to choose one initial 
point on each of the corresponding regular tori (and, in particular, to find the number of these points) 
for an arbitrary set of constant first integrals. However, such an algorithm has not been constructed 
in [3]. The following Ryabov hypothesis was also used (this is known from reports in conferences; see, 
e.g., [2]): each regular torus in P 5 is related to the subset { 0 J 2 = 0, 02 = 0}; this hypothesis has not 
been proved. 

The papers [30, 44] appeared before [32], 

In [30, 44], explicit exact analytic expressions for eigenvalues of symplectic operators were obtained; 
these expressions allow one to find the type of critical points. The paper [44] contains an exact analytic 
classification of rank-0 points with respect to the degenerateness, with respect to their types in the 
nondegenerate case, and with respect to the complete common level of first integrals, and the first 
(and unique at the present moment) rigorous proof of the nondegenerateness of rank-0 points outside 
Ryabov’s separating set. 

3.3. Smale diagrams and isoenergetic surfaces. For the Kovalevskaya problem (A = 0), the 
bifurcation diagram for energy and area integrals was constructed by Jacob [60]. Using the Smale 
construction (the reduced bundle of unit spheres over the domain where the motion is possible), 
he also found the topological type of isoenergetic manifolds that are three-dimensional levels of the 
“reduced Hamiltonian” (3.2) 

Qlh = {C e : He(C) = h}. (3.46) 

In this case, the diagram Slh of the mapping 

LxH : P 5 -> M 2 


consists of the two parabolas 


<5? : h = -l + £ 2 , d° 2 : h = l + f 
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Fig. 4. Smale diagram of the classical problem 


and the pair of curves whose equations can be written in the parametric form: 

x 2 + 4 


5%: £ = 


,3x1 

h = — -\ , x € (0, 2 ; 

4 x 


(3.47) 


4y/2x ’ 

these curves are symmetric with respect to the axis Oh. These values are attained at relative equilibria 
whose phase coordinates are 


lx 

W1 = "V2’ 

u 2 = 0 , 

V4 — x 2 

^3 = 0 r 

2y/X 

X 

o' 

II 

CN 

V4 — x 2 

ai = ~2’ 

«3- 2 


where yfx and \/4 — x 2 are algebraic radicals. The curves (3.47) touch the upper parabola at the 
points (±1,2) and pass it transversally at the points 

^± \ J 2(\/2 — 1 ), 2\/2 — 1 ^ . 

The cusps have the coordinates (±2/3 3 / 4 , \/3); they correspond to x 2 = 4/3 (note that = 4/3). 

The isoenergetic surfaces are empty in the domain h < — 1 ± l 2 . For other domains they are 
diffeomorphic to the following manifolds (see Fig. 4): 

S 3 , K 3 = (S 2 xS 1 )#{S 2 xS 1 ), S 2 xS\ RP 3 . (3.48) 

The smoothness type Q 3 h at each point (£, h) can be found by the index of the Morse function 
at its critical points lying in the preimages of bifurcation curves that arrive at the point (£, h ) along 
a vertical line from a point that lies sufficiently low and has a nonadmissible value of h (i.e., from a 
point Q\ h =0). 

Now we recall the results for the classical problems (see, e.g., [8, 38, 39, 57]; note that explicit cal¬ 
culations have never been published). The characteristic polynomials of the operator an are obrained 
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from (3.33) by passage to the limit: 


^1 : Xh{v) — ^ [/^ 2 + (1 + ^ 2 )] i 

$ ■■ Xh(v) = (m 2 ~ ^ [M 2 ~ (1 - ^)] , 

<^3 : XH(fj) = {h 2 + r 2 ) p 2 - 2 (V4 + r 4 - 2r 2 ) . 


(3.49) 


Therefore, on the lower parabola, all rank-0 critical points in the preimage have type “center-center,” 
whereas on the upper parabola they have type “saddle-saddle” on a bounded segment between two 
points of tangency with the third curve T ; these points are symmetric with respect to the axis Oh (in 
particular, rotations where the center of mass occupies the highest position for \(\ < 1 are unstable 
with respect to all variables). Further, these critical points have type “saddle-center” on unbounded 
segments outside the points of tangency (rotations where the center of mass occupies the highest 
position for \(\ > 1 are stable with respect to some variables). The preimage of each point of the 
parabola contains one critical point of such type. On the curves (3.47), the preimage contains two 
points of type “saddle-center” for r A <4/3 (on bounded segments between the cusps C and the points 
of tangency T) and two points of type “center-center” for r 4 > 4/3 (on unbounded segments from 
the cusps to infinity). The type is independent of the passing through a transversal point I along a 
smooth branch. 

Due to the mechanical nature of the Hamiltonian H(, its Morse index is equal to the Morse in¬ 
dex of the “effective potential,” i.e., the function defined on the Poisson sphere whose sublevels are 
domains where the motion is possible (domains of possible motion). The effective potential for the 
Kovalevskaya—Yehia case calculated by the Srnale scheme has the form 

rr _ , (2£ — Aa 3 ) 2 

'off)/ 2 i 2\ i 21* 

2[2(af + a 2 ) + a 3 ] 

For the coordinate-free calculation of the Morse index of the restriction of the function of three 
variables /(au, 0 : 2 , 013 ) to the Poisson sphere (1.3), we use the following assertion. 


Lemma 3. Consider the differential operator 


a = a x 


(3.50) 


that generates the second group of equations (3.3). Let ckq € S 2 = {a : |a| = 1} be a critical point of 
the restriction of the function /(a) to S 2 , which is nondegenerate in the sense of Morse. The Morse 
index of the function f at the point ccq is equal to the number of negative roots of the polynomial 


£/(m) = - det (S 2 /)(«o) - h E 

u, 


Applying this lemma to the function Ue o, we obtain 

<5? : (^) = (M - 1) [h - (f 2 + 1)], 

<5° : (m) = (M + 1) [h - ~ 1)], 

^3 • £H(ji)= L + Tr^ + r 2 ) n - ^= (V4 + r 4 - 2r 2 ) . 

We obtain the arrangement of Morse indexes and types along bifurcation curves shown in Fig. 4. The 
notation n * m means that the preimage contains n points of index m, and the symbols c + c, s + c, 
and s + s show the type of the point (“center-center,” “saddle-center,” or “saddle-saddle”). 

It is known that if the value h passes through a critical value, then the following surgeries of 
the domain of possible motion (projections of the energy level to the configuration space) occur: 
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supplement of the disk D 2 for index 0, pasting a handle for index 1 (a single disk becomes a disk with 
a hole, i.e., a ring, or produces a single disk from two disks), and sealing a hole with a disk for index 2 . 

We draw a vertical line t = const on the plane R 2 (£, h) between the points I and T. Along this 
line, the Hamiltonian and the effective effective potential have the critical values 

hi = — 1 + P < /i 2 < h -3 = 1 + f 2 < /14; 

the numbers of critical points in the preimage are 1 , 2 , 1 , and 2 , respectively, and their indexes are 0 , 
1, 2, and 2. The corresponding domains of possible motion on the sphere are as as follows: a disk, 
a disk with two holes (a sphere with three holes), a ring (a sphere with two holds), and a sphere. 
Reduced bundles of unit circles over these domains yield the manifolds (3.48), respectively. 

In the general case, calculations for A > 0 yield the following Morse indexes. 

Theorem 6 . The roots of the characteristic polynomial f,H(h) it the critical points (3.8) (in the 
notation of Proposition 2) are as follows: 

hi = -\ r (r — X) + d , 

h 2 = i( 2r - A)(r — A) — d] \( 2 r - A)(r - A)r + A d . 

2(r — A )d L J L J 

In particular, the sign of n 1 is always opposite to the sign ofd ; therefore, Hi is positive in the domain 8 \ 
and is negative in the domains 62 and 63 . The sign of p 2 is determined by the position of the point 
(?’, A) relative to the separating curvesT, ttq, 7 ^ 3 , 7724 , and 7731 . Finally, the Morse index of the effective 
potential U(__\ and the Hamiltonian H^\ is equal to 

(i) 0 in the domain 5 \; 

(ii) 1 in the domains 623 , $ 27 , and £ 31 ; 

(iii) 2 in the domains 821 , 822 , 824 , <$ 25 , 4' 26 , 823 , and 832 . 

In Fig. 5, we present the pairs of Morse indexes and information necessary for for analysis of 
modifications of the Srnale diagram for nonzero A (see below). 

Remark 6. In Fig. 5, we use the following notation for two key values of the parameter A: 

A* = l/ 2 3/4 , A* = (4/3) 3/4 . (3.51) 

The first of them is the ordinate of the node point Pq, where several separating curves meet. The 
second value is the minimum of A on the curve Iq determined by Eq. (3.43). These important values 
will be useful in the sequel. 

The evolution of Srnale diagrams is as follows. When A becomes nonzero, the curve 8 ® is transformed 
to <5i; it ceases to be a parabola, but no singular points appear on it. Each connected part of the 
curve followed by two critical points splits, and one of the segments obtained from the bounded 
segment joins with the branch of the curve <5!) lying above the point T. The curve £3 appears: each 
of its components is symmetric with respect to Oh and consists of twi infinite branches that meet at 
the cusp (the perturbation of C ). The second part of $3 and the curve 8 ® generate the curve 82 ; the 
triangle CTI on it generates a swallowtail on it. Taking into account the Morse indexes known, we 
obtain the situation shown in Fig. 6 . This type of diagram is valid for all A € (0, Ai), where Ai is the 
value of the parameter such that the rightmost cusp on the curve 82 (corresponding to the separating 
curve 7 t 2 4 ) lies on the curve 83 . The existence of such a value was predicted earlier on the basis of 
numerical simulation [27-29], but nobody can calculate it analytically. Modern system of analytical 
calculations allows one to perform this. We write the system of equations starting from (3.11) and 
(3.40) (we must choose the sign d > 0): 

£(r,X) = £(r 0 , A), h(r, A) = h(r 0 , A), r 0 = ^ (a - A 1/3 V 4 + A 4 / 3 ) . (3.52) 
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Fig. 5. Morse indexes of the effective potential 



Fig. 6. Diagram Slh for small A 
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Note that it describes not only all possible cases where the point 7^4 lies on the curve 63 , but also 
the hypothetic possibility of the case where this point lies on the other branch of the curve S 2 . We 
eliminate vq and r from this system assuming that ro / r. We proceed as follows. First, we eliminate 
the radical of d in the first two equation, substitute ro, and perform the change of variables U = A 1 / 3 . 
The two equations obtained can be cancelled by (2 r — t/ 3 + UV 4 + U 4 ) 2 (this corresponds to the case 
r = ro). We arrive at the system 

Pi3(U, r) + (U 3 - r)ps(U, r)^4 + t/ 4 = 0, 

_ (3.53) 

P25 (U, r ) + ( U 3 - r)p2o(U, r)V± + U 4 = 0 , 

where Pi(U, r) is a polynomial of degree i with respect to U. We calculate the resultant of the left-hand 
sides with respect to r in the substitution X = ( U 2 + \J 4 + U 4 ) 2 ; this corresponds to the change of 
variables 

X = (A 2/3 + y/4 + A 4 / 3 ) 2 . (3.54) 

We obtain the equation 

(X - 8 ) 2 (X 2 - 16) 4 (X + 4) 2 P 1 {X)P 2 (X) = 0, 

where 

Pi = X 4 - 24X 3 + 720X 2 - 2048X - 3072, 

P 2 = 4X 15 - 293X 14 + 7864X 13 - 70320X 12 - 831232X 11 + 26316032X 10 - 263235584X 9 
+ 1223192576X 8 - 2241200128X 7 + 323747840X 6 + 1465909248X 5 - 16521363456X 4 
- 25736249344X 3 - 74155294720X 2 - 75161927680X - 17179869184. 


Under the obvious condition X > 4, which is implied by (3.54), the equation P 2 (X) = 0 has a 
unique real root X ~ 11.2707, for which A ~ 1.1268 and the unique common real root of Eqs. (3.53) 
is r ps 0.002089 € ( 0 , A), i.e., we do not obtain points from £ 2 , 3 . If we assume that we have obtained 
a point on <5i, then we must choose the sign d < 0 in Eqs. (3.52), but then the values (r, A) obtained 
above do not satisfy these equations. Thus, all roots of the polynomial P 2 (X) are extraneous. The 
condition X = 8 leads to the separating value A* = 1 / 2 3 / 4 , which corresponds to the meeting of both 
cusps and the self-intersection point of the curve d 2 . The equation P\ (X) = 0 has a unique admissible 
root X ~ 4.3418 for which A ~ 0.02349 and a unique common real root r ~ 1.47328 € (A,+ 00 ) of 
Eqs. (3.53). Therefore, this value of A corresponds to the required case 7724 G 63 . In addition, we have 
proved that the cusp 7^4 cannot he on the other branch £ 3 . We denote the separating value A by 
Ai ~ 0.0235. The exact result is as follows: 


Ai — 



1 1 

12 P2 + 2 



3823 _i/3 

~r Pi 



+ 


16393 A 
9^2 / 



pi = 1045767 + 183872\/34, p 2 


_410 - 34407 p 1 1/3 + 9 p\ /3 . 


(3.55) 


For A > Ai, the curves 5 2 and ^3 do not have other common points (see Fig. 7). In particular, 
the domain in which Q^ h = K 3 (the superposition of the “angle” 63 and the “tail” of the curve 5 2 ) 
disappears. 

Further surgeries of the Srnale diagram are related only to evolution of the curve 5- 2 . In this case, 
degenerate rank-0 critical points generating the separating curves 7121 and tt 22 do not affect the diagram 
itself since they lie on its smooth segments. On passing through A*, the “ends of the tail” are swapped 
(see Fig. 8); this can be easily seen if one draw a horizontal line (A = const) in Fig. 5 and move it 
upward. 
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Fig. 7. A fragment of the diagram Slh for A > Ai 
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Fig. 8. Passing through A* 


Now we examine the possibility of the case where the cusp of the curve 7^3 lies on the other branch 
of the curve 82 - For this case, based on (3.37), we write the system 

£{r, A) = £(r 0 , A), h(r, A) = h(r Q , A), 

r ° = X 2 x ^ ' A = 3X 2x 3 4 ’ x4 _ 2r>3;3 - 4 7^ 37 € ^-^4/3, \/2 . 

Eliminating r'o, r, and A from this system, we arrive at the equation for X = x 4 : 

(X-2) 2 Q 1 (X)Q 2 (X) = 0, 

where 

Qi = 3X 4 + 32X 3 - 180X 2 + 96X - 64, 

Q 2 = 6903X 11 - 153216X 10 + 1489200X 9 - 9324352X 8 + 44169408X 7 

- 160186880A 6 + 425104384X 5 - 806682624A 4 + 1108361216A 3 

- 1120534528A 2 + 784072704X - 285212672. 
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The equation Q 2 (X) = 0 does not have roots in the interval X € (4/3,4] and the equation Qi(X) = 0 
has exactly one root in this interval: 


x = i 

3 


q\ = 2951 - 408v 34 


308 + —+ 3 9l + —-g 2 -8 

1l 12 


V3 ( 435 

, 12 = ( 154-3<?1 

V 1i 


1/2 


Respectively, in terms of A, we have the equation 

64A 16 + 2784A 12 - 274803A 8 + 15476896A 4 - 45349632 = 0 


with a unique real solution A 2 ~ 1.3263. Its exact value is 


A2 


Pi 

P2 

P3 


( 87 1 1 A 1/4 

[-T- 16 P 2 + 2 P3 ) ’ 

10467417865895 + 1809781698048\/34, 

\J 213478 - 1093516839p]" 1/3 + 9 p{ /3 , 

1 106739 1093516839 _i/ 3 9““^ 88449457 

V 16 + 64 Pl ~ 64 Pl + 16^2 


(3.56) 


Other changes are related to the axis Oh (see [29]). The constant £ vanishes not only at equilibrium 
positions of the body (r = 0), but also under the condition (3.43) that follows from (3.11). The 
corresponding curve (above it was denoted by £q) is shown in Fig. 2. The minimum of A on this curve 
is equal to A* = (4/3) 3//4 . Therefore, except for the points h = ±1 and l = 0, which always exist, the 
diagram has two additional points of intersection with the axis Oh for A* < A < \[2 and one point for 
A > \/2- This collection of points on the axis Oh is reconstructed at the moment when two of these 
points coincide. One can easily calculate that this occurs when A = 2y \/2 — 1 ~ 1.2872 € (A*,\/2). 
Thus, we have the types of diagrams whose significant fragments are shown in Fig. 9: 

(a) A* < A < A*;_ 

(b) A* < A < 2 yV2- 1; 

(c) 2 VV 2 - 1 < A < A 2 ; 

(d) A 2 < A < \/2; 

(e) A > y/2. 

Taking into account the Morse indexes calculated above, we obtain the arrangement of isoenergetic 
manifold shown in the same figure. 

We prolong the symmetry (3.17) to the whole phase space: 


symm : (uji,uj 2 ,a 3 ) (-uq, -w 2 , -« 3 )- (3.57) 

It states an isomorphism of phase flows on Qj h and Recall that in correspondence with (3.19), 

the extended Srnale diagram 

M$lh) = U‘5lh(A)x{A} 
a 

appears in the space A(M 2 (£, h)) = M 3 (£, h, A). It divides the space M 3 (£, h, A) into open connected 
components called chambers. Due to the symmetry symm, we also accept the agreement that the 
union of two components symmetric with respect to the plane l = 0 forms a single chamber. We 
obtain the following assertion. 
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Fig. 9. Surgeries of Slh on the axis Oh 
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Fig. 10. Chambers of the Srnale diagram 


Theorem 7. In the Kovalevskaya-Yehia case, there exist seven Smale diagrams Slh(X) that are 
structurally stable with respect to the parameter A . The separating values of the parameter A are 

0, Ai, A* = 1/2 3 / 4 , A* = (4/3) 3 / 4 , 2 \J V 2 — 1, A 2 , a/2, 

where \\ and A 2 are defined by the relations (3.55) and (3.56). The extended diagram A (Slh) divides 
the space M 3 (£, h, A) into eight chambers A, ...,H with nonempty manifolds Q^ h (X). Corresponding 
information is represented in Fig. 10 and in Table 3.1. 

Now we examine generic critical points that appear in rank-1 systems with two degrees of freedom. 
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Table 3.1 



4. Classification of Rank-1 Critical Points 

4.1. Formulas for calculation the type. We study the type of rank-1 critical points in the first 
critical subsystem, i.e., points of the set’Ali\C 0 . 

Proposition 10 (P. E. Ryabov [30]). The type of rank-1 critical points in the first critical subsystem 
is completely defined by the eigenvalues of the symplectic operator of i; where 

Fi = K - 2 vH, (4.1) 

and the undefined multiplier at the corresponding point is equal to v = p 2 . The characteristic polyno¬ 
mial of the operator a p has the form 

X,F\ ip) = p 2 - 4Ci, Ci = y - [h - y'J 2 s 2 -2^h+ s + 1 . 

Rank-1 points are degenerate if and only if 

2 s ~ (j 1 — 9 ”) 2s 2 — 2 ^h +—''j s + 1 =0 (4.2) 

and have type “center” for C i < 0 and type “saddle” for C\ > 0. 

Proof. In the coordinates (u>, ck), the Hamiltonian field generated by the function F\ at the points 
(1.25) has the form 

sgrad F\ = 2 (p 2 - v) |o, 0, yjR(r ), ry/R(r), ~R'(r), -p j ; 

it vanishes for u = p 2 . We note that for other values u, the condition sgrad F\ = 0 implies that R(r) 
and R'(r) simultaneously vanish, which leads to immovable points, i.e., rank-0 critical points. 

Recall that due to (2.13) 



Owing to Eqs. (1.25), the calculation of the polynomial XFi(p) for the manifold M\ is easy. Thus, if 
Ci 0, then rank-1 critical points are nondegenerate. Let Ci = 0. Assume that at such a point, there 
exists another combination G = uiK — 2 v- 2 H satisfying the relation sgrad G = 0. By the assumption 
rankjRT, Lf}|p4 = 1 we obtain that the undefined multipliers are proportional: 

dF = dK — 2udH\ P 4. = 0, dG = u\dK — 2u-2dH\ P A =0 V 2 = v\v. 
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This implies that v\ fi 0, i.e., the characteristic numbers of the operator ac are proportional to 
(zero) characteristic numbers of the operator ctpj and hence they are also equal to zero. Therefore, 
the integral generating a regular element in the algebra of symplectic operators does not exist at the 
points (4.2). □ 


We have presented this proof since the results of [41] (where the function (2.16) was used for 
the construction of the system A4 1 ) cannot be applied immediately. The roots of the characteristic 
polynomial also vanish on the set s = 0, which implies the relation £ = 0. The appearance of this 
extraneous singularity is explained by the degeneration of the Reinran-Semenov-Tyan-Shansky integral 
(see [41]) in L' 2 . 


Proposition 11 (P. E. Ryabov [30]). The type of rank-1 critical points in the second and third critical 
subsystems is completely defined by the eigenvalues of the symplectic operator ap 2 , where 

F 2 = K + ^2A 2 - ^ H (4.3) 

and s is the value that determine the point (1.28). The characteristic polynomial of the operator of 2 
has the form 


c 2 = --L (8AV - 1 ) 


XF 2 {P)=P -^ 2 , 
2s 2 — 2 (h -\— s + 1 


j(8AV- l)(2 AV-s + 2t' 2 ). 


L \ * / 1 0 

All rank-1 points in A4 2 are nondegenerate and have type “center.” Rank-1 points in M .3 are degen¬ 
erate if and only if 


(8A 2 s 3 - 1) 


A 2 

2s" — 2 ( h H—— ) s + 1 


L \ * / 

and have type “center” for C 2 < 0 and type “saddle” for C 2 > 0. 


= 0 , 


Proof. Setting F 2 = K + nH , we obtain from (1.27) and (1.28) 

{QxpX pn\ 2 QkY Q{n\ + 2x 3 XY) 


sgradF 2 = (2A 2 s — 1 — ns) 


\ \fs ’ 2 XS* y/~S 


Xy/s 


U 2 g(A psX + £Y) 

S 2 ’ yfsx 


where 

m = -ipX + A sY + 2x i XY, u 2 = kX + sXY(3£p - 4x 3 Y) + As 2 (2 - 3y 2 ). 
If 2 A 2 s - 1 — ns 0, then this vector vanishes only under the condition 


64£ 6 X 2 s 4 + 4 i 4 s 2 ( - 1 + 20A 2 s - 4A 4 s 2 + 48A 4 s 4 ) 

- 4^ 2 s(l - 8 A 2 s + 4A 4 s 2 + 20A 2 s 3 - 4A 4 .s 4 + 8A 6 s 5 - 48A 6 s 7 ) 

- (1 - 4s 2 ) (l - 2 A 2 s + 4A 4 s 4 ) 2 = 0. 

This condition can be obtained as follows. Note that in the case where X = 0 and Y = ±1 the vector 
sgrad F 2 does not vanish even if p = 0. Therefore, Q = 0. Expressing X from this relation, we 
substitute in relations for n\ and n 2 , calculate the resultants with respect to Y. and find the required 
relation for £ and s. Since the condition found does not become an identity on the surfaces 112,3, it 
determines a one-dimensional subset. On the other hand, it becomes an identity for the corresponding 
expressions of £ and s through r and A from (3.11) and (3.15); therefore, it corresponds to the curves 
Si, i = 1,2,3, i.e., it is valid at rank-0 critical points. Thus, at rank-1 critical points, the relation 

2A 2 s — 1 — ns = 0 

holds, and the function (4.3) is a required integral that determines the type of these points. 
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The direct calculation of the characteristic polynomial leads to the required expression for x_f 2 (m). 
Obviously, for s < 0 the value C 2 is always negative, whereas for s > 0 it can have various signs, which 
determine the corresponding type of singular points. The fact that critical points are degenerate 
for C 2 = 0 (i.e., there is no another element in the algebra of symplectic operators that has nonzero 
eigenvalues) is proved similarly to Proposition 10 (actually, this follows from the fact that the coefficient 
of K can be chosen equal to 1). □ 

We see that for all critical points (rank-0 and rank-1), the characteristic polynomial of the defining 
symplectic operator depends only on constant first integrals. Therefore, the following important 
property is valid. 

Theorem 8. All critical points of the same rank lying on the same common level of the first inte¬ 
grals L, H, andK have the same type. 

In particular, this implies that the molecule denoted by Wy in [14, 16] is impossible in the problem 
considered. Below we show that this molecule cannot be realized by a significant cause. 

We introduce some terms related to the study of critical subsystems. 

It is known (see [8]) that geometric singularities of bifurcation diagrams of integrable systems with n 
degrees of freedom correspond to rank-m (m < n—1) critical points and rank (n— 1) degenerate critical 
points. It has not yet been proved in the general case that this correspondence is bijective. Additional 
singularities appear in families of Hamiltonian systems obtained by the factorization procedure of a 
system with symmetry. The problem considered here is an example of such situation. The surface Hi 
has a singularity of the type “self-intersection,” which appears in all planar sections of it except for 
the sections £ = const (the curve s = 0 in the representation (2.7) and the singular parabola in the 
representation (2.10)). Among sections by a fixed constant of area, the section £ = 0 is singular due 
to the presence of a singular “multiple” parabola. 

Definition 8. A nondegenerate rank-1 critical point is said to be multiple if the bifurcation diagram 
£ of the moment mapping J = LxHxK is not a smooth two-dimensional surface in a neighborhood 
of its image. 

Recall that the rank of a critical point is its rank in the reduced system with two degrees of freedom; 
this number is one less than its rank with respect to the mapping J. 

Since the structure of a neighborhood of a nondegenerate rank-1 critical point in Qf h is similar to 
the structure of one of standard atoms and can preserve its structure under small variations of h, we 
conclude that such a multiple point cannot be detected by a local analysis of singularity in the reduced 
system with two degrees of freedom on Pf. We can only conclude that to some points on arcs of the 
bifurcation diagram of the system on Pjr, several critical circles correspond. In this case, the decay of 
a multiple point can occur under a perturbation of £, i.e., under the passing to another system with 
two degrees of freedom. In the general integrable cases of problems of the dynamics of a rigid body, 
the zero level of the cyclic integral L is always singular in this sense. There are other causes for a point 
to become a multiple point. For example, as in the problem considered, rank-0 points can lie on the 
same level with rank-1 points. Then this “multiplicity” is irrelevant since rank-0 points also affect the 
surgery diagram. Finally, the system can also contain “decaying” topologically instable atoms. This 
effect was described in [8]; it also cannot be detected by local analysis. Note that topologically instable 
atoms appeared in specific mechanical systems of the classical dynamics of rigid bodies (see [40]), but 
the decay itself did not occur. The decay of a topologically instable atom was analyzed by numerical 
simulation in [22] for a system on the sphere (the so-called Dullin-Matveev case, [58]), which is a 
particular case of an integrable system constructed by Yehia [70]). Various versions of similar surgeries 
of fibrations on three-dimensional isoenergetic levels were described in [66]. Recently, decaying atoms 
were detected in subsystems of an irreducible problem on the motion of a rigid body in a double field; 
however, this effect is closely related to degeneration of the induced symplectic structure (see [63]). 
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Definition 9. A key set in P 5 is the union of all rank-0 critical points and degenerate and multiple 
rank-1 critical points. We denote the key set by /C. The intersection JCj = K, H Mj is called the key 
set of the critical subsystem Mj. 

Definition 10. Let and T be first integrals (possibly, particular) of the critical subsystem Mj that 
are independent almost everywhere. The image of the key set ICj under the mapping <hxT : Mj —>• M 2 
is called the (<L, 'L)-diagram of the subsystem Mj 

Definition 11. If an integral mapping Mj —>• M fc exists, then its image is called an admissible domain ; 
notation Vj. Points of the admissible set is said to be admissible. 

In other words, values of integrals are admissible if they correspond to some motions (the integral 
manifold is nonempty). It will be clear from the context which mapping is considered and to which 
space admissible points belong. 

In [48, 50], ( S,H )- and (S, L)-diagrams of critical subsystems were studied. Here S' is a particular 
integral that appears from the Lax representation, whose constant is a parameter on surfaces in 
Eqs. (2.7) and (2.8). If we introduce the function with undefined Lagrange multipliers to describe 
critical subsystems by the formulas (2.16), (2.17), and (2.18), then this integral naturally appears as the 
coefficient of the integral K. For the classification of bifurcations that appear when the point (£, h, k ) 
crosses the surfaces n y in M 3 (f, h, k), it is convenient to choose a plane of constants of functionally 
independent integrals of critical subsystems whose admissible points are in a bijective correspondence 
with points of the corresponding bifurcation surfaces (then, for example, information on the number 
of rank-0 points or rank-1 circles in the preimage becomes definite). 

We also consider the problem of classification of bifurcation diagrams Y,hi<{£, A) of the mappings 
H xK of reduced systems on Pf (see [46, 50]) and bifurcation diagrams TiLi<{h, A) of the mappings 
LxK restricted to four-dimensional isoenergetic levels Qj t = Fd~ l {h) in P 5 (see [48]). These diagrams 
are the sections of the bifurcation diagram E(A) of the general integral mapping J = LxHxK by the 
planes i = const and h = const, respectively. Due to the free physical parameter A, one can consider 
the problem on the search for the so-called separating set in the corresponding plane (X,£) or (A, h). 
We state a general assertion that allows one to state an algorithm of construction of separating sets. 

Let F : P 5 —>• K be a general integral of the system. We fix / € M and consider the mapping 

J(F,f) = F-!(/) : F 1 (/) ^ 2 > 

where M 2 is the plane of values of a pair of integrals that supplement F to a complete involutive 
triple of integrals that are independent almost everywhere. Let &(F,f)( A) be the bifurcation diagram 
of the mapping J(F,f)- Clearly, if the pair of integrals U, V supplements F to an involutive triple of 
functionally independent integrals, we can naturally identify &(F,f){ A) with the bifurcation diagram 
Tjjjvifi A) of the restriction of the mapping UxV to the submanifold F~ 1 (f). In particular, below we 
consider the case where F = L (in this case, A) = T,hk(£, A)), and where F = Ff (in this case, 

a (H,h)( A) = T, LK (h, A)). 

Proposition 12 (M. P. Kharlamov [43, 45]). Consider the set Of in the plane (A,/) such that the 
type of the diagram a^ F f\(X) changes after passing through it. The set Of consists of pairs (A,/), 
where f is the critical value of the restriction of the function F to the key set 1C for a given A. 

If necessary, one can state this statement rigorously, clarifying the notions of diagrams of the same 
type and critical values of functions on a stratified manifold. In particular, we assume that the 
integral F is “rationally” selected and has no critical points on regular levels of the moment mapping. 
One of possible approaches using the concept of equipped admissible intervals was proposed in [50]. 
In practice, it means the following. For each critical subsystem Mj, we consider its (G, i ? )-diagram, 
where G is an integral on Mj, which is functionally independent of F. Critical values of F on 1C 
correspond to nodal points of this diagram (the images of degenerate rank-0 critical points), extreme 
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values of the coordinates / on smooth domains that are the images of the set of degenerate or multiple 
rank-1 critical points, and various self-intersections of smooth domains of the diagram. 

The set is called the separating set in the classification of bifurcation diagrams c r (F,/)(A). For 
brevity, we call it the F-atlas. 

4.2. Specification. The first critical subsystem. For the first critical subsystem, none xof the 
pairs of functionally independent integrals (S', H) or (S, L) solves the problem of providing a bijection 
of the selected plane and the surface ITi. Indeed, on the (s, /i)-plane, each admissible point satisfying 
the condition 



has two preimages on the surface IIi, whereas in the (s, ^(-representation, the whole admissible segment 
of the parabola 

/ \ 2\2 \2 

k = l+[h--j , e = 0, h> -, (4.4) 

corresponding to the value s = 0 on the surface IIi, is mapped to the same point. Thus, this 
correspondence is either two-valued on a set of full measure, or there is a “bad singularity,” i.e., a 
point possessing an infinite number of preimages lying on the surface. Note that since the surface IIi 
contains a self-intersection line, a bijective parametrization does not exist. However, one can provide 
a parametrization in which each point of the self-intersection line naturally corresponds to a pair of 
points on the parameter plane. For this, we choose a particular integral P = oj\ on Adi instead 
of S. Then Eqs. (1.30) yield a bijective correspondence of admissible domains on the (p, /i)-plane 
and the surface Hi, with the following obvious sole exception: any interior point of the segment of 
the parabola (4.4), i.e., for h > A 2 /2, corresponds to two points with opposite signs of p. Next, we 
consider the diagrams and key sets of the subsystem in terms of ( p,h ); taking into account the 
problem on the classification of bifurcation diagrams T*hk(£, A) and £la _ (^, A), we also present the 
description of the ( S,H )- and (S, L)-diagrams for M\. 

To find topological invariants, such as Fomenko graphs, one needs to know another indicator, namely, 
the Morse-Bott index of K integral at nondegenerate rank-1 critical points (i.e., the Morse index of 
the restriction of K to a small domain transversal to the trajectory consisting of rank-1 points). It 
gives information on the behavior of critical circles on the isoenergetic surface (A) if the value of 
the additional integral increases. In the space of these three integrals, we can trace events that occur 
when a line parallel to the axis Ok intersects the corresponding bifurcation surface Hi. The form of 
the functions (4.1) and (4.3)(1.3) allows one to perform the corresponding calculations explicitly. 

At points of the manifold, we find a couple of vectors that define a domain transversal to the 
critical circle (at the rank-1 critical point). We choose this domain so that it is orthogonal to the 
vectors gradT, gradL, grad#, and sgrad#. Note that on any trajectory (1.24), (1.25), the variable r 
oscillates between the roots of the polynomial (1.25); hence we can choose a point xq on the trajectory 
at which R(r) = 0. Then the tangent plane to the transversal domain is spanned by the vectors 

vi = ( 0 , 1 , 0 , 0 , 0 , 0 ), v 2 = ^A + r, 0 , -4 Ip, 2p(X - r), 0 , 2 (h - p 2 - 7 — 

These vectors are obtained from the equations for u = (u \,..., uq): 

u ■ grad T = 0, u ■ grad L = 0, u ■ grad H = 0, u ■ sgrad H = 0 

under the conditions ( U 2 , u 3 ) = (1,0) and (u\ . U 2 ) = (0, A + r). 

At the point xq, the function K achieves the conditional extremum on the common level Q\ h of 
the functions F, L, and H in P e . In particular, x$ is a critical point of the function with undefined 
Lagrange multipliers 

K\ = K — 2i'H — k\L — KyP . 
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Obviously, the part of this function that does not contain the Casimir functions L and T coincides 
with (4.1), i.e., v = p 2 . We directly calculate that = 4 p and X 2 = 1. 

The type of a conditional extremum is determined by the restriction of the second differential of the 
function with Lagrange multipliers to the tangent space of the restriction manifold (see [1]). This fact 
was used in the study of the phase topology of Hamiltonian systems by A. A. Oshemkov (see [23]). 

The restriction of d 2 K\ to the linear span of the vectors v\ and V 2 , under the condition R = 0, is a 
diagonal matrix with elements (eigenvalues) 


hi = 4 


h - (A 2 + p 2 ) + A r - -r 2 


= 2 


2s — (A — r)' 


A 2 

h2 = -8 ( h - — - ‘ip 2 
\ 2 

= -32 ( h - — - -s 
2 2 


A 2 


h-^-3p 2 -- 2 (X + rY 
(X + rf 


In particular the product 


H 1 /X 2 = 16( — 2h + A 2 + 6 p 2 ) 1 — (2 h — A 2 — 2 p 2 ) (A 2 + p 2 ) 
= -32(2 h - A 2 - 3s) [(2 h + A 2 - 2s)s - 1 


= -64 


-s — ( h — 


A 2 


2s 2 - 2 h + 


A 2 


s + 1 


is independent of r and is determined by the location of the point (s, h ) relative to the degeneracy 
curves of rank-1 critical points given by Eq. (4.2). Therefore, on nondegenerate trajectories, the values 
of hi and h 2 (they are called the Morse-Bott indexes) do not vanish and hence therefore do not change 
their signs. 


Proposition 13. If the integral K increases on the isoenergetic level Qj h , we have the following 

bifurcations at points of the critical subsystem M.\ on nondegenerate critical circles'. 

(1) for elliptic trajectory (type “center”): the appearance of a torus if n\ >0 and h 2 > 0 (the atom A 
with the edge up ); the disappearance of a torus if hi < 0 and h 2 < 0 (the atom A with the edge 
down)', 

(2) for a hyperbolic trajectory at a critical level of K for hi > 0 and h 2 < 0: the atom B with the 
“outer” edge ( “head”) up and a pair of “inner” edges ( “feet”) down', if H\ < 0 and h 2 > 0, then 
we have the atom B with an “outer” edge ( “head”) down and a pair of “inner” edges ( “feet”) up ; 

(3) for the two hyperbolic trajectories at a critical level of K if two pairs (hi, P 2 ) coincide: two atoms B 
whose directions of the “outer” edges is determined by the same rule (both “heads” up for fi\ > 0 
and h 2 < 0 or both “heads” down for ^2 > 0); 

(4) for the two hyperbolic trajectories at a critical level of K with different combinations of signs in 
the pairs (hi,h 2 ) ; two A* atoms. 


Proof. The proposition is obvious for elliptic trajectories (the function K has a local minimum or 
maximum on the transversal domain). 

It can be shown that for hyperbolic trajectories the vector v\ is directed into the exterior domain 
of the “figure-eight”: the direction of the axis Ou 2 is responsible for the transition from a critical 
surface to its enveloping torus. This can be easily obtained by analyzing, for example, projections of 
the integral manifolds on the plane Olo\lo 2 (this analysis was partially performed in [16]). Clearly, 
similarly to the classical problem (see, e.g., [4, 17]) projections never have discontinuities in the 
direction of the axis Ovj 2 - Now, if p\ > 0, then this means that the function K on the transversal 
domain increases in the direction toward the “outer” circle and decreases in the direction toward the 
pair of “inner” circles. 
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If the sets of signs in the pairs (ji \, /i-f) are different on two hyperbolic circles, then, suggesting the 
existence of two -B-atoms with opposite directions of “heads,” we obtain a bifurcation of three tori 
into three. As was noted in [14, 15], the number of tori at a regular level may be equal only to 1, 2, 
or 4 (for all chambers in M 3 (£, h, k ) \ X, the number of regular tori can be strictly calculated below by 
using only obvious A-type atoms). Therefore, in this case, we would have a bifurcation of four tori 
into four tori. However, as we will see below, there is no such chambers in this problem. If we assume 
that at the point considered, there is a C2-atom, then the analytic solution of problem (1.34)-(1.35) 
would contain a heteroclinic trajectory, but such a trajectory does not exist. Hence, there are two A* 
atoms at this point. □ 


The criterion of the existence of motions in the system Adi is obvious (see [50, Proposition 1]). 
Since there is no “double” points on the ( p , h)-plane, we can easily calculate the number of critical 
circles on the appropriate level of integrals. 

Proposition 14. For given p and li, real solutions of (1.25) exist if and only if R(r) > 0 for some 
r € R. If all roots of R(r) are simple ( i.e ., there are no rank-0 critical points on the integral level 
considered), the number of periodic solutions is equal to the number of intervals on which R(r) is 
nonnegative. 

Remark 7. Recall again that, due to the formulas (1.30), after the transition to the space M 3 (H, h, k ), 
any pair of points with opposite p from the preimage of the parabola (4.4) is mapped into the same 
point. In [14], this fact was noted as a special case. 

Remark 8 . Proposition 14 allows the following obvious generalization. If i?(?’) has a multiple root 
(as was shown above, in this case it is a unique root), the number of periodic solutions is the same as 
the number of intervals between simple roots in which R(r) is nonnegative. Now, taking into account 
Proposition 9, we obtain complete information on critical motions on the general level of common 
integrals that contain a rank-0 critical point of the corresponding class, since the number and location 
of roots of the polynomial R(r) can be easily examined for all classes, based on the dependences h(r, A) 
and p(r, A) defined by the formulas (3.11) and (3.12). We present this information in the third column 
of Table 6.1; it will be used in the further study of the topology of full levels of integrals for rank-0 
points. 

Introduce the following notation: 


<P±(r) 


. . . 2r — A 

r 1 

1 

r(A r)± D 

r — A 

-e- 

H- 

II 

to I 

-r ± - -D 

L r — A 


D = \d\ = y/4 + r 2 (r — A) 2 > 0. 


(4.5) 


We denote by hc( A) the value of the energy at the degenerate critical point defined by (3.38). 
The index C is introduced since in the sequel we will mark the image of this degenerate point in the 
diagrams of critical subsystems by C. For a given A > 0, the value of he is determined in accordance 
with (3.38) and (3.11) as follows: 


3x 4 


- 2.x 3 A -4 = 0, 
x€ (-</473,0) 


hc = - n x 2 + -,. 


x u 


(4.6) 
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Theorem 9. The (P, H)-diagram of the critical subsystem consists of the following sets: 

(5i : p 2 = 'if-{r), h = <p-(r), r<E[0,A); 

5 2 : p 2 = ip + {r), h = <p+(r), r<E(-oo,0]; 

5 3 : p 2 = ip + (r), h = <p+(r), r£{ A,+oo); 

Ai : h- ( 3 p 2 + =0, y < h < hc( A); (4.7) 

A„: h-(p + y)- 2(A^ +p 2) - 0 - 

4 : h - + y ^ =0, p € M. 

TTie boundary of the admissible V\ domain consists of the connected curve 4 and two connected 
components of the curve S3. The curve S 2 is connected and breaks the admissible domain into two 
subdomains. The points of the domain below the curve S 2 correspond to one critical circle, whereas 
the points of the domain above the curve S 2 correspond to two critical circles. 

Proof. Since all rank-1 critical points belong to A4i and the values of the functions P and H at these 
points are given by the formulas (3.8) and (3.11), we conclude, taking into account the notation (3.20), 
that the images of these points are the curves Si, i = 1, 2, 3. The curves Ai and Ao are obtained from 
Proposition 10 and (2.13). The curve 4 corresponds to the singular parabola (4.4), which, as was 
noted above, is the image of the set of multiple points. It remains to clarify the conditions of the 
existence of motions in the system A4 1 . These motions exist if and only if there exist intervals on 
which the polynomial R(r) defined by the formula (1.26) is nonnegative. Since its discriminant set 
consists exactly of the curves Si, we see that such intervals exist between the curves 4 and S3, whereas 
below S\ (in the direction of the axis OK) and above S3 such motions do not exist. Below the curve S 2 , 
the polynomial R has two real roots and hence there is only one interval (and one critical circle), 
whereas above this curve there are four roots and hence two intervals (and two critical circles). The 
part of the curve Ai on which rank-1 critical points are degenerate and which is contained in the 
admissible domain is bounded from below by the natural boundary h = A 2 /2 and from above by the 
point of intersection with the curve S 3 whose ordinate is the value hc(X) described above. □ 

Theorem 10. For nonnegative X, surgeries of the ( P , H)-diagrams of the critical subsystem M.\ occur 
for the following values of the parameters: 



Proof. Due to the choice of particular integrals that provide a bijective correspondence between the 
set key and its image (except for the cases noted above), if A changes, surgeries of the diagrams occur 
simultaneously with surgeries of the set of rank-0 degenerate points and multiple points. Since all 
curves of degenerate and multiple points are known (see Fig. 2), these surgeries are as follows: if 
A = 0, then the point 7T2 i disappears (the corresponding curve goes to infinity) and the points tt 22 
and 7T24 merge; if A = A*, then all curves of degenerate points of the class S 2 merge at the double 
point Pq: if A = 1, then the point 7T2i disappears (the corresponding curve terminates); if A = A* 
(the minimum of A on 4, see Remark 6), then multiple points appear; if A = y/2, then the point 7723 
disappears (the corresponding curve terminates). There are no other surgeries. □ 

In Fig. ll(a)-(e) the (P, P)-diagrams of the subsystem M 1 stable with respect to the parameter A 
are shown: 

(1) 0 < A < A*; 
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Fig. 11. (P, P)-Diagrams of the system Mi and its enlarged fragments 


(2) A* < A < 1; 

(3) 1 < A < A*; 

(4) A* < A < \/2; 

(5) A > y/2. 

The connected components of the admissible domain in the space {(p , h, A)} formed by the extended 
(P, P)-diagram are denoted by ai-ai 2 - We recall the existence of the symmetry (3.57). Domains that 
pass one into the other under this mapping are denoted by the same symbols. 

In terms of the types of critical points, the transition through the value A = A* does not imply 
any qualitative changes; the points of the domains aio and ai 2 and their common boundary £q are 
the same. For A = A*, the curves 62 and intersect; this demonstrates the exit point of the domain 
aio U ai 2 to the manifold L = 0 and allows one to compare bifurcations and atoms with the case L = 0, 
which was examined in detail by P. E. Ryabov and P. V. Morozov. Moreover, from the local point of 
view, critical points of the domains a\ and 02 and and their common boundary are also the same. 
However, points of the bifurcation diagram T, in M 3 (£, h, k) generated by these domains differs by a 
nonlocal feature found in [11], which, in terms of the parameters (1.33), consists of the presence or 
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absence of regular tori on the complete joint levels of common integrals (two-frequency motions). This 
will be discussed below. 

On flat sections h = const or £ = const (i.e., isoenergetic bifurcation diagrams or bifurcation 
diagrams of reduced systems with two degrees of freedom on Pf) of the surfaces Ilj, smooth arcs 
correspond to the domains a,j , where the type and number of singular points and the character of 
bifurcation are fully preserved, i.e., the atom that appears along a segment of small length that 
transversely crosses Ilj at the corresponding point. 

For additional information, see Fig. 11(f), which shows the diagram of the first critical subsystem in 
the classical problem (A = 0) corresponding to particularly remarkable motions of the fourth Appelrot 
class. It consists of two parabolas h = ±1 + p 2 , the curve obtained from £3 and the part of 82 


p = ± 


(^ITr 4 - r 2 ^ , h = ^ ^2\/ 4 + r 4 + r 2 ^ , 


and the limit curves 


1 


A 0 : h = - ( 2 p + — ) , p G M\{0}, Ai : h = 3p , p G 


0 , 




Analyzing all characteristics calculated for the domains aj, we collect information in Table 4.1. 
Here we introduce the following notation for the atoms A and B: the index “+” means that the 
number of tori increases when K increases, whereas the index ” respectively indicates decreasing. 
In other words, A_ means that the edge terminates at the top and A + means that the edge terminates 
at the bottom; F>_ is the “outer” circle of the “figure-eight” (“head”) at the top, B + is the “outer” 
circle of the “figure-eight” (“head”) at the bottom. We see that all domains except for 04 are related 
with the classical Kovalevskaya problem {£ = 0) and to the problem with gyrostatic moment at zero 
area constant (i = 0). For these cases, bifurcations are found in the papers [38, 46] (for a detailed 
presentation for the Kovalevskaya case, see [39, 40]), and marks on the corresponding molecules are 
calculated in [7, 20]. Thus, in the column “Analogs” we provide references to the notation of parts [7, 
20, 40] or paths crossing the corresponding parts [38, 46]. If an analog exists, all atoms were described 
in [38, 39, 46]. In these papers, a language for the description of atoms was proposed, including the 
orientation of asymmetric atoms (the concept of atoms did not exist at that time). Comparing them 
with analogs considered in [7, 20], one can specify many marked molecules. 

Now it is possible to clarify the relationship between the parameters (1.33), conditions (1.36), and 
the classification of points of the first critical subsystem given here. For brevity, we denote by £ 0 , Ai, 
and Ao the values of the left-hand sides of the equations (4.7) of the corresponding curves (this will 
not lead to confusion). Then, in the notation (1.33), we have 


L 



L2 


A! 

2 ’ 


L3 


Ap 
2 ' 


(4.8) 


From (4.5) and (4.7) we see that on the curve <5i, there is a unique relationship 


h = Hi(p), pGK, 


and on the curve S3, for nonzero p , there is a unique relationship 


h = H 3 (p), p G M\{0}. 

Moreover, for all p / 0, we have H 3 (p) > H\ (p). The curves 5i form the discriminant set of the 
polynomial (1.26). If 

h G (- 00 , F/i(p)) U ( H 3 (p ), + 00 ), (4.9) 

this polynomial does not have real roots and, therefore, motions on Adi are impossible. We show that 
under the condition (1.30), motions in the domains (4.9) are also impossible on the whole manifold G 4 
(we recall that G 4 is the complete preimage of the surface Hi in P 5 under the moment mapping J). 
For this, it suffices to note that at the point where the restrictions of the functions H and P = uj\ to 
the manifold G 4 are dependent, the rank of J does not exceed 2 , so this point is critical of rank 0 or 1 


284 



Table 4.1 


Domain (life time) Number 
of circles 


a i (0 < A < +oo) 1 


02 (0 < A < +oo) 


03 (0 < A < +oo) 


04 (0 < A < +oo) 


05 (0 < A < 1) 


o 6 (0 < A < y/2) 


a 7 (0 ^ A < A*) 


a 8 (0 < A < A*) 













Morse-Bott 

indexes 

Exit on 

a = o/e = 0 

Atom 

Analogs 

(- -) 

Yes/Yes 

A_ 

2,3 [40, Fig. 6.3] 
ai,a 2 [46, Fig. 2] 
71,74 [7, Fig. 11] 
a 2 ,<23 [20, Fig. 1] 

(- -) 

Yes/Yes 


3,3' [40, Fig. 6.3] 
a 2 [46, Fig. 2] 

71,74 [7, Fig. 11] 
a 2 ,a 3 [20, Fig. 1] 

(+-) 

Yes/No 


9 [40, Fig. 6.3] 

75 [7, Fig. 11] 

(+ +) 

No/No 

41+ 

Absent. 

(4 )> (-b) 

Yes/Yes 

2 A* 

6 [40, Fig. 6.3] 
a 4 [46, Fig. 2] 

72 [7, Fig. 11] 

<Si,to [20, Fig. 1] 

(- +) 

Yes/Yes 

B+ 

5 [40, Fig. 6.3] 
b 2 [46, Fig. 3] 

73 [7, Fig. 11] 

Pi [ 20 , Fig. 1 ] 

(4 —)> (4—) 

Yes/No 

2 

D [38, Fig. 2] 

76 [7, Fig. 11] 

(4- 4-), (+ +) 

Yes/No 

2A+ 

E [38, Fig. 2] 

7r [7, Fig. 11] 

(+ +)> (-) 

No/Yes 

+ 

o 5 [46, Fig. 2] 
a 5 ,a 6 [20, Fig. 1] 

(- -), (- -) 

No/Yes 

2A_ 

c 3 ,c 4 [46, Fig. 4] 
a 3 ,a 8 [20, Fig. 1] 

(-4), (-h) 

No/Yes 

2B+ 

64 [46, Fig. 3] 

[20, Fig. 1] 

(- -), (- -) 

No/Yes 

2A_ 

d 2 ,d 3 [46, Fig. 5] 

«3, a 8 , <+ 9 , 0+0 

[20, Fig. 1] 














































































and belongs to M\. Therefore, the maximal and minimal values of H for fixed P on G 4 coincide with 
these values on A4 1 : 

max H = max H = H 3 (p), min H = min H = HAp). 

G 4 n{P=p} MiC\{P=p} G 4 n{P=p} Min{P= P } 

In this regard, we introduce the additional parameters 

L 4 = h — Hi(p), L 5 = h-H 3 (p ) 

and the following notation for domains on the (p, h)-plane: 

al = {(p, h) : L a < 0}, aj = {(p, h) : L 5 > 0}. 

In particular, the facts proved above imply the following assertion on the existence of degenerate 
rank -1 critical motions in the first subsystem. 


Proposition 15. The admissible domain of the bifurcation diagram £ of the moment mapping J 
consists of the following segments of singular points on the surface IIi —the images of degenerate 
rank -1 critical motions: 

(1) the whole curve along which the surfaces IIi and II 3 touch: 


= ±\l-(l- 2 XTs), 


h = 


1 - A 2 s + 2 s 2 
2 s 


0 < s < 


2A 2 


( 2 ) the cusp of the surface IIi between its intersections with the curve £ 3 : 


Ai : £ = ± 


3\/3 


K) 


2 \ 3/2 


A 2 


<h< h c { A), 


or, in terms of the parameter s, 


,_ 3 A2 

h ~ 2 S + T’ 


\l\<tc{ A), 


s € [ 0 , sc (A)], 


(4.10) 


where the functions /ic(A) and CcW that determine the extrema of the corresponding coordinates 
on 5 3 for A > 0 and the value of the s-coordinate sc(A) at the extremum point are obtained from 
Eqs. (3.11) and (4.6) by using the auxiliary parameter: 

h C = z x + 


( .3.1 ; 4 - 

- 2 Ax 3 -4 = 0, 


(4 - x 4 ) 3 / 2 

I x G 

-^ 473 , 0 ) ^ ' 

£ c = 

4x 3 


s c = 


4 — x 4 
2 x 2 


The information in Table 4.2 supplements the results obtained in [14, 16] concerning the classifica¬ 
tion of motions on G 4 . 

The existence of of asymptotic motions or regular tori was proved in [14] by the analysis of explicit 
quadratures (1.34)-(1.35). Thus, we see that regular tori may be located on a critical level of the first 
integrals only in the neighborhood of type-A atoms, whereas asymptotic motions are, of course, form 
a part of hyperbolic atoms. The nonlocal distinction between pairs of points in the domains separated 
by the curve £q, i.e., points of the domains 010 and a 12 , consists of the fact that at the level of common 
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Table 4.2 


Domain 
(see [ 14 ]) 

Conditions 

Domain of A 4 i 

Periodic motions 

Asymptotic motions 

Regular tori 

I 

L\ > 0 , L 2 > 0 , 
L 3 > 0, -L5 > 0 

02 

No 

No 

No 

I 

L\ > 0, L 2 > 0, 
L 3 > 0, L§ < 0 

os (0 < A < A*), 

04, ag (0 < A < +00) 

“center” 

No 

Yes 

II, VI 

L\ > 0 , L 2 < 0 , 
L 3 < 0 

a 2 (0 < A < +00), 

010 (A > A*) 

“center” 

No 

Yes 

III 

L\ > 0 , L 2 < 0 , 
L 3 > 0, -L5 > 0 

02 

No 

No 

No 

III 

L\ > 0 , L 2 < 0 , 
L 3 > 0, L5 < 0 

03 

“saddle” 

Yes 

No 

IV 

L\ > 0 , L 2 > 0 , 
L 3 < 0 

o.5 (0 < A < 1 ), 
a 6 (0 < A < a/2), 
an (A* < A < +00) 

“saddle” 

Yes 

No 

V 

L\ < 0 , L 2 < 0 , 
L 3 < 0, L 4 > 0 

or (0 < A < +00), 

012 (A > A*) 

“center” 

No 

No 

V 

L\ < 0 , L 2 < 0 , 
L 3 < 0, L 4 < 0 

oT 

No 

No 

No 

VII 

l 2 = 0 

Ai 

degener. 

Yes 

No 

VIII 

0 

II 

CO 

Ao 

degener. 

Yes 

No 


first integrals for the domains 02 and an, regular tori exist, whereas for the domains a± and ai 2 , there 
are no such tori. This difference in terms of parameters of analytical solutions is found in [11], 

Due to (2.13) and (1.30), 

£ = -p^/i-y-p 2 ^ , s = h~Y~P 2 , 

it is easy to classify (S', i/)-diagrams and (S, L)-diagrams of the first critical subsystem. The separating 
values A remain the same (they are defined by surgeries of sections of the key set). In Fig. 12, we 
present the (S, i/)-diagr arris for the folloiwng cases: 

(a) 0 < A < A*; 

(b) A* < A < 1; 

(c) 1 < A < A*; 

(d) A* < A < \/2; 

(e) A > \/2; 

(f) A = 0. 

We indicate all points on the planes (p, h) and (s, h ) that are important for the construction of the 
//-atlases of bifurcation diagrams of systems on Qj t (points of the extreme values of /.-coordinates on 
the key sets). The corresponding notation is given in Figs. 13 and 14. Points associated with the 
curves <5i, 62 , and d '3 are denoted by the letters A, B, and C, respectively (with indexes if necessary). 
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Fig. 12 (continued). (S, ff)-diagrams of systems and their enlarged fragments 
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Fig. 13. Singular points of the (P, 4/)-diagram of the system A4 1 



Fig. 14. Singular points of the (S,H )-diagram of the system A4i 


The letter D is used for points on A j (here j = 0,1; in the third system, A 3 appears). To see all 
singular points, it suffices to consider the cases (a) 0 < A < A* and (b) A* < A < V 2 in these figures. 

In Fig. 15 we present the (S, L)-diagrams of the subsystem A4\ for nonseparating cases: 

(1) 0 < A < A*; 

( 2 ) A* < A < 1; 

(3) 1 < A < A*; 

(4) A* < A < V2; 

(5) A > x/ 2 ; 

( 6 ) A = 0. 

Critical points that are relevant to the construction of an L-atlas of bifurcation diagrams for systems 
on Pf are marked here. Note that, in comparison with the diagrams containing H, the additional 
point H 4 appears: it is an extremum of the ^-coordinate on the image of degenerate points Ao- 

We clearly list all values of the parameters and integrals (general and particular) at marked critical 
points. We obtain the image of degenerate rank-0 critical points cutting the separating set in Fig. 2 
at the set level A. The points on the axis r = 0 are considered only on the curve 82 - The curve 8 \ 
does not contain degenerate points but contains an extremum of h, namely, the global smallest value 
of energy reached at the point 

A 2 

A : r = 0, p = 0, s =-1 - —, h = -1, £ = 0. (4.11) 
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Fig. 15. (S, L)-diagrams of M\ subsystem (detailed elaboration) 
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Fig. 15 (continued). (S, L)-diagrams of A4 1 subsystem (detailed elaboration) 
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Introduce the notation 


Qo — A (r — A) + D, 


Qi = (r — A)(2r — X) — D, 


Q 2 = r(r - A) (2r - A) + A D, Q 3 = r(A — r) + ZA 

On the curve 5 - 2 (r < 0) we find 


dh 1 

dr 2(A — r) 2 -D 
dp 2 _ 1 

dr = _ 2(r — A) 2 D 


Q 1 Q 2 , 
Q 2 Q 3 , 


dr 2D^ 3 ’ 


1 


d7 2 


dr 8 (r — A) 2 D 


QoQiQ2<53- 


In particular, this implies 

dh 

dp 


2 r Qi 


d7 


r-A^QT ds 2v / 2r(r - A) 3 / 2 


Ql\/Q3- 


(4-12) 


Thus, the extrema of /i(r) on (52 are as follows: the point Ri (at r = 0), the degenerate point B 4 (the 
curve 7T23) 

(A-r)(A-2r)-D = 0, (4.13) 

and the cusp B 3 (the curve 7 ^ 4 ) 


r( A — r)(A — 2 r) + AD = 0. 


(4.14) 


The touching of the curves 5 2 and Ao generates the point B$ (the curve of degenerate points 7T2 i , 
Eq. (3.21)), and the intersection of the three curves <52, Ao, and Ai occurs at the point B 4 (the curve 
of degenerate points 7 ^ 2 , r = —A). In accordance with (3.11), the intersection of 5 2 with the curve of 
multiple points is determined by the equation 


\(r-X) + D = 0. (4.15) 

Solutions of Eq. (4.13) are written in the parametric form (3.37) by using the substitution (3.27); 
solutions of (4.14) can be found explicitly. Equation (4.15) can be reduced to the equation 

(r — A ) 3 (r + A) + 4 = 0 


under the condition r < 0. Using the same substitution (3.27), we obtain the necessary parametric 
representation of the coordinates of the points Bq and Bj for A > A*. In addition, according to (4.12), 
critical points of the curve d 2 on the plane ( s ,£) are generated by the condition Q 3 = 0 , which implies 
that p = 0 and 7 = 0, i.e., this is again the point A. Recall that, due to the identity £ = —sp , the whole 
curve 7o corresponding to the value s = 0 “collapses” on the plane (s,7) to the origin; in particular, 
the axis 0£ consists of inadmissible points except for (0,0). Therefore, in Fig. 15, it is shown by the 
dot line. Thus, the curve 5 2 generates the following set of points and values of the first integrals: 


A 2 1 

B 1 : r = 0, p = 0, si = 1 - —, S 2 = h=l, 7 = 0; 


B 2 = 7T22 : < 


r = —A, 


Si = Vl + A 4 — A 2 


” = ± \ 

S2 = Si 


' yiTA 1 - A 2 


h = | y/l + A 4 - A 2 , l = =f-^= (Vl + A 4 - A 2 ) 7 ; 


(4.16) 

(4.17) 
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B 3 = 7T24 : < 


r 2 


•si = 2 


4 + C / 2 - [/) , 
V4 + C / 2 - 2C7 


± i^( v ^ 5 - f, ) ;V2 - 


- U 


2 > 


P = 


S 2 = 


" 2\/2 


T 7 <t/=An 


2C7 


h = - , 
4 L 


(4 + U 2 ) 3/2 - U (6 + U 2 ) 


e = T- 


V4 + U 2 - 2U 


y/2 (\/4 + c / 2 - [/) 


1/2 ’ 


#4 — 7T23 : < 


x 4 - 4 


r = 


2 x 3 ’ 


p = ± 


v/4^ - 


2 x 


A = 


= 


h = -x 2 + 


x u 


= =F 


3x 4 - 4 
2 x 3 ’ 

4 — x 4 
2 x 2 ’ 

(4 — x 4 ) 3 / 2 


«2 = 


x € 


f r = A - A" 1 / 3 , 

1 


B 5 = 7T21 : < 


Si = 


2A 2 / 3 : 


4x 3 

p = ±\ l ^\/1 — A 4 / 3 , 
S2 = si, 


2(3x 4 - 4)' 

(</473,V2 


A G (0,1], 


A 2 


1 


#6,7 : < 


h ~ “Y + a2/3 + 2A 2 / 3 ’ 

12 1 




%/1 — A 4 / 3 
2A 1 / 3 ; 


r = —x + 


A — -x H—o > A*, 


x° 


p = ± 


\/^ 4 — 4 


x J 


si = 0, s 2 = 


2 x 2 


4 + x 4 ’ 


, 1 2 2 2 

n — ~ X —77 — —. 

8 x 2 x 6 


-Bo : x € 


\/ 2 , 


4 = 0 , 

B 7 : x € 


V 2 ^ 


+oo 


Here si and s 2 are the values of the parameter s on the surfaces Hi and ^ 3 , respectively. 
(4.18) at the cusp B 3 can be simplified by introducing the parameter z: 


B 3 — 7T24 : < 


Si = 


/1 \ 3/2 

- z\ 

0 < z < 1 , 

J 

-y/z( 1 - z 2 ), 

p = ±z 3/2 , 

3z 2 - 1 

z 

2 z 3 ’ 

82 - 2(1 -z 2 

1 , 9 , 

3z 2 - 1 

-z(3 + z ), 

l — zn 

1 2 z 3 / 2 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


Formulas 


(4.22) 
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On the curve £ 3 , the minimum point of h corresponding to the degeneration curve 7 r 3 i is determined 
by the same equation (4.13), and the parametrization of this point is given by Eqs. (3.38). We have 


r = 


x 4 — 4 
2x 3 ’ 


A = 


3x 4 — 4 
2 x 3 


vr 3 i : < 

V4 — x A 

4 — x A 

s 2 = 

X 6 

P T 2 x ’ 

Sl ~ 2 x 2 ’ 

2(3x 4 - 4)' 


, 3 2 2 

h = -x T —, 

1 8 x b 

(4 - x 4 ) 3 / 2 
^ 4x 3 ’ 

x G 

-^473,0) 


(4.23) 


Finally, the extremum values of h on the curves Ao and Ai are reached at the points 

1 


D 1 


D 2 


d 3 



S2 = Si, 

i = 0 , 


A > 0; 


A > 0; 


si = 


y/ 2 ’ 


\ 2 

h = V 2 --, 
2 ’ 


= +1 


'1 - V2X 2 
2 V 2 


A 2 < 


V 2 ' 


(4.24) 


(4.25) 


(4.26) 


In addition, we see from (4.10) that, as was mentioned above, there is an extremal value of t on Ao 
at the point (see Fig. 15) 

1 + 6A 4 1 1 

4A 2 ’ 4A’ Sl 4A 2 ’ A > 0. (4.27) 

S 2 = si, P = +A, 


D 4 : 


h = 


4.3. Specification. Second and third critical subsystems. For the second and third critical 
subsystems, we rewrite the equations of the surfaces ( 2 . 8 ) in the form 

n 2 U n 3 = j/i = 2 ^ 2 + i- - y(1 - 4s 2 ), k = -M 2 \ 2 + -^2 - y (1 - A 2 s)(1 - 4s 2 ) 

As above, s < 0 for II 2 and s > 0 for II 3 . This implies that the pair (s,£) determines a unique point on 
the corresponding surface and this correspondence is bijective. So, we can consider the (S, L)-diagrams 
for the subsystems M . 2 and M. 3 . 

To obtain a simple criterion for the existence of solutions to systems (1.27)—(1.29), we represent 
them in algebraic form. We perform the substitution 

1-C 2 


x = 


1 + C 2 


Y = 


2C 

1 + C 2 


We have single-valued dependences 


l 2,xp(, 

S 1 + C^ 


cj 3 — at 2 x 


1 -C 2 


Ol\ — 


Aa(l - C) + 2^C(1 + C 2 ) - 


.3 a2 


x(l + ( 2 ) 2 


a 3 = 


1TC 2 ’ 

£(1 - C 2 ) - 2 A ps( 
M1 + C 2 ) 
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and the formulas with radicals 


oj 2 — ~ 


y/2(l + C 2 ) V xs 


z(0, «2 = - 


2\/2xc rr 

(1 + C 2 ) 2 V XS 


z(a 


where 

Z(C) = (x - 2 A.s 2 )C 4 + 4£psC(l + C 2 ) T 2 x(l - 4 x 2 s)C 2 T (x + 2 As 2 ). 
The dynamics is governed by the equation 


dt 2y/2 V XS 


Further, due to (1.27), we set 

c = ^ 

C = iz 


P = P+, 
P = i P-, 


P 2 > 0; 

P 2 < 0 


(we assume that p + and p_ are nonnegative). We obtain the solutions in the following form. For 
p 2 > 0 , we have single-valued dependences 


£ 2xp + z 

^1 “j j 2~ ’ 

S 1 + Z z 

_ Aa(l - z 4 ) + 2£p+z(l + z 2 ) - 8 x 3 z 2 
ai “ x(l + z 2 ) 2 ! 


CJ 3 — A + 2>^' 


1 + ^ 2 ’ 


£(1 — z 2 ) — 2A p+sz 


x{l + z 2 ) 


and the formulas with radicals 


Similarly, for p 2 < 0 we obtain 

l 2xp_z 

w i =-- 2 " ’ 

s 1 — z z 


1 7 /.A 2 V 2 XZ 

„.„Z+\Z), / -1 oy 


(1 + z 2 ) 2 V >rs 


Z+(z). 


CU 3 — At 2^< 


1-z 2 ’ 


Aa(l - z 4 ) - 2£p_z(l - z 2 ) T 8 x 3 z 2 


£(1 T z 2 ) T 2A p-sz 


x(l — z 2 ) 2 


x(l — x 2 ) 


and the formulas with radicals 


oj 2 = — 


\/ 2 (l — z 2 ) 


1 V ^ _ 2\/2xz 


(1 — £ 2 ) 2 V xs 


Z„(z). 


Z + (z) = (x — 2Xs 2 )z 4 + 4£p + sz(l + z 2 ) + 2x(l — 4x 2 .s)2: 2 + (x + 2As 2 ), 

Z_(z) = — (x — 2 As 2 )z 4 T 4£p_sz(l + z 2 ) T 2x(l — 4x 2 .s)2 2 — {x T 2As 2 ). 

The dynamics is governed by the corresponding equation (the upper sign is for p 2 > 0 and the lower 
sign is for p 2 < 0 ): 


Z±(z). 


dt 2y/2 V xs 

From these expressions we obtain the following criterion. 


(4.28) 


Proposition 16. For given values s and £ corresponding to a level in M 2,3 that does not contain 
rank-0 critical points, the number of critical circles in the systems M. 2,3 is equal to the number of 
trajectories in the phase space 1x1 of the corresponding equation (4.28), where Z± is a straight 
line supplemented by the point Z±. Thus, if the corresponding polynomial Z± has 2m of real roots 
(rn = 0 , 1 , 2 ), then the critical levels s, £ contains m critical circles, except for the case where m = 0 
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and the leading coefficient of the polynomial is positive. In this case, there are two critical circles 
(z runs over all M), and the variable ui 2 preserves its sign on each of them. 

The last case (m = 0) is considered as a special case: 

p 2 > 0, x - 2\s 2 > 0, Z+(z) >0 Vz € M. (4.29) 

It is impossible to write out both Morse-Bott indexes explicitly for .M 2,3 systems. However, it is 
possible to obtain simple formulas for their calculation and, which is more important, to distinguish 
clearly a vector on the transversal area to the trajectory at hyperbolic points on which a gap of the 
“figure-eight” cannot occur (except for the special case (4.29)). This allows one to determine the 
direction of type-B atoms when the integral K increases. 

As above, we select the orthogonal complement to the vectors gradT, gradT, grad H, and sgrad IT 
as a transversal area for the critical circle (at a rank -1 critical point). 

We consider the case that does not satisfy (4.29). On any trajectory, the variable z oscillates 
between the roots of the corresponding polynomial Z(z), including infinity. Consider a point xq on 
the trajectory such that Z(z) = 0. Then the vectors gradT, gradT, and grad If are orthogonal to the 
plane Ouqa^ and the vector sgrad H lies in this plane and has the form 

sgrad H = (0, b 2 , 0, 0, 65 , 0). 

The constrained extremum of the function K at the common level of the functions T, T, and H in P 6 
is a critical point with undefined Lagrange multipliers 

o 1 \ 9 2 x 2 

2A 2 -] H + 2sL 2 - r. 

s) s 

Obviously, the part of this function that does not contain the Casimir functions T and T coincides 
with (4.3). Let £ = d 2 K 2 (x 0 ) be the matrix with elements Cij. Since the vector sgrad H belongs to 
the kernel of C core, we have 

£22^2 + £25^5 = 0 , £25^2 + £55^5 = 0 . 

In particular, since 62 ^ 5/0 and the coefficients £25 and £55 have the form 

r - 8 * pC r - on 2 

L25 — ~T+C 7 ’ 

we can perform some simplifying: 

= £25 = 32 x 2 C 2 = £25^2 = 4 x ( b 2 

^~' 22 £55 (1T C 2 ) 2! 5 “ £55 _ P(1 + C 2 )' 

As the first vector of the transversal area, we choose the vector lying in the same plane OLO 2012 an d 
orthogonal to sgrad H: 

vi = ( 0 , - 65 , 0 , 0 , b 2l 0 ). 

Then for computer calculations, the vector V 2 can be easily found as the null space of the matrix 
consisting of the five vectors: grad T, grad L , grad H , sgrad H , and v\. Although its analytic expression 
is quite complicated, all necessary calculations can be easily performed in the dynamic mode. In 
addition, the matrix of the restriction of the quadratic form C to the transversal area is diagonal in 
the basis (rq,^), so the Morse-Botta indexes have the following form: 

l-i 1 = Cv 1 • Vi, /l 2 = CV 2 ■ v 2 . 

The analysis of projections of integral manifolds to the plane Ouqw 2 shows that the critical surface of 
the hyperbolic circle never has discontinuities in the direction of the axis Ovj 2 in a neighborhood of 
the systems A 42 , 3 - 

In particular, this means that for type-P atoms, the vector v\ is always directed to the exterior 
domain of the “figure-eighth.” Therefore, if p,\ > 0, then the function K increases in the direction of 
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Fig. 16. Special case 


the “head” of the atom (“head” up), whereas if pi < 0, then the K function decreases in the direction 
of the “head of the atom (“head” down). From the expressions obtained we have 

9 9 [16x 2 c 2 + p 2 (i + C 2 ) 2 1 2 9 

Ml = £22^5 — 2 £ 25 & 2&5 + £. 55^2 = 2- 2/1 1 7214 -^2 • 

F (1 + W 

Thus, the direction of edges of asymmetric atoms is determined by the sign of p 2 . Since in Ad 2 all 
rank -1 points have elliptic type and p 2 > 0, we obtain the following statement. 

Proposition 17. If the integral K increases on an isoenergetic level Qj fl , then at points of the critical 
subsystem Ad 2 for each nondegenerate critical circle, the bifurcation A + of the torus birth occurs. 

In the system Ad 3 the situation is more complicated since we know from the classical problem that 
in this case CVtype atoms can appear (see [39]). For the nonsingular case, it was proved above that 
for two trajectories at the same integral level, the distribution of signs of the Morse-Botta indexes is 
the same. Let us consider the special case (4.29). On the plane ( s,i) plane, the domain 


AV( 2 s 2 - 1 ) < t 


< |(1 - 2A 2 s) 


corresponds to the positivity condition of these constants. In Fig. 16, these conditions are fulfilled 
in the domains 1-3. The curve kq corresponds to the boundary case where the leading coefficient 
of the polynomial Z + vanishes. According to the agreement adopted above, in all spaces of integral 
constant, the curve 82 is the image of the corresponding set of rank -0 critical points; certainly, it is 
the discriminant set of the polynomial Z + . The notation Ao and A 3 has the same sense as above. In 
the domains 2 and 3, the polynomial Z + has respectively two and four real roots. In the domain 1, 
there are no real roots, so the case (4.29) is realized. Obviously, the domain 1 has an intersection with 
the axis I = 0 at 

s € (0, min{l/ 2 ,1/2A 2 }); (4.30) 

since > 0 on the whole axis, it is convenient to take z = 0. Then one can easily find vectors 
determining the transversal area to the periodic solution: 

v / 2 s( 1 - 2A 2 s) 


Vi = —(1 + 2A 2 s), 0, 0, 0, 


VI + 2 s 


■ [(1 + s)A 2 s — l], As(3 — 2A 2 s) 


V2 = 0 , 


Xy/ s(l + 2 s) 
■v/2(l - 2A 2 s)' 


1 , 0 , 0 , 0 


The matrix /C of the quadratic form on this pair of vectors is a diagonal matrix whose eigenvalues are 
as follows: 


4s 


Mi = 


1 + 2 s 


[l + A 2 s 2 (5 — 2A 2 s)] 2 , 


1 


M 2 = -( 8 A s d - 1). 
s 
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It is clear that fi\ > 0. The vanishing of p 2 corresponds to the set A 3 where rank-1 points become 
degenerate (the vertical straight line in Fig. 16). On the interval (4.30), we always have P 2 < 0. Thus, 
in the special case, the Morse-Bott indexes have opposite signs (this fact is obvious since trajectories 
are hyperbolic) and, by the above, the distribution of signs is the same on both trajectories of this 
level, namely, (+, —) for chosen basis. We note that the calculation of the Morse-Bott indexes does 
not allow one to distinguish the atoms 2 B and C 2 by local analysis. 

We obtain the following proposition. 

Proposition 18. If the integral K increases on the isoenergetic level Q\ h , at points of the critical 
subsystem .M 3 on the nondegenerate critical circles we have the following bifurcations: 

(1) for elliptic trajectories (type “center”) we have the appearance of a torus for p 2 > 0 (the atom A + ) 
and the disappearance of a torus for p 2 < 0 (the atom T_); 

(2) for a hyperbolic trajectory at a critical level of K: the atom B_ for p 2 > 0 (an “outer” edge up 
and a couple of “inner” edges down ) and the atom B + for p 2 < 0 (an “outer” edge down and a 
pair of “inner” edges up)-, 

(3) for two hyperbolic trajectories at a critical level of K: two atoms B the direction of whose “outer” 
edges is defined by the same rule (both “heads” up for p 2 > 0 or both “heads” down for p 2 < 0)or 
the atom C 2 . 

Item 4 of Proposition 13 (i.e., for two hyperbolic trajectories at a critical level of K , we have different 
combinations of signs in the pairs (p ,\, ^ 2 )) is impossible for the subsystem .M 3 namely due to the fact 
that the sign of pi coincides with the sign of p 2 , i.e., is the same for all critical circles on the same 
level of integrals. 

Conditions of the existence of motion in the subsystems .M 2,3 were obtained in [48, 50] in terms of 
the integral constants (s, h) and (s,£), respectively. From Eqs. (1.28) and (1-29) it follows that 

sgn(p 2 ) = sgn(F 2 ), sgn(p 2 ) = sgn(sf/ 2 ). 

Let T* = Y for real p and Y* = iY for purely imaginary p. Then the curve 0 in the plane (X, Y*) 
defined by the equation X 2 + Y 2 = 1 (the necessity of this equation is implied by (1.27)) is a circle or 
a hyperbola, whereas the curve Ti defined by the equation Q 2 (X,Y) = 0 is an ellipse for all p 2 / 0. 
We obtain the following result (see [48]). 

Proposition 19. Real solutions of the system (1.27) for given s, h, and £ related by the equations of 
the surfaces 112,3 exist if and only if the following conditions hold: 

( 1 ) if s < 0 , then the circle To and the ellipse Ti have a common point ; 

( 2 ) if s > 0 and p 2 > 0 , then the circle To does not lie entirely in the interior of the domain bounded 
by the ellipse T 1 ; 

(3) if s > 0 and p 2 < 0, then the hyperbola To and the ellipse Ti have a common point. 

Indeed, one can immediately verify that for s < 0, the point X = 1, Y = 0 lies outside the ellipse. 
The motion on To is performed along a segment where Q 2 < 0, i.e., inside the ellipse. Therefore, 
segments of the circle inside the ellipse exist if and only if the intersection To n Ti is nonempty. If 
s > 0 and p 2 < 0 , for an actual motion we have Q 2 < 0 , i.e., at least one point of the hyperbola lies 
inside the ellipse, which again means that To FI Ti 7 ^ 0 . Finally, if s > 0 and p 2 > 0, then there exists 
a point on the circle at which Q 2 > 0, i.e., a point outside the ellipse. This does not hold only if the 
circle To lies entirely inside the ellipse Ti. 

Obviously, this assertion and Proposition 16 are consistent. 

Moreover, the number of critical circles at a fixed level of first integrals is equal to the number of 
segments of the circle To lying inside the ellipse Fi for s < 0 and to the number of segments of the 
hyperbola To for s > 0 and p 2 < 0. For s > 0 and p 2 > 0, the number of critical circles is equal to the 
number of segments of the circle inside the ellipse, except for the special case (4.29). In this case, the 
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whole circle lies outside the ellipse, i.e., there is no intersection points; there exists one segment but 
two trajectories due to the possibility of choosing the sign of the nonzero radical Q. 

Recall that D = y^4 + ?’ 2 (j— A ) 2 > 0. In addition to (4.5), we set 

0±{r) = ^j^[r{r- X)^D], 

V±(.r) = ^ [A(r - A) ± D] \/tp±(r). 

As follows from (3.11) and (3.15), the functions <p±(r), 9±{r ), and r)±(r ) are expressions of the integrals 
H, S , and L at rank-0 critical points that belong to the subsystems M 2 and M 3 . 

The following theorem corresponds to [50, Proposition 4], 

Theorem 11. The (S, L)-diagram of the critical system M 2 consists of the following sets: 

<5i : s = 6L(r), £ = ±r/-(r), r G [0, A); 

5 3 : s = 6 + (r ), £ = ±r?+(r), r € (A, +00). 

The outer boundary of the admissible domain T >2 is the connected curve S\. The curve 63 consists 
of two components and divides T >2 into three subdomains. Points of the subdomain containing the 
values £ = 0 correspond to a critical circle ; points of the other two subdomains bounded by the curve 63 
correspond to two critical circles. There are no surgeries of the diagrams with respect to the parameter 
except for the limit case A = 0. 

In Fig. 17, together with the (S', L)-diagram of the second critical subsystem (for the general case 
(a) A > 0 and for the limit case (b) A = 0), the domain b\ with one critical circle and two domains 62 
symmetric with respect to £ = 0 with two critical circles are shown. In the domain marked with an 
asterisk, To n Ti 7^ 0; this means that there is no motions here. According to Proposition 11, all 
rank-1 points have the “center” type. The key set here contains only rank-0 critical points. As usual, 
we denote by the same symbols the images of singular points generated by the extremal values of first 
integrals on the key set for the same points of the preimages in Mi and M 2 (in our case, these are 
the points A and C ). The properties of the corresponding atoms are collected in Table 4.3. 
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Table 4.3 


Domain (life 
time) 

Number 
of circles 

Morse-Bott 

indexes 

Exit to 

A = 0/ 

£ = 0 

Atom 

Analogs 

bi 

(0 < A < +00) 

1 

(+ +) 

Yes/Yes 


1 [40, Fig. 6.3] 
ai [46, Fig. 2] 
ai [7, Fig. 11] 
ai [20, Fig. 1] 

b 2 

(0 < A < +00) 

2 

(+ +)> (+ +) 

Yes/No 

2A + 

Transition III—>-VI 

[40, Fig. 6.1e] 
a 2 [7, Fig. 11] 


The results obtained for critical subsystems A4 1 and A4 2 imply the following simple description of 
the admissible domain T> in the space R 3 (£, h, k) of constants of general integrals, which is the image 
of the phase space P 5 under the moment mapping (2.6). We assume that A is given. 

Theorem 12. The admissible domain V = J(P 5 ) is a simply connected set whose exterior boundary 
consists of the images of the domains a\ and a 12 of the first critical subsystem and the domains b\ 
and b 2 of the second critical subsystem. 

Proof. We fix a point (£, h) for which Qf h 7^ 0. The set of such points is a simply connected domain 
on the plane 0£h bounded by the curve <5i below (see Fig. 6) and homeomorphic to the closed half¬ 
plane. All Qj h are compact and connected; therefore, K(Qj h ) is a segment retracting over hi to a 
point. In particular, T> is a reduced bundle of segments over the closed half-plane homeomorphic (but 
not diffeomorphic) to the closed half-space. 

We fix h and consider the compact invariant subset Qt = H~ l (h) in P 5 and the bifurcation diagram 
£ Li<{h ) of the mapping 

LxK\ Qt :Qi^R 2 . 

Obviously, £la'(^-) is the image under the moment mapping of the intersections of Qf with the critical 
subsystems A4 t . We consider the closed domain on the surface IIi bounded by the image of the band 
of the plane ( p , h) between hi and £q (i.e., the image of the domains ai and 012) and the image of the 
whole subsystem .M2 on II2 (i.e., the image of the domains b\ and 62). The image of the subsystem .M3 
contains points inside this domain (for example, for l = 0). In the space M 3 (£, h, k), the image of AI3 
is connected; therefore, if there exists a point in J(M3) outside this domain, then the transversal 
intersection of the image J/.M3 fl Q\) and the closure of the union of the curves ai U 012 is nonempty 
(it is easy to understand from the geometry of IIi and II3 that the common part, which is the tangent 
curve, does not go beyond the surface II3). This intersection contains a part of the curve $2 on the 
boundary between ai and 012; however, by the above, on a small domain transversal to All in the 
preimage of the corresponding points of the diagram of the system A4 1, all rank-0 and rank-1 critical 
points have elliptic type and hence going beyond the boundary of this domain is impossible. □ 

Different types of envelopes of intersections of the set £ C M 3 (£, h, k) with planes on which the 
coordinate h is constant are shown in Fig. 18: 

(a) -1 <h< A 2 /2; 

(b) A 2 /2 < h < h c ( A); 
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(c) h > h c { A) 

(without sections of the surface II3 lying entirely inside). For h = — 1, the section of the admissible 
domain contracts to a point. In the case (c), the image of the domain aj is shown by the dotted line; 
as was shown above, for points of this domain, there is no real motion. Thus, the admissible domain 
has the shape of the “bow” of a boat whose “deck” is smooth at its beginning and then an “edge” 
appears in the middle. 



Fig. 18. Intersection of the admissible domain with the planes h = const 



Fig. 19. The “bow” of the admissible domain 


Now we turn to the classification of domain in the image of the subsystem .M3. A complete 
description of the admissible domain on the planes of constants of the integrals (S, H) and (S', L) is 
given in [48, 50]. The following theorem corresponds to [50, Proposition 5]. 
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Theorem 13. The (S, L)-diagram of the critical system M .3 consists of the following sets: 
d 2 :s = 0 + (r), l = ±ig + (r), rG(-oo,0]; 


A ° :^ = ±y|(l 
As :S = 2 ^’ 


- 2 A 2 s), 0 <s< 

'i G M, 

< . p . > 2 A 2 / 3 - Vi + A 4 / 3 

, ~ + A 4 / 3 — A 2 / 3 ) 1 / 2 ’ 


A < A*, 
A > A*. 


The admissible domain V 3 does not contain the following components of the completion to the diagram 
in which critical motions do not exist: the domain adjacent to the axis s = 0 and bounded the branches 
of the curves Ao and 5 - 2 for all A; the domain bounded by the curve d 2 between two points of intersection 
with the axis £ = 0 (r^O) for A > A*. 

Surgeries of diagram types in the domain A > 0 occur for the following values of the parameter: 0, 
A*, 1, A*, and \/2. 


Inequalities for r, s, and £ that determine the admissible domain for points of the key set are 
immediately implied by the criterion stated in Proposition 16. 

For critical points of the subsystems .M 2 , 3 , we separately formulate an assertion on the existence of 
degenerate periodic solutions, which is similar Proposition 15. 


Proposition 20. The admissible domain of the bifurcation diagram X of the moment mapping J 
contains the segment of the cusp of the surface II 3 , which is the image of degenerate rank-1 critical 
motion: 

[/i*,+oo) for A < A*, 


where 


Aq 


h* = - 77 —(3 - A 4 / 3 ), 


e = ± 


h — h* 


h G 


[h**,+oo) for A>A*, 


If* = - 


(4 + A 4/3 ) 3/2 - A 2/3 (6 + A 4/s ) = ^ (\/4 + A 4 /3-A 2/3 ) 2 (2v / 4 + A 4 / 3 + A 2/3 ) , 


or, in terms of the constants s and l, 

1 


Aq : 


s = 


2A 1 / 3 


G 


> 


2A 2 / 3 - ^4 + A 4 / 3 
\/2(\/4 + A 4 / 3 — A 2 / 3 ) 1/2 ' 


A < A*, 
A > A*. 


(4.31) 


The assertion in the coordinates (s,£) is a part of Theorem 13. The transition to the plane (£, h ) is 
performed by the formulas (2.8), which imply that on the surface II 3 

1 A 2 

h = 2 £~ + /io(s), ho(s) = — —-—|- 2A 2 s 2 . (4.32) 

It is easy to see that the minimal value of ho(s) is attained at s = 1/(2A 1//3 ), i.e., at the intersection 
point of the minimum curve h on M 3 with the image of the cusp A 3 ; this value is 


min ho(s) = h*. 

s> 0 

Then for A < A* and £ G I, we have h G [/i*,+oo). The boundary point is a new singular point 
on the key sets, which is important when considering 17-atlases, since it is an extremum point of the 
/i-coordinate on the image A 3 of the family of rank -1 degenerate points: 

1 A 2 / 3 

D 5 :£ = 0, *2 = ^, h = ^-(3-A 4 / 3 ), A G [0, A*]. (4.33) 
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(c ) l<^<(4/3) 3 4 


Fig. 20. The (S', L)-diagrams of the subsystem M.%: detailing 
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Fig. 20 (continued). The (S, L)-diagrams of the subsystem .M 3 : detailing 
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In the case where A > A*, calculating the value (4.32) on A 3 for the minimal admissible value of \i\ 
from (4.31), we arrive at the value h**. These values correspond to the singular point B% of the key 
sets obtained earlier. 

In Fig. 20, the (S, L)-diagrams of the critical subsystem .M 2,3 are shown: 

(a) 0 < A < A*; 

(b) A* < A < 1; 

(c) 1 < A < A*; 

(d) A* < A < \/2; 

(e) A > \/2; 

(f) the limit case A = 0. 

Here we should recall the notation (3.51). The dotted lines represent the degeneracy curves Ao 
and A 3 , except for a part of the curve Ao, which is the exterior boundary of the admissible domain T >3 
(the solid line). As above, in domains marked with an asterisk (here cT and C 2 ), critical motions are 
impossible. The diagrams also contain the detailing of key sets in accordance with the classification 
of rank 0 critical points given above and also with a clear decomposition of the smooth curves Ao 
and A 3 into smooth sections (for Ao, the decomposition has already been obtained earlier in the 
subsystem .Mi). Recall that this decomposition is generated by common points of these curves with 
the image of rank-0 points and intersection points. A relevant analysis will be performed in the next 
section. 

The assertions of Theorem 13 that regard curves involved in the diagram are consequences of the 
definition of the diagram. Surgeries of the diagram with respect to A are calculated directly (for detail, 
see [50]) or by the general property: they correspond to separating values of A for sections A = const 
of the family of separating curves for rank-0 points in the plane (r, A) and the set of extrema t on 
the images of rank-1 degenerate points. To prove the part of the theorem regarding the existence 
of motions, it suffices to analyze the distribution of roots and signs of leading coefficients of the 
corresponding polynomials Z± and apply Proposition 16. A computer visualization of the curves To 
and T 1 yields a visual proof in accordance with Proposition 19. The fact that in the domain cf adjacent 
to the axis Ol (see Fig. 20(a)) motions are impossible is almost obvious since s / 0. The absence 
of critical motions in the domain C 2 (see Fig. ??(d), (e)) was proved in [29] by the analysis of points 
lying on the axis t = 0 . 

In Fig. 20, we also show singular points appearing on the (S, L)-diagrams of the subsystem .M 3 . 
We can immediately verify that all of them have already been seen in the subsystem .Mi, except for 
the point D 5 defined by Eqs. (4.33) for the use in the 17-atlas. 

We reformulate Theorem 13 in terms of the integrals S and H (see [48]). 

Theorem 14. The (S, H)-diagram of the critical system M .3 consists of the following curves: 




htan{s), 

1 


Ao : 

h = 

0< *<^2 

5 

5 2 : 

h = 

<P+(r), 

s = 0+(r), 

r G (— 00 ,0] 



1 

jh> h*, 

A < A*, 

CO 

< 

s = 

2A 2 / 3 ’ 

j/i > h**, 

A > A*; 

h m [ n 

h = 

ho{s), 

s € /(A), 



where 

. . 1 — A 2 s + 2s 2 

«tan (S) = -X-- 
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and the range of s on the curve h m ; n is defined as follows: 

' (0, +oo), 

/(A) = < (0, s 0 ] U [s°, +oo), 

. (0, s 0 ] U [1/2,+oo), 


A < A*, 

A* < A < V2, 
A > V2. 


Here so(A) and s°(A) are abscissas of the tangency points of the curves 5-2 and h m , n that exist if A > A*. 
They satisfy the inequalities so(A) < s°(A) for A > A* and s°(A) < 1/2 for X* < A < y/2. 

The exterior boundaries of the admissible domains are as follows: 

(1) the curve h m i n within the limits s€/(A); 

( 2 ) the curve Aq within the limits 


s € 


(0,1/(2A 2//3 )], 

(0, y/l + A 4 — A 2 ], 


A < A*, 
A > A*; 


(3) the curve §2 within the limits 


[so, A 

[so, 1/2], 


A* < A < \/2, 
A > a/2. 


Surgeries of diagram types in the domain A > 0 occur at the following values of the parameter: 0, A*, 
1, A*, and y/2. 



Fig. 21. The connection of s and A at the points Bqj 


The curve h m ; n corresponds to the value t = 0, and its occurrence is naturally associated with the 
fold that appears on the covering of the plane (s, h) by the surface II3, since h is expressed through l 2 . 
A complete analysis of {S, i7)-diagrarris of the subsystem M 3 and its computer visualization was 
performed in [31, 48]. The values so and s° correspond to singular points Bq and By marked on the 
(S, L)-diagrams in Figs. 20(d) and (e). These points exist for A > A* (note that they coincide when 
A = A*) and at the moment A = y/2, the points B 7 and B\ merge and then pass into the inadmissible 
domain (the corresponding value of r would be positive on 82 , which is not so). For proofs and 
explanations, see [48]; they easily follow from the constraint equation for s and A at these points. The 
constraint equation can be obtained by eliminating the auxiliary parameter x from the corresponding 
parametric representation (4.21) of the values S 2 and A at these points (see Fig. 21): 


A 


vT + 2s - VI - 2s 
(2s) 3 / 2 ’ 


s € (0, -]. 
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Table 4.4 


Domain (life time) 

Number 
of circles 

Morse-Bott 

indexes 

Exit on 

A = 0/t = 0 

Atom 

Analogs 

Ci (0 < A < + 00 ) 

1 

(- -) 

No/No 

A_ 

Absent 

C 2 (0 < A < + 00 ) 

1 

(- +) 

No/Yes 

B+ 

a 5 [46, Fig. 2] 

03 [20, Fig. 1] 

c 3 (A* < A < + 00 ) 

2 

(- -), (- -) 

No/Yes 

2 A_ 

64 [46, Fig. 3] 
a 7 [20, Fig. 1] 

C 4 (0 < A < + 00 ) 

2 

(d - )> (d ) 

Yes/Yes 

c 2 

8 [40, Fig. 6.3] 

04,65 [46, Fig. 2, 3] 
02 [7, Fig. 11] 

7 [20, Fig. 1] 

c 5 (A* < A < V2) 

2 

( b)> ( b) 

No/Yes 

2 B+ 

63 [46, Fig. 3] 

04. [20, Fig. 1] 

c 6 (0 < A < 1) 

1 

(+ +) 

No/Yes 

A + 

o 3 , 0,4 [46, Fig. 2] 
a 4 [20, Fig. 1] 

c 7 (0 < A < 1) 

1 

(+-) 

Yes/Yes 

B_ 

7 [40, Fig. 6.3] 
a 3 [46, Fig. 2] 

01 [7, Fig. 11] 

02 [20, Fig. 1] 

c 8 (0 < A < A*) 

2 

(+ +)> (+ +) 

No/No 

2A+ 

Absent 

c 9 (0 < A < A*) 

2 

(d -), (d -) 

Yes/No 

2 

E [38, Fig. 2] 

03 [7, Fig. 11] 


Application of the results obtained to points of the domains cj-cg in the ( s , f)-image of the subsystem 
M 3 leads to the description of characteristics and atoms collected in Table 4.4. As we can see, all 
domains except for c\ and cs, when considering the expanded diagram in the space (s, £, A), are related 
with the corresponding domains of problems studied earlier (A = 0 or £ = 0), so only the atom direction 
is added here. In particular, the presence of the atom C 2 in the domain C 4 and the two atoms B in 
the domain C 9 is explained in [38-40](in another notation). The presense of the two atoms B in the 
domain C 5 follows from the results of [46]. In the new domain ci and cs, according to what was proved 
above, critical circles have elliptical type, their number is calculated according to the criteria proposed, 
and the direction is defined by the Morse-Bott indexes. 

4.4. Classes of rank-1 degenerate points. 

4-4-1. The set Ao- The first class of rank-1 degenerate points consists of points in the preimage of 
the tangent line of bifurcation surfaces. The image of this set in the space of integral constants is 
denoted by Ao; this is the tangent line of the surfaces IIi and II 3 . The equation of this line is 

2.s 2 -2f/i + yjs + l = 0 (4.34) 
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in the parameters s and h and 


2AV - s + 2 £ 2 = 0 (4.35) 

in the parameters s and £. Together with the realness condition, we obtain the following representation 
for Ao: 

£ = ± \l^ 1 ~ 2X2s )’ h = s + Ys~ Y’ sG ( 0 ’^A2_- 

The intersection of this set with the image of the set of rank-0 critical points (i.e., with the subsets eh, 
i = 1,2,3, of the space of integral constants) divides it into several parts with completely different 
behavior. Indeed, outside this intersection, a neighborhood of a point of the set Ao is mapped by a 
diffeomorphism as a stratified manifold and the topology of the corresponding preimage is preserved. 
Let us recall the equations of the sets <5*: 




A(r — A) + d 


—r + 


r — A 


h = —-r(r - A) + \- (L 

2 v ; 2(r — A) (4.36) 


d 2 = 4 + r 2 (r — A) 2 , r € (—oo, 0] U [0, A) U (A, +oo). 


The correspondence of the numbers of subsets, the intervals of values of the constant r that defines 
the relative equilibrium, and the signs of d are as follows: 


Si : 

re [0, A), 

d < 0; 

S 2 : 

r G (-oo,0], 

d > 0; 

As : 

r G (A, + 00 ), 

d > 0. 


The parameter s of the surfaces II, , which is a particular integral for the critical subsystems A4 t , takes 
the following values at the points (4.36): in the subsystem (for all three curves <5,;) 

1 


S 2 


A(r — A) + d 


and in the subsystems M 2 (the curves hj and 5%) and M 3 (the curve 62 ) 

r — A 


s = 


4A 


r (r — A) — d 


(4.37) 


(4.38) 


Now we use the expression on M 1; then at surgery points of Ao, the equation (4.36) must be com¬ 
patible, for example, with (4.34). We obtain 

2 


r z - 2 X Z + 


Hence either 


or 


A (r — A) + d 
r = —A 


r — A 


d = 0. 


(4.39) 


2 + (r — A) 2 (r 2 — 2A r + 2A 2 ) — (r — A) (r — 2A)d = 0. 
The last equation can be rewritten in the form 


1 

2 L 


(r — A )(r — 2A) — d 


12 


= 0. 


(4.40) 


In particular, this means that it described the tangency point of Ao with one of the curves <5*. We 
examine which curves eh correspond to the values (4.39) and (4.40). For (4.39), it is obvious that 
r < 0, so that it is an intersection point of Ao with 62 - A consequence of Eq. (4.40) is 

A(r - A) 3 + 1 = 0, 


which implies 


r = A 


1 

AW 


(4.41) 
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At the same time, according to (4.40), the equation 

sgn d = sgn(r — A )(r — 2 A) 

holds, which, under the obvious inequality r < A, is possible only for r < 0, i.e., on the curve r < 0. 
Thus, the tangency point (4.41) of Ao and 62 exists only when A < 1. 

Thus, taking s and A as independent parameters on Ao (if we replace s by h or by 7, then the 
correspondence in the domain is not bijective), we obtain the following separating curves in the 
square {(A, s) : A > 0, s > 0}: 

(1) the upper boundary of admissible values corresponding to the point D\ in the key sets: 



(2) the curve obtained from (4.36) and (4.39) and corresponding to the point B§ in the key sets: 

if>\ : s = \/l + A 4 — A 2 , A > 0; 

(3) the curve obtained from (4.37) and (4.41) and corresponding to the point in the key sets: 

^ = -=2^73. A 6(0,1]. 

The corresponding domain is shown in Fig. 22, where we denote the appearing subdomains in accor¬ 
dance with Fig. 15 and 20, where each such subdomain generates a segment of the key set. These 
subdomains play the role of nodal points in the bifurcation diagrams T,hk {£, A). Note that the sep¬ 
arating curves for Ao correspond to structurally unstable diagrams Hhk(£, A): on the plane (A, 7) 
the calculation of 7(A) by the formula (4.35) maps the curve ipo into the axis 7 = 0, the curve (p\ 
into the curve 1122 of the separating set ©£,, and the curve ip 3 into the curve 7T2 i in accordance with 
formulas (3.45). 

As we see in Fig. 22, the points A02 and A03 are related with the classical Kovalevskaya problem 
(A = 0), whereas the points A01 and A04 also exist in the Ryabov-Morozov case (7 = 0). So, at first 
glance, it may seem that the description of the relevant loop molecules can be performed within the 
frameworks of [7, 20]. However, this is not quite true. The topology of all levels of first integrals 
in punctured neighborhoods of these points was examined in [39, 46], but this does not answer the 
following question: How must one join different families while traversing the test points? Depending 
on this, the manifold in the preimage of a small circle centered at the point can have various number 
of connected components. Let us denote by Jo« the union of connected components of the preimage of 
point A 0 i that contain critical points of the moment mapping (i.e., critical circles consisting of rank-1 
degenerate points). 

Note that segments of Aoi inside the critical subsystem A4i are not bifurcation segments: they 
correspond to the degeneration of critical motions only relatively the complete system. In particular, 
the number of critical circles in the preimage Ao* is the same as in the preimage of the neighboring 
domains on both sides of a,j. According to the diagrams of the subsystem M \ and Table 4.1, we 
establish that Aoi separates 05 and ag (two critical circles in J01), A02 separates ag and an (one 
critical circle), A03 separates as and as (two critical circles), and A04 separates ag and an (two critical 
circles). Thus, the surface J02 is connected. For the rest of the points on the level Jo*, there are two 
critical circles, so one cannot restore loop molecules by the available data. In [7, 20], the hypothesis was 
implicitly used that the loop molecules of these points have two connected components, but there are 
no motivation for this hypothesis, not to mention proofs. In fact, a “principle of maximum likelihood” 
was used, according to which molecules to be found are among the already known molecules. We 
prove that the answer obtained is correct. 

Theorem 15. The surfaces Jqi, J 03 , and J 04 consist of two connected components. 


310 




Fig. 22. Separating curves for points of Aq 


Proof. For A 7^ 0, we use the quadratures obtained in [11]. The set Ao corresponds to the following 
values of the parameters: 

L 3 = 0, Li = | > 0. (4.42) 

O 

Therefore, the system is reduced to two equations 

V 2 = i (4.43) 

where 

F(z) = z(z - 2A) - ^ , f{z) = (| - z 2 ) F(z), 

the variable 2: is real, and the auxiliary variable 77, which is not equal to infinity only on asymptotic 
motions, is introduced by the equality 

V= vwr 

where the variable y is also real. From (4.42) and (4.43) we conclude that rj 2 > 0, so the variable r] 
is also real, and, in addition to the condition f(z) > 0 of the existence of a real interval where the 
variable 2: oscillates, we obtain the condition F(z) > 0 of the existence of asymptotic motions. Setting 


we obtain 


l = 0 , 


s = 


1 

2A2’ 


F(z) 




m 



F(z). 
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The variable z oscillates on the intervals 


1 

2A’ 


A 


1 

2A 



' 1 1 ' 

5 

2A ’ A+ 2A 


under the condition A < 1 and on the intervals 


1 

1 ' 


2A ’ 

2 A 

1 


1 

2A’ 



under the condition A > 1. In both cases, F(z) > 0 in the first interval and F(z) < 0 in the 
second interval. Therefore, asymptotic motions exist for the first critical circle and do not exist for 
the second. Consequently, the second critical circle exhausts its connected component of the critical 
integral surface. Thus, for points Aoi and A 04 , the preimage contains two connected components with 
critical circles. 

For the point A 03 we can set A = 0. Then 


n2> = z 2 -I(s-b), /(*) = (| - T) F(z), s> 1. 

Thus, z oscillates in the symmetric intervals 



and for both critical circles, asymptotic motions exist. Clearly, different intervals of oscillation of 
z cannot have the same limit when they tend to distinct critical circles. However, this still must 
be strictly proved based on quadratures obtained in [ 11 ]; this requires certain technical calculations. 
As a result, we obtain that on a connected component of the integral surface, the variable z cannot 
change its interval of oscillation, so there are two such components (the same number as the number 
of critical circles). We also present another proof of the existence of two connected components. Note 
that, despite the strong degeneracy, the formulas derived from the equations of first integrals do not 
provide any foreseeable solution to this problem. We use equations obtained by Kovalevskaya and the 
results of [42]. Recall that in the Kovalevskaya separation variables si and s 2 the equations of motion 
have the form 

(s2~ si)^-=iy/2S{sij, (-2-«i)^ = -i>/2Stej, 

where 

S(x) = (x — h + Vk)(x — h — Vk)ip(x), cp(x) = x(x — h) 2 + (1 — k)x — 2 £ 2 , 

and .si > s 2 . S. Kovalevskaya also found expressions for all phase variables through separation variables 
in the form of single-valued functions of si and S 2 and a set of algebraic radicals: 


Ri'f — \f^H 67 ? i — 1 ) 2 , 'y — 1 ,...., 5 , ( 4 . 44 ) 

where e 7 are roots of the polynomial S(x). We denote by e 2 , and e 3 the roots of ip(x) in the 
ascending order if they are real and set 64,5 = h \fh. It is known (see [4, 17]) that in the case where 
<p(x) has three real roots, the separation variables vary in the domains 


si € 


maxjei, e 2 , e 3 , e 4 }, e 5 


s 2 € 


- 00 , min {d, e 2 , e 3 , e 4 , e 5 } 


and, as was shown in [42] by the method of Boolean functions, the expressions for phase variables 
in the general case can be reduced to single-valued functions of the following products of algebraic 
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radicals in the real-valued form: 


Vsi - e\ 




e2 ~ ^2 
es — s 2 




e3 — s 2 
es — «2 




64 ~ ^2 

es — S2 


Vsi - e 2 




ei - s 2 
es — s 2 




e3 ~ ^2 
es — «2 




e4 — s 2 

es — ^2 ’ 


(4.45) 


a/si - e 3 




e2 ~ s 2 
es — ^2 




e2 ~ ^2 

65 - S 2 




e4 — s 2 
e5 - S2 


(the radicals containing quotients of differences do not change their sign for S2 = —oo). In the case 
where the separation variable is reflected from the boundary of its oscillation domain, the correspond¬ 
ing radical should be considered to change its sign (the radical corresponding to the value s 2 = —oo 
has already been excluded by the procedure of reduction to formulas (4.45)). If any of the boundary 
values is a multiple root, then the corresponding radical is also considered to change sign, since the 
asymptotic trajectory together with its boundary points is a connected set. 

Now we consider the point A03 for A = 0. In this case, 


1 1 

ei = - < e 2 = e 3 = e 4 = s < e 5 = s + - 
s s 

(recall that s > 1). Then, denoting the triple root by e*, we obtain 


si€[e*,e 5 ], s 2 € [-oo,ei]. 

Among the expressions (4.45), the second and third ones coincide, amd the product of the pair 
of identical radicals becomes a single-valued function that does not influence the sign of algebraic 
expressions. The following algebraic expressions remain: 


Vsi - ei 




e* — s 2 
es — s 2 ’ 


y/ S\ — e* 




ei - S2 
es - «2 


the second one changes its sign on the connected component of the integral surface and hence the 
number of connected components is independent of it. The first expression along the closure of any 
trajectory preserves its sign selected at the initial moment. Consequently, the surface J03 has two 
connected components. The theorem is proved. □ 


4 . 4 . 2 . Ai set. The following class of rank-1 degenerate points consists of points lying in the preimage 
of the cusp of the surface IIi. For integral parameters on this set Ai, taking into account the conditions 
of existence of motions, we have 

£ = ±—^=s 3 ^ 2 , ^ = 2 S ~^~ 2 ' (4.46) 

where s* = s*(A) is the largest real root of the equation 

9s 4 + 2A 2 s 3 - 24s 2 - 24 A 2 s + 4(4 - A 4 ) = 0, (4.47) 

which exist and is positive for all A > 0. The intersection of the set Ai with the image of the set of 
rank-0 critical points (i.e., with the subsets <5* of the space of integral constants) and with the set Ao 
examined above divides the set Ai into several parts with completely different behavior. Obviously, 
the curve <5] is not involved in these intersections. We can immediately verify that all points of the 
intersection Ai n Ao are contained in Ai 0 S 2 - The condition of intersection Ai n ( 5 2 U 63 ) can be 
represented either by Eq. (4.47) or by the equation 

s 2 + 2A 2 s -1 = 0. (4.48) 
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On the plane (A, s) we obtain the separating curves 


Vh : 


V2 : 


A 

A 


V / ?+4-2 s 

( v ^-,) 1/2 ’ S> 5 

\j\ G - ")’ se(o,i] ° 


s = \A + A 4 — A 2 , A > 0 . 


(4.49) 

(4.50) 


Taking into account (4.46), we see that the image of the curve ipi on the plane (£, A) is the union of 
the curves 7^3 and 7731 , which can be represented by the graph of the single-valued function 


s = (V 2 £) 2/3 , £ > 0 . 

The image of the curve V 2 is the curve 7^2 in the representation (3.45). 

The upper boundary s* of values of s is the branch of the curve (4.49) for s > 2/y/2> obtained as 
the intersection of Ai with 63 (the point C on the key sets). Another branch of the curve (4.49) for 
values s € (0, 2/\/3] and the curve (4.50) is the intersections of Ai with 62 (respectively, the points £>4 
and B 2 on the key sets; note that and C merge for A = 0). The separating set and the notation of 
the corresponding points on the diagrams of the subsystem A4 1 and on the bifurcation diagrams are 
presented in Fig. 23. The points An and A 13 are related with the particular case A = 0, the point A 14 
exists in the particular case £ = 0 , and the point A 12 does not have such analogs. 


Vs 2 + 4 — 2s 
(Vs 2 + 4 - s) 1 ^ 
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As above, the number of critical circles in the preimage of An can be found by the diagrams of the 
subsystem M\ and data for adjacent domains from Table 4.1. We see that An separates 0,2 and 06 
(one critical circle in the preimage), A 12 separates 03 and 04 (one critical circle), A 13 separates 07 
and as (two critical circles), and A 14 separates aio and an (two critical circles). In the case where 
the preimage of the point An contains only one degenerate circle, i.e., the corresponding critical 
integral surface and the loop molecule are connected, we establish below the number of sets involved 
in bifurcations for accurate description of the topology of neighborhoods of these points. This easily 
follows from the results of adjacent chambers and numbers of families in these chambers. If there 
are two circles, then for justifying the structure of loop molecules, it is necessary to determine the 
number of connected components of the critical surface, which is the preimage of An- By default, 
in [7] it was assumed that in the case where A = 0, the loop molecule of the point A 13 consists of two 
connected components, and also it was considered in [ 20 ] that the loop molecule of the point A 14 also 
consists of two components for t = 0. We prove the following general assertion using quadratures found 
by Gashenenko. Similarly to the previous case, we denote by Jn the union of connected components 
of the preimage of the point Aij that contain critical points of the moment mapping. 

Theorem 16 . The surfaces J13 and J14 consist of two connected components. 

Proof. By the definition of point classes, we assume that (A, s) belongs to the corresponding open 
subdomain in Fig. 23. In terms of (A, s), the polynomial f(z) from (1.32) becomes 

f(z ) = - A ) 2 + \ sz ( z ~ 2A ) - \{s 2 ~ !)• 

Obviously, its discriminant is the product of the polynomials on the left-hand sides of (4.47) and (4.48). 
The polynomial f(z) has four real roots in the domains A 13 and A 14 (it suffices to consider its traces 
on the coordinate axes when the corresponding equation is solved explicitly). As was noted above, the 
variable z cannot change its oscillation interval on the connected component of the integral manifold; 
therefore, there are exactly two such components. □ 

In general, according to (4.8), the points Ai correspond to the case L 2 = 0 in the Gashenenko 
solution. Then from (1.33) and (1.34) we obtain 

L\ > 0, F(z) = (z — A ) 2 > 0; 
therefore, for both critical circles there exist asymptotic motions. 


4 .4-.3. The set A 3 . The last class A 3 of rank-1 degenerate points consists of points lying in the 
preimage of the cusp of the surface II 3 . For the integral parameters on this set, taking into account 
the conditions of existence of motions, we have 


h = h* + 2 £ 2 , 


fel, 

\i\ > t, 


A < A*, 
A > A*, 


1 

* 2A 2 / 3 ’ 


k = -4A 2 £ 2 + k*, 


where 


A* = 


2y/2 


t = 


2A 2 / 3 - \J 4 + A 4 /3 


> 0 (A > A*), 


3 3/4 ’ x /2(v'4 + A 4 /3-A 2 / 3 ) 1/2 

h* = ^A 2/3 (3 - A 4/3 ) , k* = (A 4/3 - l) 3 + 1. 


The intersection of the set A 3 with the image of the set of rank-0 critical points (i.e., with the 
subsets Si of the space of integral constants) and with the sets Ao and Ai studied above divides this 
set into several parts with completely different behavior. Obviously, the curves <5i and A 3 are not 
involved in these intersections since they have no common points with II 3 . We can directly verify that 
all points of intersections A 3 0 Ao and A 3 n Ai are contained in A 3 n 62 - Since s is constant on A 3 
for fixed A, as the parameters defining the point A 3 we can take A and l. Then the conditions of the 
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intersection A3 D 62 (the point B3 of the key sets for all A and the point B$ for A £ [ 0 , 1 ]) are reduced 
to the equations of two separating curves 6\ and 62 and allow one to classify the diagrams Y<hk(£, A): 


0 \ ■ £=^Vl~\ 4/3 , 0 < A < 1 , 

* 2 : // A ^ 0 

>/ 2 (V 4 + A 4 / 3 _ AVs) 1 / 2 

Clearly, the equality £(A*) = 0 holds on 62 and 

2*/»-v/4T*g =f 

V 2(\/4 + W 5 -A 2 / 3 ) 1/2 

for A > A*. Therefore, the branch of the curve O2 for A > A* is the boundary of the domain A3. 
Separating curves and the notation of points A3* corresponding to Fig. 20 are shown in Fig. 24 . 
Comparing it with ( 3 . 45 ), we see that the curve 9 \ is exactly 7T2 i and the curve O2 is exactly 7^4. 






Fig. 24 . Separating curves for points of A3 


The number of critical circles in the preimage of A3j can be found by the diagrams of the subsys¬ 
tem 7W3 and the data for neighboring domains from Table 4 . 4 . We see that A31 separates ci and C2 
(one critical circle in the preimage), A32 separates cs and eg (two critical circles), A33 separates cq 
and C7 (one critical circle), and A34 separates C3 and C5 (two critical circles). 

As for the sets Ao and Ai, if the preimage of the point A3, contains only one degenerate circle, then 
the corresponding critical integral surface J%i is connected. We prove the following assertion, which 
will be useful in the sequel. 
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Theorem 17. The critical manifold J 32 consists of two connected components. 

Proof. As is seen from Fig. 24, the points of the class A 32 are related to the classical Kovalevskaya 
problem A = 0. We pass to the limit in the image in the parametric equations of IT 3 , assuming 
s = 1/(2A 2 / 3 ) —>■ 00 . As a result, h = 2£ 2 , k = 0, and the points A 32 correspond to the values £ > 1. 
The critical limit manifold is known (see [40]), namely, J32(0) = 2S' 1 . Therefore, for small A, this 
manifold is disconnected, and since there are two critical circles in the preimage, the number of 
connected component is also two. □ 

Note that points of the class A 34 are related to the case £ = 0. Loop molecules for this case were 
examined in [ 20 , 21 ], but the existence of two connected components in molecules for these points was 
assumed without proof. Later, we will prove this hypothesis. This implies that the critical manifold 
J 34 consists of two connected components. 


5. Topology of Reduced Systems 

5.1. Separating sets and bifurcation diagrams. In accordance with Proposition 12, the set Ql 
on the plane of the parameters (X,£) that classifies the bifurcation diagrams Yihk(£,X) of reduced 
systems on Pf consists of pairs (A,£), where £ is the critical value of the restriction of the function L 
to one of the key sets /Q of the critical subsystems A4i, i = 1,2,3, for fixed A. Obviously, this is the 
value of £ at the points (4.11), (4.16)-(4.21), and (4.23)-(4.27). As was noted above, if X£ 0, then 
only the point (4.27) does not correspond to degenerate rank-0 critical points. We recall that, due to 
existing obvious symmetries with respect to the parameters A and £, this set ©l is considered in the 
first quadrant: 

A > 0, £ > 0 ; 

moreover, the semi-axes A = 0 and £ = 0 are included into the separating set by default. The points 
A, B 1 , D 1 , and D 2 give the values £ = 0, A G R, and for A = 0 and any other £ there are multiple 
rank -1 critical points, corresponding, for example, to the merging of the subdomains C 2 , cq, c% with the 
domain b\ (in the subsystem .M 3 , we must also consider the possibility s = 00 ). Obviously, multiple 
fragments also appear in the subsystem A4 1 . A complete analysis of the corresponding diagrams for 
A = 0 is performed in [40]. 


Theorem 18 (P. E. Ryabov). In the classification of bifurcation diagrams T,hk(£,X), the separating 
set 0 l consists of the axes {A = 0} U {£ = 0} and the following five curves determined by the explicit 
single-valued functions: 


7 1 : A = X\(£) 

72 : £ = 4(A) 

73 : £ = 4(A) 

74 : £ = 4(A) 

75 : £ = 4(A) 


|v4 2 + 4 —2s| 

(Vs 2 + 4 — s ) 1/2 
|2A 2 / 3 — V 4 + A 4 / 3 | 
V%(V 4 + A 4 / 3 — A 2 / 3 ) 1/2 5 

42 (ViTa 1 _ a 2 ) 3/2 


1 

2 AV 3 


\J 1 — A 4 / 3 , 


1 

4A’ 


s = (2£ 2 ) 1/3 , £>0] 

A > 0; 

A > 0; 


0 < A < 1; 
A > 0. 


We emphasize again that except for the new curve 75 , all other curves correspond to (3.45): 


7l = 7T23 U 7T 3 i, 72 = 7T 2 4, 73 = 7r 2 2 , 74 = ^21 • 
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The first quadrant is divided into 18 domains in which the diagrams /<(■£. A) considered as one- 
dimensional stratified manifolds are identical. The separating set and the numeration of domains 
appeared are shown in Fig. 25. In what follows, the word “domain” means the connected component 
of the complement of a separating set on the parameter plane, whereas connected components of 
complements to diagrams are referred to as chambers. 

In Figs. 26-31, the diagram T,hk(£, A) for each of the 18 domains of the plane (X,£) is constructed 
in two representations. The first representation contains notation of smooth segments and singular 
points in accordance with the diagrams of critical subsystems. In the second representation, atoms 
of bifurcations that occur in the system with two degrees of freedom on Pf are indicated on smooth 
segments of the diagrams. Arrows on asymmetrical atoms indicate the direction of increasing the 
number of connected components of regular integral manifolds. 

Combining the diagram T,hk(£, A) in the extended space 

A{R 3 {£, h,k)) =M 4 (A ,£,k,h) 
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(see notation (3.19)), we trace the evolution of chambers, i.e., open domains that appear when 
£ hk{ •£, A) divides the plane ( k,h ). For this purpose, we introduce the expanded diagram 

A(£) = U {(A ,i,k,h) : (k,h) € Hhk(£, A)}- 
{\,e) 

The open connected component of the complement of A(£) in M 4 is called an extended chamber. We 
do not distinguish between expanded chambers that differ only by signs of A and £. The integral 
manifolds corresponding to two points of the same expansion chamber are naturally diffeomorphic 
and can be obtained from one another by a smooth isotropy (possibly supplemented by a symmetry 
that changes the signs of A and £). In the sequel, we will not make notes about the signs of A and £. 
It follows that two chambers with given A and £ should be considered the same if they are obtained 
as sections of the same extended chamber. Using the figures, it is easy to see that there are only eight 
different chambers; moreover, one of them (namely, the chamber containing arbitrarily large negative 
values of h) is inadmissible, i.e., the integral manifolds in the preimage of its points are empty. Thus, 
there are seven admissible chambers we are interested in. In the figures, they are marked by Roman 
numerals I—VII. 

Summing up information presented above, we construct 18 bifurcation diagrams Yihk{£, A) of sys¬ 
tems on Pf that are stable with respect to parameters. These diagrams were initially presented in [29]. 
In the construction of these diagrams, the problem of visualization of small domains becomes impor¬ 
tant. Therefore, the diagrams Yihk(£, A) in our figures are severely distorted, and we do not intend to 
preserve the structure of levels h = const, which will allow one to see the rough isoenergetic Fomenko 
invariant in the sequel. Isoenergetic diagrams (see [31, 48]) are more convenient for this; we perform 
a classification of these diagrams below. It is important to see all arcs of diagrams that correspond to 
domains on the planes of the diagrams of critical subsystems and all nodal points. This will be used 
for classification of loop molecules. 

5.2. Topological analysis. The first question of topological analysis is how many Liouville tori 
are contained in the preimage of a regular point in an integral manifold? If we do not know exactly 
the topological type of a critical integral manifold with hyperbolic singularities, we cannot answer 
this question for each specific system on Pf. In a system with two degrees of freedom, the nature of 
a bifurcation that occurs at the intersection of rank -1 critical points is uniquely determined by the 
type of the point only in the case of elliptic singularities, i.e., atoms A. For this problem, only two 
particular cases A = 0 (see [38-40]) and £ = 0 (see [29, 46]) are examined. In the general case (see [14, 
15]), the assertion on the number of tori contained in regular manifolds is based on the numerical 
simulation of their projections onto the plane of the first two components of the angular velocity (an 
analog of Zhukovsky domains). However, considering the entire picture in the four-dimensional space 
R 4 (A,£, k,h) of integral parameter (the extended space), we immediately obtain the answer without 
any additional studies. In fact, we reject all hyperbolic atoms whose exact form is unknown. For 
an extended chamber, an edge of the diagram T,hk(£, A) with a fixed notation generates a (three- 
dimensional) wall. For a selected chamber, an atom A on an edge or on a wall is called an entrance if 
its is directed inside the chamber (i.e., a Liouville torus is generated when entering the chamber) and 
an exit in the opposite case. In the extended space, we consider three sections corresponding to the 
domains 1, 2, and 2' (see Fig. 32). The bottom (with respect to the direction of the h- axis) point of 
the diagram corresponds to one nondegenerate rank -0 critical point of the type “center-center,” which 
is the lowest relative equilibrium for given (\,£). Therefore, as was noted above, there is no motions 
in the unnumbered chamber below the diagram. Consequently, in the chamber I we have exactly one 
torus (each of the walls aq and 6 i contain one entrance). The diagram for the domain 1 shows that 
there exist one entrance from the chamber I to the chambers II and III (through the walls c§ and C 4 , 
respectively). Similarly, from the diagram for domain V we see that there is exactly one entrance from 
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Fig. 26. Diagrams Y,hk {£, A) for the domains 1-3 









































Fig. 27. Diagrams T,hk(^ A) for the domains 4-6 


321 
































322 


Fig. 28. Diagrams T,hk( A) for the domains 7-9 



























Fig. 29. Diagrams T,jjk{£, A) for the domains 1/-3' 
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Fig. 30. Diagrams Y>hk(1, A) for the domains 4'-6' 
























Fig. 31. Diagrams T,jjk{£, A) for the domains 7'-9' 
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the chamber I to the chambers IV and V (through the walls <24 and ci, respectively). Therefore, in 
the chambers II-V, each point corresponds to two Liouville tori. 

The diagram for the domain 2 shows exactly two entrances from the chamber II to the chamber VI 
(the wall as), and the diagram for the domain 2 yields exactly two entrances from the chamber III to 
the chamber VII (the wall as). Therefore, in the chambers VI—VII we have four Liouville tori. The 
following theorem clarifies the results of the papers [14-16], which contain the erroneous assertion 
on the existence of six rather than seven different chambers. In particular, for small nonzero values 
of A, the existence of only five chambers is ascertained in [14-16], although, as one can see from the 
diagrams for the domains 1 and 2 in Fig. 26, for arbitrarily small A, the chamber III of the domain 1 is 
transformed to the chamber VI of the domain 2. This means that in the space R 3 (£, h, k ) for arbitrarily 
small fixed A, the diagram of the total moment mapping divides this space into six rather than five 
connected components with nonempty integral manifolds. If A increases, then for A > 2 -3 / 4 the sixth 
chamber disappears and the seventh chamber appears. 

Theorem 19. Integral manifolds in the preimages of points of chambers of the space of integral con¬ 
stants are as follows: 

(i) T 2 for the chamber I; 

(ii) 2T 2 for the chambers II-V; 

(iii) 4T 2 for the chambers VI-VII. 

In the study of fine topology of the system, one of the main objects, along with chambers and the 
number of tori in them, is the so-called families of Liouville tori. For brevity, we denote the extended 
phase space by P: 

P = A(H 5 ) = [J -P/x{(A ,£)} = J5 2 (q)xM 3 (lj)x{A} = 5 2 (q)xM 3 (lj)xM(A). 

(A/) A 

All integral manifolds A) can be embedded in P as A)x{A}. Naturally, the following 

definition arises. 

Definition 12. We exclude from P all connected components of critical levels containing critical 
points of the moment mapping. The connected component of the remaining set is called the family of 
Liouville tori. 

If the system has no minimum (maximum) regular tori (and there are no such tori in the problem 
considered: this can be easily seen from the classification of diagrams equipped by atoms), then the 
intersection of the preimage of any chamber with any family can consist of no more than one torus. In 
some problems, one can speak of several tori of the same family in the preimage of points of a chamber 
(see [7]) or, what is the same, of families consisting of several tori (see [20]). In this case, the definition 
of a family is as follows: each family consists of all tori that are affected by the same bifurcations on 
the boundaries of the chamber (see [20]). Such families exist in the classical Kovalevskaya case A = 0 
(see [7]) and in the particular case where \I / 0 (see [20]). However, it is easy to see that in the 
general case where A l 0, there are no such families. Indeed, in the chambers II and III, there are 
two tori and a wall with one entrance; therefore, only one torus is affected by a bifurcation. Between 
the chambers IV and V, there is a wall with one entrance and and one exit; therefore, two tori of 
these chambers are affected by clearly different bifurcations. In the chambers VI and VII whose points 
correspond to four tori (not only in the extended space but also for any fixed A and I for which these 
chambers exist), there exist one or two walls with two entrances from neighboring chambers; moreover, 
if there are two walls, then distinct pairs of tori are affected by bifurcations on them. Hence, in this 
case, there are no four or even two tori that are affected by the same bifurcations. In particular, this 
implies that Theorem 6 from [18], which claims the presence of seven families, two of which consist 
of two tori, is invalid. Families in [18] are understood in the sense of the definition given in [20]. We 
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Fig. 32. Number of tori in chambers 


emphasize once again that for any fixed A and l such that \t 7 ^ 0, there are no chambers such that 
all their tori are affected by the same bifurcations on the walls of these chambers. 

Now we find the number of families and their distribution in chambers. We indicate the following 
property of the Kovalevskaya-Yehia case: in the preimage of a point of the bifurcation diagram, regular 
tori cannot coexist with hyperbolic atoms. In other words, if an edge of the diagram corresponds to 
hyperbolic atoms, then all tori of neighboring chambers are involved in bifurcations. Note that this 
property does not take place, for example, in the Goryachev-Chaplygin case. Therefore, in order to 
have the same family in two chambers, it is necessary (but certainly not sufficient) to have a path from 
one of them to another consisting of entrances or exits. If there is only one family in some chamber, 
then the existence of an exit into another chamber is a sufficient condition of the presence of this 
family in this second chamber. 
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Table 5.1 


Chamber 

Family 

I 

{1} 

II 

{1 + 2} 

III 

{1 + 3} 

IV 

{1 + 4} 

V 

{1 + 5} 

VI 

{1 + 2+ 4 +6} 

VII 

{1+ 3 + 5 + 7} 


Let us consider the diagrams for the domains 1, l 7 , 3, and 3 7 ’ (see Fig. 33). Families are marked by 
Arabic numerals and, as is customary for sets, are denoted by braces (in contrast to the numbers of 
domains on the (A, f)-plane). We denote with {1} the family starting in the chamber I. It can be seen 
from the diagram for the domain 1 that there is exactly one entrance into the chambers II and III from 
the chamber I. Therefore, the family {1} also presents in these chambers. Can tori generated at the 
entrances into these chambers belong to the same family? On diagrams with type-A atoms indicated, 
we can see that there is no a path between the chambers II and III that passes through the chamber I 
and contains only type-A atoms (if a path passes through the chamber I, then new families disappear). 
Therefore, the second families in these chambers are distinct; we denote them by {2} and {3}. From 
the diagram for the domain 1' it follows that the family {1} presents also in the chambers IV and V. 

Can tori generated at entrances of these chambers belong to the same family? Obviously, they 
cannot, since the common wall has simultaneously an ebtrance and an exit, so the intersection of this 
wall preserves only the family {1}, and there is no other path between the chambers IV and V that 
does not pass through the chamber I and contains only type-A atoms. Similarly we can establish that 
families generated on a common wall of the chambers IV and V cannot coincide with any of families 
introduces earlier. Therefore, we denote these new families by {4} and {5}. From the diagram for the 
domain 3 we can see that there are two entrances to the chamber VI with four tori from the chambers II 
and IV. Therefore, the families {1}, {2}, and {4} are preserved in the chamber IV and one new family 
appears; we denote it by {6}. Similarly, from the diagram for the domain 3 7 we can see that there 
are two entrances to the chamber VII with four tori from the chambers III and V. Therefore, the 
families {1}, {3}, and {5} are preserved in the chamber VII and one new family appears; we denote 
it by {7}. We join the numbers of families that simultaneously present in the chamber by a plus sign. 

Theorem 20. In the Kovalevskaya-Yehia case, there exist exactly seven families of regular tori. Their 
distribution in the chambers is shown in Table 5.1. 

In the next section, we establish topological invariants of nodal points of bifurcation diagrams. 

6 . Fine Topology of Nodal Points 

6.1. Definitions of rough and fine invariants (molecules). We briefly recall the notion of a 
topological invariant of a three-dimensional manifold equipped on which a Liouville foliation. A de¬ 
tailed description can be found in [7, 8, 35] and in many other subsequent works, where Fomenko- 
Zieschang invariants are described. 

Let Q be a smooth, three-dimensional invariant manifold of an integrable Hamiltonian system sgrad H 
with two degrees of freedom that does not contain rank-0 critical points and degenerate rank-1 critical 
points of the integral mapping T = HxK. This manifold can be obtained by setting a closed path in 
the ( h , fc)-plane that does not pass through nodal points and transversely crosses all one-dimensional 
segments of the diagram (see [23]). Identifying each connected component of the integral manifold 
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or a critical integral surface with a point, we obtain a one-dimensional object, which can be turned 
into a graph whose edges are unions of points corresponding to components without critical points 
and vertices are points obtained from singular sheets (i.e., components that contain critical points of 
the moment mapping). Then we mark vertices by the bifurcation symbols and obtain a graph W(Q) 
with equipped vertices (knots) called the Fomenko graph or a rough molecule. The pair (Q,W(Q))) 
determines the topology of the foliation on Q up to a rough Liouville equivalence. Let us recall the 
definition. 

For a Hamiltonian system v = sgrad H and a three-dimensional invariant manifold, we denote 
by v(Q) the restriction of the field v to Q. Assume that two systems v' and v" are given. The 
restrictions v'(Q') and v"(Q") are said to be roughly Liouville equivalent (see [8]) if there exists a 
homeomorphism of bases of Liouville foliations (i.e., of simple graphs without additional marks), 
which can be lifted to a homeomorphism of Liouville foliations over a neighborhood of each point, ft 
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is clear that in this case a homeomorphism of graphs preserves the marks of bifurcations at vertices. If 
there exists a global homeomorphism of v'(Q') and v"(Q") preserving Liouville foliations, then these 
restrictions are said to be Liouville equivalent. To emphasize the distinction, the Liouville equivalence 
is sometimes called complete, or fine, or exact Liouville equivalence. 

Thus, Fomenko graphs are rough topological invariants (i.e., invariants of the rough Liouville equiv¬ 
alence). It is easy to see that a Fomenko graph can be determined by a geometrical representation 
or as a sequence of bifurcations along a path. Having a sequence of bifurcations along some fairly 
informative set of paths in the ( h , fc)-plane, one can obtain a Fomenko graph for an arbitrary path. 
Such description of the rough phase topology of systems with two degrees of freedom for classical 
problems of rigid-body dynamics was proposed in [38] (a detailed description is presented in [40]). 

As a rule, in problems obtained by reduction by a cyclic coordinate with the moment integral L, 
Fomenko graphs can be constructed for isoenergetic manifolds Qj h of reduced systems. In such 
problems, the term “Fomenko graph” is used for graphs W(Qg h ), i.e., rough isoenergetic topological 
invariants. However, for a more accurate description of the total topology of a system and isoenergetic 
invariants, the method of loop molecules proposed in [7] is very useful. Consider a nodal point of the 
bifurcation diagram of a system sgrad H and surround it by a sufficiently small closed path that does 
not contain other nodal points inside. Naturally, as in the general case, this path must intersect all 
one-dimensional segments of the diagram transversely. Consider the full preimage of this path. Each 
component Qi of the preimage is a manifold satisfying the requirements stated above, and we can 
construct the corresponding graph W(Qf) for it. The resulting set of these graphs for the preimage 
of the chosen path is called the (rough) loop molecule of the nodal point considered. It defines the 
topology of the foliation on Q = U iQi up to the rough Liouville equivalence. 

To obtain a complete (exact) topological invariant, the Fomenko graph must be equipped with 
information about gluing matrices on edges of the graph of tori that came to meet from the bifurcations 
of the atoms at the ends of edges. In [7, 8], a construction of admissible bases of cycles on boundary 
tori of atoms of all types is described. As a result, when meeting on the edge, such bases considered 
as elements of the one-dimensional homology group of the torus are connected by an integer-valued 
matrix with determinant of —1. Since this matrix is determined up to an admissible change of bases 
on the boundaries of the edge, one can introduce its invariants, which are called the r-mark (a rational 
number or oo) and the £-mark (taking values ±1). Thus, each edge of the graph is equipped with a 
pair of numbers (r, e). In addition, the global topology of the Liouville foliation on Q is determined 
by a set of n-marks that are assigned to families of hyperbolic atoms. A family is formally defined 
as follows. Cutting a Fomenko graph equipped with (?’, e)-marks on edges by all the edges with finite 
values of r, we obtain a set of connected graphs. A connected graph without elliptic atoms is called 
a family (atom family, in contrast to tori family). To each atom family, we assign an integer number 
called an n- mark, which is calculated by the set of gluing matrices; in fact, it describes an obstruction 
for spreading of sections of the S^-foliation in the preimage of the family at the boundary tori to the 
whole preimage (by construction, the preimage of each family is a Seifert foliation). The Fomenko 
graph with all three types of marks is called a marked molecule or a Fomenko-Zieschang invariant and 
is denoted by W*(Q). It is shown in [8] that a marked molecule determines the class of fine Liouville 
equivalence of systems on three-dimensional manifolds. For this reason, marked molecules are also 
called fine topological invariants. 

In the sequel, we examine rough and marked loop molecules of nodal points of bifurcation diagrams. 
A considerable amount of information is obtained from the classification of loop molecules of low 
complexity described in [6, 8, 56]. The main problem is the correct tracing of families meeting at 
bifurcation points. Knowledge of the rough topology, i.e., joint levels of first integrals, both regular 
and critical, cannot always uniquely assign marks of families. But, as one can see below, after this 
problem is solved, it becomes possible not only to construct a rough molecule but also to apply results 
of the general theory and studies of particular cases to the construction of all marked loop molecules. 
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6.2. Loop molecules of rank-0 points. We describe the topology of full integral levels of relative 
equilibria (rank-0 critical points). These results is presented in [44] and is based on [47], where 
Tables 4.1, 4.3, and 4.4 were initially obtained. 

Let x be a point of nondegenerate relative equilibrium and t = L{x). Consider the integral map¬ 
ping (2.10) of the reduced system: 

Jt = HxK : Pf -> M 2 . 

We denote the complete preimage of the point J({x) (i.e., a critical integral surface) by Z{x) and a 
sufficiently small saturated neighborhood of the surface Z{x) that does not contain relative equilibria 
with other values of the mapping by U(x). The surface Z(x) can consist of several connected 
components. As was noted above, the component of the point x always contains exactly one relative 
equilibrium. At the same time, an analysis of critical subsystems shows that along with components of 
relative equilibria, components containing exactly one critical circle may also exist. The number and 
topology of these components can be found by the results for the critical subsystem M. \ according to 
Remark 8. Next, we consider the point x in each of two critical subsystems containing it and analyze 
information on neighboring domains in the image of the critical subsystem. These domains produce 
arcs of the bifurcation diagram of the mapping Jt in a neighborhood of the point x. The structure of 
the critical set in the preimages of these arcs and surgeries in U (x) at their intersections are presented 
in Tables 4.1, 4.3, and 4.4. After this, the type of loop molecules of relative equilibria themselves 
and marked critical periodic trajectories lying at the same level can be uniquely established by the 
exhaustive description of loop molecules of nondegenerate singularities of low complexity (see [8]). 
In addition to components containing critical point, the surface Z(x) can contain regular tori filled 
with doubly periodic trajectories. Their number can be uniquely found by the form of the bifurcation 
diagram of the mapping Jt supplemented with atoms on arcs. 

The result of topological classification is given in Table 6.1. For loop molecules, we indicate only 
r-marks with the sole purpose to distinguish molecules of relative equilibria “center-center” (mark 
r = 0) from molecules lying on the same level of elliptic periodic trajectories (mark r = oo). In fact, 
all marks (including e- and n-marks) are set automatically in accordance with [8]. 

On fragments of the bifurcation diagram of the mapping Jt in a neighborhood of a relative equi¬ 
librium, we indicate atoms arising on critical circles at the intersection with arcs of the diagram; only 
atoms of the types A, B, A*, and C 2 appear here. For nonsymmetric atoms, arrows indicate the direc¬ 
tion of increasing the number of tori. This number itself in regular domains is set in the frame. When 
describing singular trajectories (singular components of the integral manifold), pec is a fixed point of 
the type “center-center,” pcs is a critical manifold of a fixed point of the type “center-saddle” (the 
figure-eight), pss is a critical manifold of a fixed point of the type “saddle-saddle” (two figure-eights 
with a common central point and four rectangles glued to them and filled with asymptotic trajectories 
of regular points; the rule of gluing is completely determined by the corresponding loop molecule). 
Further, Sjj denotes a periodic trajectory of elliptic type that exhausts the corresponding connected 
component, and S ] H is a surface of the periodic trajectory of hyperbolic type corresponding to an atom 
of the type B (the direct product of a figure-eight and a circle). This classification by the number of 
classes and the type of loop molecules of the complete preimages of the values of the moment mapping 
at relative equilibria differs from the results presented in [18]. Since the classification principle has 
not been clearly stated in [18], we do not compare the results. We only note that after the publication 
of the work [44], the description of classes of rank-0 points and their loop molecules (still without a 
precise definition of classification principles and necessary analytical studies) was particularly ordered 
in [32], 

According to information from Table 6.1 and a remark about marks on loop molecules, to complete 
a fine classification of all loop molecules of nondegenerate rank-0 points it suffices to mark the edges 
of all molecules found by the numbers of families of regular tori. For molecules of points of the types 
“center-center” and “center-saddle” (with molecules lying on the same integral level of critical circles 
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Table 6.1 (continued) 


Class 
of points 

Number of 
components 
in the preimage 

Singular 
trajectories 
in the preimage 

Fragment of 
the ( H , I\ )-diagram 

Molecules 
in the preimage 

Regular 

tori 

822 , 826 

1 

PSS 

0 _ 

t B 

-► 2B 

0 

f 2B 

B U B U B 

0 B iV 

0 

T 

m ~ 

1 

-► B 

r 

0 

^21 

1 

PSS 

0 J 

t B 

— 2A 

.0 

c 2 

- -A*-J 

B i/2 ^C 2 

0 

T 

m ~ 

-►B 

0 


that consists of nondegenerate rank-1 points), this problem is solved directly based on the form of the 
adjacent fragment of the equipped bifurcation diagram T,hk(£, A). 

For example, consider the point from the class 623 , which can be found on the diagram for the 
domain 6. The atoms on the adjacent segments are as follows: B ^ as ], 2-B[ a7 ] and A[ Cg ], 2A[ Cg ]. Here 
and in the sequel, we use the notation of atoms supplemented by the notation of the edge of the 
diagram (in brackets) at which the bifurcation considered occurs. Since the rank-0 point itself has 
the type “center-saddle,” its loop molecule (i.e., a connected molecule corresponding to the connected 
component that contains the rank-0 critical point) describes the bifurcation 

A [c 6 ] B la 7 ] 2 Aesb 

moreover, the family {2} appears in the atom A[ C6 j and the families {2} and {6} disappear in the 
atoms 2y4[ cg j. Thus, in the atom By 7 j that belongs to the connected molecule of the point <^ 23 , the 
family {2} is transformed into the families {2} and {6}. The smooth segment B^ — B j a7 j overlays 
the corresponding diagram; this segment corresponds to a nondegenerate critical saddle circle. As is 
seen from the fragment of the diagram, its left atom B [ a3 j separates the chambers I and IV, so the 
family {1} is transformed into the families {1} and {4} through a critical circle in the atom U[ a3 ], and 
hence in the atom B[ ay y which belongs to the second component of the preimage of the loop molecule 
(i.e., the component that does not contain the rank-0 point). 

In the sequel (for the molecule of the points 622 an d A 13 ), we will need information on the surgery 
of families in the atoms B^ a7 y For this reason, we state this as a separate proposition. 

Proposition 21. In the two atoms B on the edge 07 , the following bifurcations hold: on one atom, 
a torus of the family {1} is transformed into two tori of the families {1} and {4}; on the other atom, 
a torus of the family { 2 } is transformed into two tori of the families { 2 } and { 6 }. 

Information on loop molecules in preimages of nondegenerate points of the types “center-center” and 
“center-saddle” is collected in Table 6.2. For convenience, the first column, along with the notation 
of the class of points, contains the number of a domain on the plane (A,t?) whose diagrams contain 
points of this class (certainly, such a point can be contained on the diagrams for other domains; for 
example, points of the class exist on all diagrams). 
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Table 6.2 


Class 

of points 


Class 

of points 


Molecule 


Molecule 


Si 

( 1 ) 


p]p A[fl > ] 


^25 

(10 


A [c,] pp- A[a 9 ] 


A [a 4 ] pp“ A[a 9 ] 


x" 

°28 


(60 


A K]pp- A[cj] 


A [a 12 ]pp- Ap] 


^32 

( 6 ) 


A[6 2 ] pp A[a 4 ] 


A[Z> 2 ] pp- Ap] 


<^24 

(3) 


A[fl 8 ] pp- A[ C J 


A [c 6 ] p^ A [c j 


A[aJ pp- A [a 4 ] 




28 


(60 


A [c 3 ] pp- A[ fll0 ] 


A[aJ pp- A[ a , 0 ] 


pp- A[cJ 


£31 

(9) 


{4} 


.A[6 2 ] 


{5} 


A[a 12 ] 


{1} 


A[6,] - B [a 3 ] 


,{i} 


S 27 

(60 


{ 1 } 


A [a,] - Bp] 


Tl} 


A [b 2 ] 


A[a 12 ] 


<^23 

( 6 ) 



A[a 10 ] 


S'o 


27 


(60 



B[fl 3 ] B [«0 

{i} 


B[ C J B[c 2 ] 

m 


334 

































Consider loop molecules of the remaining three classes of nondegenerate rank-0 critical points of 
the type “saddle-saddle.” 

Let us take the point 622 ■ It is contained, for example, in the diagrams for the domains 2 and 3. 
We already know that the surgeries of families on the adjacent edges of the bifurcation diagram are 
as follows: 

b [C7 ] '■ { 1 } | -t { 1 } + { 2 }; 

B[a 3 ] '■ {1} | -t {1} + {4}; 

2-B[a 7 ] : {1} {1} + {4}, {2} !-)• {2} + {6}. 

The remaining edge C 9 separates the chambers IV (with the families {1} + {4}) and VI (with the 
families {1} + {2} + {4} + { 6 }). Unfortunately, we cannot establish the rule of marking the four 
edges connecting 2 L>[ a7 i and 2 B^ cg j by families without considering degenerate rank -1 critical points. 
However, according to Theorem 17, we know that the loop molecule of the point A 32 has two connected 
components; therefore, each of these components has a well-known form, namely, the family in the 
“head” again comes to one of the “feet.” But the point A 32 is a cusp at which the hyperbolic edge C 9 
and the elliptic edge cs come in the diagram of the domain 7. Hence we immediately see that in the 
atoms -B[ Cg ], the surgeries of families are as follows: 

{1} 1 y {1} + {4} V { 6 }, {4} ^ {2} + {4} V { 6 } 

(the notation {4} V { 6 } means an uncertain alternative). Thus, the following logical problem appears. 
There are two males Mi and M 2 (the families in the “heads” of the atoms 2-£>[ a7 ]) and two females F\ 
and F 4 (the families in the “heads” of the atoms 2L>[ cg ]). It is known that each male has exactly one 
common child with every female; we denote these children by D 1 , D- 2 , D 4 , and Dq (the families on 
the edges of the loop molecule of the point S 22 between <27 and ag). It is also known that D\ and D 4 
are children of M 1 and D 2 and Dq are children of M 2 . Further, we also know that F± is the mother 
of D± and F± is the mother of D 2 . We must identify the pair of parents (namely, the mother) of each 
child. The answer is now obvious: the parents of D 2 are (M 1 , ip), the parents of D 2 are (M 2 , F 4 ), the 
parents of D 4 are (M \, F 2 ), and the parents of D$ are The bifurcation in the atoms -B[ C9 ] is 

important for the sequel, so we describe it in the following assertion. 

Proposition 22. The surgeries of families in the two atoms H[ Cg ] are as follows: 

{l}^{l} + {6}, {4}^{2} + {4}. 

Finally, we obtain the first molecule in Table 6.3 for the point 622 - 

Consider the point 626 - Here we do not know the rule of marking by families of the four edges that 
connect 2B\ ail \ and 2L>[ C 5 i. The edge an separates the chambers V and VII, i.e., the families {l} + {5} 
are transformed into the families {1} + {3} + {5} + {7}. Note that if the values of the parameters 
belong to the domain 4', then the edge an has an exit to the point An, at which, according to 
Theorem 16, the loop molecule has two connected components, so that the families in the “heads” of 
the atoms L>[ an ] come to the “feet” of the same atoms. In particular, the bifurcation of the families 
in the atoms B[ an ] are as follows: 

{1} 1 y {1} + {3} V {7}, {5} 1 y {5} + {3} V {7}. 

On the other hand, from the loop molecule of the point 8 ' 2 7 (see Table 6.2) we know that the bifurcations 
in the atoms B[ C5 i are as follows: 

{lW{l} + {5}, {3}^{3} + {7}. 

Since, according to the atoms set on the edges of the bifurcation diagram, the number of regular tori 
in the adjacent chambers, and the general classification (see [ 8 ]), the neighborhood of the point 825 
has the type BxB, the alternative in the atoms is solved unambiguously. 
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Table 6.3 


Class 
of points 

Molecule 

Class 
of points 

Molecule 

^22 

(2) 

B[ c ] B[c 9 ] B[c 9 ] 

B [a 3 ] 

<^26 

(2') 

l ^ b ki^) 

B[c ] B[. ; ] B[ c ,] 

I L 2J Bee r = 0 / / 

{ ^- b K] 

821 

(1) 

f w Vi 

1 /Xn ' 

B [c 7 ] C 2[ c J 

\r = 0 

{1 * B [aj 




Proposition 23. In the two atoms B^ ail j, the surgeries of families are as follows: 

{l}^{l} + {3}, {5}^{5} + {7}. 

As a result, for the point 8 26 , we obtain the second molecule in Table 6.3. 

For the last nondegenerate point 621 (the diagram for the domain 1), the question is how one can 
join families in the bifurcation 2A* a ^. The edge 05 separates the chambers II and IV, i.e., the pair of 
the family {1} + {2} is transformed into the pair {1} + {4}. The rule of the correspondence of families 
will be uniquely deduced below in the study of the loop molecules of the point A 01 , which is always 
one of the endpoints of the edge 05 (the diagrams for the domains 1-6). We state it now without proof 
(for the proof, see p. 337). 

Proposition 24. In the two atoms A* ag p the surgeries of families are as follows: 

{ 2 } ^{ 4 }. 

Again, we use the atoms on the edges of the bifurcation diagram, the number of regular tori in the 
adjacent chambers, and the general classification (see [ 8 ]) to determine the third molecule for a point 
821 in Table 6.3. 

This concludes the study of loop molecules of nondegenerate rank-0 critical points and the structure 
of their full preimage under the moment mapping. We emphasize that Proposition 24 has not been 
proved. 

6.3. Loop molecules of rank 1 degenerate points. The list of marked loop molecules of rank- 
1 singularities is given in [32]. Although the classification in [32] is purely visual (by the shape of 
bifurcation diagrams) without an analytical description of separating values of the parameter, this list 
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8=1 


Fig. 35. Loop molecule for Aqi 


is exhaustive. The marks of molecules are mainly derived from already existing marks by the addition 
rules for marks (see [8]), but in some cases they are calculated by the Topalov formulas (see [34]), 
which requires knowledge of the topology of the loop manifold of a degenerate rank-1 critical trajectory 
(in particular, the number of connected components, which was not previously found proved in all 
cases). In this section, we present the classification of loop molecules with full argumentation. 

6.3.1. Points Aq. We construct loop molecules for degenerate closed trajectories in the preimage of 
the set Ao- 

Consider the point Aoi, taking, for example, the diagram for the domain 1. The atoms on the 
edges of the diagram in a neighborhood of this point are indicated in Fig. 26 and families of tori in 
chambers in Fig. 33. All this information is collected in Fig. 34(a). 

The graph of a molecule with vertices-atoms and edges-families can be restored by this information 
and gluing rules for families of two atoms A*\ this is shown in Fig. 34(b). However, we see that 
there are only two options of gluing, and if we glue the point E with F' and F with E', then we 
obtain a connected molecule, which is impossible, since, according to Theorem 15, the surface Joi 
consists of two components. Therefore, we must glue the point E with E' and F with F'. Then we 
obtain a molecule consisting of two components, as is shown in Fig. 34 (c). This, in particular, proves 
Proposition 24.. 

Since the assertion from [20] about the rough structure of the molecule has been rigorously proven, 
we can also use marks calculated in [20] (the point Aoi corresponds to the point Z 3 in [20]). The final 
result for the point Z 3 is shown in Fig. 35 There is a difference with [20] in the gluing of families. Thus, 
in the second component of the atom A*, tori of the same family converge. This cannot be established 
within the manifold £ = 0, which was studied in [20]. In the present paper, families are treated in the 
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Fig. 36. Construction of a loop molecule for A 02 


A 


[« 2 ] 


{ 1 } 





Fig. 37. Loop molecule for A 02 


sense of the extended phase space. Therefore, two distinct families of tori in a one-parameter set of 
systems with two degrees of freedom (. = 0, A € M turn out to be the same family in the two-parameter 
set of such systems. 

Consider the point A 02 on the diagram for the domain 1 . As in the previous case, the atoms on the 
edges of the diagram in a neighborhood of this point are indicated in Fig. 26 and families of tori in the 
chambers in Fig. 33. All this information is collected in Fig. 36(a). Here rough molecules and families 
on the edges can be restored uniquely (see Fig. 36(b)). Marks on edges can be found in [7, 8 ], where 
only the r-marks are presented; however, all e-marks are equal to + 1 , according to the corresponding 
result of [19]). The resulting molecule is shown in Fig. 37. The discrepancy between the numbers of 
families with [7] is due to the same reason as in the previous case: the impossibility of identifying 
some of the families in the framework of a particular case (here the particular case is A = 0). We will 
make no further remarks about the number of families when comparing molecules. 

Consider the point A 03 . It can be found in the diagram for the domain 2. The atoms on the edges 
of the diagram in a neighborhood of this point are shown in Fig. 26, the chambers and the number of 
tori in them in Fig. 32, and the families are shown in Table 5.1. All this information is collected in 
Fig. 38(a). Gluing of families in the atoms A* of the edge a 5 is described by Proposition 24, which was 
proved above (see the construction of the molecule for the point A 01 ). The families {4} and { 6 } arising 
on the atoms A of the edge as come to the atom B of the edge eg in accordance with Proposition 22. 
The choice of alternative gluing of pairs of the points E , F and E', F' is determined uniquely: the 
point E cannot be connected with F' since a connected molecule appears, which is impossible due to 
Theorem 15: it must consists of two components. Thus, at this stage we have a rough molecule shown 
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Fig. 39. Loop molecule for A 03 


in Fig. 39(c). The marked molecule is shown in Fig. 39. Since the point A 03 exists also for A = 0, 
the marks on the edges of the molecule can be taken from [7], taking into account the results of [19] 
regarding the e-marks. 

Consider the last of the points on Ao, the point A 04 ; it can be found on the diagram for the 
domain 2'. The atoms on the edges of the diagram in a neighborhood of this point are indicated 
in Fig. 29, the chambers and the number of tori in them in Fig. 32, and families corresponding to 
chambers are presented in Table 5.1. All this information is collected in Fig. 40(a). Here the only 
possible ambiguity is related to the surgery of families of the two atoms B of the edge an. Above, 
we have proved Proposition 23 on bifurcations of families on the edge an, using a forecast for a fairly 
simple molecule of the point A 14 . We prove it once again for other reasons. Assuming that at least 
one of the families {5} and {7} generated on the edge C 3 comes to the “foot” of the atom £?[ ail ], which 
transforms two tori into one torus of the family { 1 }, then the molecule obtained is connected, which 
contradicts Theorem 15. Then none of the points G' and H' can be glued with the points E or F 
(Fig. 40(b)). Therefore, the unique possible result is shown in Fig. 40 (c). The molecule for this point 
also exists for £ = 0 ; therefore, marks on its edges can be taken from [ 20 ] (the corresponding point 
was denoted by 24 in [20]). As a result, for the point A 04 we obtain a loop molecule shown in Fig. 41. 

6.3.2. Points A±. The following class of degenerate rank-1 points consists of point lying in the 
preimage of the cusp of the surface Hi. Marks of the loop molecule of this type for cusps on the 
diagrams were described in [20], namely, r-marks are equal to 00 on the edge B — B and are finite on 
the edge A — B\ all e-marks are equal to +1. Therefore, for the point An (see, e.g., the diagram for 
the domain 1), we have the result shown in Fig. 42. Here the r -mark on the edge A — B is equal to 
zero (see [ 8 ]). The vanishing of r-marks on the edge A[ a4 ] — R[ a3 ] for the molecule of the point A 12 
shown in Fig. 43 and lying on the diagram for the domain 3 must be proved. We apply the method 
used in [20]. Recall the addition rule for marks (see [ 8 ]). 
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Fig. 40. Construction of a loop molecule for Aq 4 



Fig. 41. Loop molecule for A 04 


Lemma 4. Let some family be involved in bifurcations of three walls u\, U 2 , and U 3 of the same 
chamber and for the transition u\ — >• U 2 we have r = 00. Then the r-marks of the transitions u\ —> U 3 
and U 2 —>• U 3 coincide. 

Proposition 25. In the loop molecule of the point A 12 on the edge joining the atom A with the saddle, 
the r-mark is equal to zero. 

Proof. We use Lemma 4. For the point A 12 , we consider the transition in the family {4} from -B[ a3 ] 
to A\ a A. Using the bifurcation diagram, we see that these two walls can coexist with various other 
walls. For example, in the domain 3, we have the transitions 

B [a 3 ] 2 Acs] Am]' 

However, the family {4} is not involved in the bifurcation on the edge cs- The point A 12 is also present 
in the domain 1 /; however, we cannot indicate a decomposition of the desired transition into the sum 
of two transitions. 

Consider the diagram for the domain 9. The point A 12 itself does not lie on this diagram, but there 
are transitions in the family {4} of the following form: the transition _B[ a3 ] —>• 4l[fc 2 ] belongs to the 
molecule of the nondegenerate point <$31 with mark r = 00, and the transition A^ —>• B^ belongs to 
the molecule of the nondegenerate point <$32 with mark r = 0. Therefore, the missing r-mark is equal 
to 0 . □ 

Consider the point A 13 . At this point, the following edges of the diagram converge (the domain 2): 
the edge as with two atoms A, in which tori of the families {4} and { 6 } appear, and the edge 07 with 
two atoms B, in which each of tori of the families {1} and {2} is transformed into a pair of tori selected 
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(a) (b) 

Fig. 42. Diagram and loop molecule for An 




(a) (b) 

Fig. 43. Diagram and loop molecule for A12 


from the families {1}, {2}, {4}, and {6}. The endpoint of this edge is the nondegenerate point $23; 
in the study of this point, we have proved Proposition 21 on surgeries of families in two atoms B^ a7 1. 
As a result, we obtain a loop molecule consisting of two components for A13; it is shown along with 
a fragment of the diagram in Fig. 44. 

Note that the fact that the loop molecule of a degenerate circular orbit consists of two components, 
even in the case where the diagram contains a return point, must be proved since there is no complete 
classification of such molecules. Earlier, we have proved this with the help of an analytical solutions 
in Theorem 16. But now it turns out that a molecule consisting of one component could arise only 
in the case where two tori of families that appear on the edge as meet in the same atom B of the 
edge 07, which contradicts Proposition 21. Thus, the presence of two components in J13 is also proved 
topologically. Therefore, marks can be taken from [8]. 

For the point A14, the same question about the number of connected components in a loop molecule 
appears. The presence of two components was proved earlier with the help of analytical solutions in 
Theorem 16. On the other hand, in this case two atoms B are realized on the edge an (see the diagram 
for the domain 4'). The endpoint of this edge is the point A04 examined above; thus, we know that 
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(a) (b) 

Fig. 44. Diagram and loop molecule for A 13 



(a) (b) 

Fig. 45. Diagram and loop molecule for A 14 


on one atom B a torus from the family {1} is transformed into two tori of the families {1} and {3}, 
whereas on the other atom a torus from the family {5} is transformed into two tori of the families {5} 
and {7}, which is indicated in the components of the loop molecule in Fig. 41. For the point A 14 we 
obtain the fragment of the diagram and two components of the loop molecule shown in Fig. 45; we 
obtained this result by using topological arguments. Marks are taken from [20]. 

6.3.3. Points A 3 . The last class of rank-1 degenerate points consists of points lying in the preimage 
of the cusp edge of the surface II 3 . 

As for the set Ai, if the preimage of the point A 31 contains only one degenerate circle, i.e., the 
corresponding critical integral surface and the loop molecule are connected, it suffices to find the 
numbers of families that are involved in bifurcations. Therefore, for the points A 31 and A 33 (see, e.g., 
the diagrams for the domain 1 and F), we obtain the result shown in Figs. 46 and 47. Marks of the 
loop molecule of the point A 33 are known from [7, 20]. For the point A 31 , the infinite r-rnark is a 
consequence of the general statement from [ 20 ] cited above; the fact that the r-mark equals zero on 
the edge A^ C1 ] — _B[ C2 ] requires further argumentation. 

Proposition 26. In the loop molecule of the point A 31 on the edge connecting the atom A with the 
saddle, the r-mark is equal to zero. 
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(a) 


(b) 


Fig. 46. Diagram and loop molecule for A 31 




(a) (b) 

Fig. 47. Diagram and loop molecule for A 33 


Proof. Again, we apply the addition rule for marks in the case where one of the r-marks is infinite 
to the point A 31 and transitions in the family {5} from -B[ C2 ] to Ay Cl y We can also consider various 
diagrams. For example, in the domain 5, there are transitions A[ Cl ] —>• A[ ag ] —> B^ C2 y However, both 
marks are finite: for the first transition, from the molecule of the point S 25 we have r = 0 and for the 
second transition, from the molecule of the point A 03 we see that r = 1/2. Therefore, the “sum” of 
marks is defined nonuniquely. The point A 31 also lies in the domain T, but, as for the point A 12 , we 
cannot indicate a required decomposition of the transformation in the sum of two transformations. 

Consider the diagram for the domain 9h The point A 31 itself does not belong to this diagram, but 
there are transitions in the family {5} of the following form: the transition H[ C2 j —>• A[ Q12 ] belongs to 
the molecule of the nondegenerate point 5'f 7 with mark r = 00; the transition Ar ai2 ] —>• A ci i belongs to 
the molecule of the nondegenerate point 6" 8 with mark r = 0. Therefore, the missing r -mark is 0. □ 

Consider the point A 32 . We take a fragment of the diagram with marked edges and the correspond¬ 
ing atoms from the diagrams for the domain 4. For this point, we must find the number of connected 
components of the molecule (there is no analytical solution for the full integral surface). However, 
a molecule is connected only if the families { 2 } and { 6 } generated on the edge cs come to the same 
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6=1 6=1 


(a) (b) 

Fig. 48. Diagram and loop molecule for A 32 



6=1 6=1 

(a) (b) 

Fig. 49. Diagram and loop molecule for A 34 


atom B of the edge eg. This contradicts Proposition 22. Thus, the molecule of the point A 32 consists 
of two components, and Proposition 22 determines the gluing rules for families. Marks are calculated 
in [ 8 ]. The result is shown in Fig. 48. 

For the point A 34 (the diagram Tihk(^i A) for the domain 2'), the question on the number of 
components and the gluing rules for families is solved on the edge C 5 containing two atoms B. The 
endpoint of this is edge is the point 825 , which is the image of a nondegenerate rank -0 critical point. 
At the same time, if we consider the diagram T,hk(^, A) for the domain 3', we can see that the edge C 5 
also has an exit to the point which is the image of another nondegenerate rank -0 critical point. 
Thus, as was shown above, under the bifurcation of families in the two atoms H[ Cs 1 , a torus of the 
family {1} is transformed into two tori in the families {1} and {5}, whereas a torus of the family {3} 
is transformed into two tori of the families {3} and {7}. Since tori of the families {5} and {7} are 
generated on the edge C3 and these tori come to distinct atoms B of the edge C5, the molecule of the 
point A 34 has two connected components. The gluing rules for families are stated uniquely. Having 
the proof of the existence of two components, we can take the marks from [20]. The result is shown 
in Fig. 49. 
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7. Isoenergetic Bifurcation Diagrams 

7.1. Separating sets and isoenergetic diagrams. In this section, we classify bifurcation dia¬ 
grams Z LK {h, A) of the mappings LxK restricted to isoenergetic levels Q\ = P _ 1 (/i) in the complete 
phase space P 5 (see [48]]). The topology of these levels can be obtained trivially. Since H on P 5 is a 
Morse function having exactly two critical points 


u> = 0 , a = (± 1 , 0 , 0 ), 


we have 


Qh 


' 0 , h<- 1 , 

< S' 4 , -l<h<l, 

S 2 xS 2 , h> 1. 


We obtain the separating set @h, or the P-atlas, immediately from Proposition 12, taking into 
account the results of all facts obtained above about critical subsystems: this set is the image on the 
(A, h)-plane of the above-mentioned singular points (4.11), (4.16)-(4.21), (4.23)-(4.27), and (4.33) of 
the key sets of critical subsystems. Recall that the point P 4 is introduced in the classification of the 
diagrams Hhk(£, A) and does not cause surgeries of the diagrams on Q^ , but an additional point P 5 
appears, which is substantial. As above, it suffices to consider A > 0 to provide illustrations in the 
plane (h, A 2 ). 


Theorem 21 (M. P. Kharlamov, I. I. Kharlamova, and E. G. Shvedov). The separating setQn in the 
classification of bifurcation diagrams T*LK{h, A) consists of the following 13 curves in the plane (A ,h): 


£0 : h 

Ci : h 

C 2 : h 
C 3 : h 
C 4 : A 
C5 : h 
?6 : A 
?7 : A 
Cs : A : 
C9 ■ h- 
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In this formulation, separating curves are listed in the sequence of the generating points A, B 1 -B 7 , 
C , -D1--D3, and -D5. It is easy to see that £4 and are combined into one smooth curve and £,q and £7 
are combined into a curve with a cusp. 

Separating curves with their numbers are shown in Fig. 50. In general, they divide the half¬ 
plane A > 0 into 34 domains, one of which, namely, the domains where h < —1, corresponds to 
empty isoenergetic manifolds Qf L . The other domains are marked with the numbers 1-33. Due to 
the sufficiently complicated structure, fragments of the admissible set with complete enumeration of 
domains are presented in Fig. 51. For convenience, this is also done on the (h, A 2 )-plane. 

Figure 52 shows the equipped isoenergetic diagrams. They are symmetric with respect to the axis 
l = 0 , so it is possible to mark all smooth segments simultaneously (for which bifurcations and gluing 
rules for families have been previously established) on the right-hand side and indicate the chambers 
and bifurcations on the left-hand side explicitly. This yields a complete presentation of the isoenergetic 
two-dimensional invariant of the restriction of the system to the appropriate set Qj L ', moreover, this 
is convenient for the construction of all isoenergetic one-dimensional invariants (Fomenko graphs) of 
reduces systems on P 4 . 

The problem of computer visualization of “ultra small” domains in the construction of bifurcation 
diagrams is well known. Keeping in mind the the use of the diagrams Tilk(}i, A) for the construction 
of Fomenko graphs (to do this, one must correctly see all sections of the admissible domain by vertical 
lines t = const that differ by their structure), we try to keep these sections in the figures where 
possible, i.e., to use distortions of the diagrams by “hberwise” diffeomorphisms. Obviously, this is 
not always possible without loss of clarity. In such cases, the requirement of clarity was of primary 
importance. Major deviations from hberwise (i.e., preserving the sections l = const) diffeomorphisms 
are observed in the diagrams for the domains 27-33. For example, the edge C 4 actually lies on the 
curve, which is a continuation of the edge eg and, therefore, it should be “shorter” than 05 . As a 
result, the absolute value of the ^-coordinate of the point A 02 is smaller than the ^-coordinate of the 
point A 03 , but this is not so in figures. Moreover, the relation between the ^-coordinates for the pairs 
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Fig. 51. Fragments of the set Qh and generated domains 


of the “tails” (A 31 , An) and (A 31 , A 14 ) can be violated, and so on. At the same time, of course, the 
equipped isoenergetic diagram itself remains valid as a rough topological invariant. 

For an accurate construction and analysis of all sections of the diagram £ by fixed h and £, software 
packages in the Wolfram Mathematica system were created (see [31]); they provide an opportunity to 
position the point (h, A) in the interactive session with sufficient accuracy with respect to the separating 
set Qh, selection of the required value of £, observing all relevant phenomena in the key sets of critical 
subsystems, and output the Morse-Bott indexes that define the bifurcation by information from the 
tables above for all singular points on the section t = const of the diagram £ Lx{h, A) in accordance with 
the formulas obtained above. A similar software package for the diagrams T,hi<(£, A) was described 
in [50]. Of course, results of such computer simulations can only be considered as auxiliary; they 
provide an illustrative “constructor of topological invariants.” All results on the classification of 
Fomenko graphs should be strictly justified analytically. This is done in the following section. 
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Fig. 52. Equipped diagrams T*LK(h, A) 
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Fig. 52 (continued). Equipped diagrams A) 
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Fig. 52 (continued). Equipped diagrams A) 
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Fig. 52 (continued). Equipped diagrams A) 
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Fig. 52 (continued). Equipped diagrams T,LK(h,\) 
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7.2. Classical Kovalevskaya case. Let us consider the limit case A = 0 (the classic Kovalevskaya 
case). In Fig. 50, we see the separating values of h: 

- 1 , 0 , 1 , y/2, y/3, 2 . 

Initially, the (L, I\ )-diagrams for the Kovalevskaya case were considered by P. Richter and his disciples 
who used powerful graphic processors for visualization of diagrams, bifurcations, and Fomenko graphs 
in the 1990s. As a result, a movie was filmed in which the typical motion of bodies for various 
families of two-dimensional tori of the reduced system was shown. The mathematical description of 
backgrounds of this movie based on the results of [7, 39, 40] is presented in [67]. In particular, in [67], 
detailed isoenergetic diagrams were presented that were calculated on the (£, \/fe)-plane. Separating 
values of h in the work were not explicitly listed, but they can be derived from the formulas provided 
there. These values were also obtained as limit values in the research on isoenergetic diagrams of the 
Kovalevskaya top in a dual field (see [49]). The diagrams X_Lic(h, A) corresponding to nonseparating 
values of h (the exit to the axis A = 0 from the domains 1, 3, 9, 26, 27, 30, and 31) are shown in 
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Fig. 53. Equipped diagrams 0) (Kovalevskaya case) 


Fig. 53. Here, due to the vertical distortion, the fiber-wise (i.e., preserving the sections t = const) 
diffeomorphism of the domain of the diagram is saved in sufficiently complex cases 27o, 30o, and 31o- 
We also note the following interesting feature. In [67], for the domain 30o, two diagrams were 
presented. Obviously, the corresponding system on Q\ l are Liouville equivalent, but they have different 
sets of Fomenko graphs since the projection of the point 622 on the axis Oh can lie either on the edge 
eg + b\ or on the edge cs + 62 . 


8 . Topological Invariants 

8.1. Diagrams of Smale—Fomenko. In order to classify Fomenko graphs on Q\ h { A), we must 
superpose the projection of the image of the set of degenerate rank -1 critical points on the plane of 
the Smale diagram. Indeed, the transition through such points determines surgeries in graphs that are 
not related with change of the topology of Q\ h {\). Theoretically, splitting of certain atoms without 
appearing degenerate points is also possible (see [ 8 ]), but, as was shown above by listing all atoms 
appearing in this problem, this phenomenon does not occur here. We write the equations of the image 
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Fig. 53 (continued). Equipped diagrams 0) (Kovalevskaya case) 


of the set of degenerate rank -1 critical points taking into account conditions of the existence of motions 
of required types in accordance with Propositions 15 and 20. 

Recall the notation (3.51) and information on the images of rank-1 degenerate points. The col¬ 
lection of curves that classify Fomenko graphs contains the following curves (see Theorem 13 and 
Proposition 20): 

(1) the whole curve of touching of the surfaces IPi and II 3 : 

U = ±^(l-2\*s), 

° ' 1 - A 2 s + 2s 2 

h = -, 

2 s 

( 2 ) the cusp of the surface IR between its intersections with the curve S 3 (i.e., between the points 
that are the images of the class 7 r 3 i of rank -0 degenerate points turning into the points £>3 on the 


0 < s < 


2A 2 ’ 


( 8 . 1 ) 
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key sets): 

A i: <-= ± 5 - 7 s (ft- T ) , T <h<hc(X), 

where the function hc( A) is determined by (4.6); 

(3) the cusp of the surface II 3 within the following limits: 


where 


A3 : 


t 



h G 


[h*,+oo) 
[/*■**, +00) 


for A < A*, 
for A > A*, 




( 8 . 2 ) 


(8.3) 


Definition 13. The Smale-Fomenko diagram S' LH is the union of Slh with the image of the set of 
rank-1 degenerate points under the mapping LxH. 


Such diagrams were first considered by A. A. Oshemkov [23] for classical problems of rigid-body 
dynamics. 

Since the set of intersection points of the curves A, with the principal diagram consists of rank-0 
degenerate points, its transformations with respect to A are already taken into account. To classify 
extended diagrams, we need to consider the values of A corresponding to surgeries of Ao U Aj U A 3 . 
Consider the evolution of self-intersections and double crossings of these curves (there are no triple 
intersection, Ao n Ai n A 3 = 0). The curves Ai and A 3 are uniquely determined by the single¬ 
valued functions h(£) and hence do not have self-intersections. The curve Ao has a self-intersection 
corresponding to the values 


s± = 


1 ± x/1 - 8 A 4 


h = 


1 — A 4 




(8.4) 


4A 2 ’ 2A 2 

where s± belong to the interval (0; 1 / 2 A 2 ) (if they exist, i.e., if A < A* = 1 / 2 3 / 4 ). 

We represent the condition of intersection Ao fl A 3 in the form 

2 s(h - h* - 2£ 2 ) = 0 , 

where h and i are defined by (8.1). Performing the substitution Z = A 2 / 3 s, we arrive at the equation 
(Z + 1)(2Z — l ) 2 = 0 whose unique positive root Z = 1/2 is double; so we obtain a tangent point with 
the coordinates 

1 r.—o/q „,o/q ,oi „ . 1 n o/o ,r,/o /ri 

(8.5) 


h = 


2 1 


A " 2 / 3 + 2A 2/3 - A 2 


= ±-a/ A -2 / 3 — A 2 / 3 ; 
2 


this point exists if A < 1. Since 1 < A*, we must verify that h > h* at the point found. For A < 1 and 
h determined by (8.5) we have 

1 - A 4 / 3 


h - h* = 


2A 2 /3 


> 0 . 


It is easy to see that the points (8.5) are the images of rank-0 degenerate points corresponding to the 
curve 7 T 2 i (see Proposition 3). 

We represent the condition of intersection Ao FI Ai as follows: 

12 \ 3 


2 s d 


4 h- 


y 


- 27f 


= 0 , 
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( 8 . 6 ) 


where h and i are determined by (8.1). We have 


(8s 2 - 2 A 2 s + l) (s 2 - 2 A 2 s + l) 2 = 0. 


The roots of the first factor are 


then we have 


s = \ [A 2 ± \/A 4 - 8] ; 

O L J 


a 2 = - 


T[v±yv3i] ;J <o. 


i.e., these solutions are extraneous. The second factor in ( 8 . 6 ) yields a tangent point and provides a 
positive root 

s = \J A 4 T 1 — A 2 , 

for which the points 

h = lv A 4 + l-A 2 , e = ±Jl(V a^+T-a 2 V 3 


exist for all A; they are the images of rank-0 degenerate points corresponding to the curve 1 x 22 (see 
Proposition 5). 

Consider the intersection Ai n A 3 . The system (8.2), (8.3) without restrictions has three solutions 




The pair of points (8.7) always exists. Obviously, at these points 

h - h* = A 2 > 0, h- h** = i (4 + A 4/3 ) 3 A 2/s - \[i + A 4 / 3 . 

The last difference must be positive for A > A*, but it is negative only if A < 1/2 3 / 4 = A* < A*. 
Therefore, these point satisfy all constraints. Note that, in accordance with (8.4), these points he on 
the the same levels of £ as self-intersection points of the curve Ao (if the latter exist). 

At the points (8.9), under the realness condition A < 2 4 / 3 , the value of l 2 is negative, so this solution 
is extraneous. 

Consider the solution ( 8 . 8 ) for all real A < 2 4 / 3 . Let us verify the conditions for h on the curve Ai. 
Excluding x from (4.6) and setting h = he in ( 8 . 8 ), we obtain the equation 

1 + 48A 4/3 + 12A 8/3 - 8 A 4 - 6^3(2 - A 4 / 3 )(A 2/3 + 4A 2 ) = 0. (8.10) 

By the substitution X = A 4 / 3 we reduce it to the equation 

(■2X-lf(2X-1-3.2'/*+3-2^) (iX + X - JL - ^ ’ + |L (i + -L + 2^ =0. 
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We see that (8.10) has exactly two positive roots: a triple root A = A* and a simple root 

/ 45 99 161\ 1//4 

A 3 =(^-^-—) ~ 0-0287, A, < A 3 < A*. (8.11) 

The condition h < he holds for A > A 3 . The root A = A* is irrelevant to this inequality; its 
appearance is related to a singular point (3.37) defined by the same equation (4.6), but with positive x. 
We verify the conditions for h on the curve A 3 at the points ( 8 . 8 ). We have 

h-h* = \\ A 2 - 9A 2/3 + 3 
4 

This value is nonnegative for A < (3/2 ) 3 / 4 ~ 1.3554 and hence for A < A* « 1.2408. Let A > A*. Then 
the difference 

h - h** = j 3A 2 + 3A 2/s + 3^3(2- A 4 / 3 ) - (4 + A 4/3 ) 3/2 

must be nonnegative. The substitution X = (A 2 / 3 + yj A + A 4 / 3 ) 2 reduces the condition h > h** to 
the inequality 

(X - 8 ) 3 (A 3 - 12X 2 - 32) < 0, 

whose solution is the interval 

X € [ 8 , 2(2 + 2 2/3 + 2 4/3 )]; 

then A G [A*, A 4 ], where 

r 1 1 4 ^ 4 

A 4 = — (25 - 27 • 2 1/3 + 9 • 2 2/3 ) « 1.2740. ( 8 . 12 ) 

Thus, a pair of intersection points of the curves Ai and A 3 defined by Eqs. ( 8 . 8 ) exists for A 3 < A < A 4 . 
Finally, we obtain the following assertion. 

Theorem 22. In the Kovalevskaya-Yehia case, there are ten structurally stable Smale-Fomenko di¬ 
agrams S' LH ( A). The separating values of the parameter A are 

0 , A 1} A 3 , A*, 1 , A*, A 4 , A 5 = 2^2-1, A 2 , y/2, (8.13) 

where X\, A 2 , A 3 , and A 4 are defined by Eqs. (3.55), (3.56), (8.11), and ( 8 . 12 ). 

The expanded diagram 

MS'lh) = U(£l#(A) x {-M 
a 

generates 29 chambers in the expanded space R. 3 (£,h, A). Moreover, the Smale chambers C, B, E, and 
H have no additional partitions whereas the Smale chambers A, B, E, and H generate the Smale- 
Fomenko chambers Ai-Ai 3 , Bi~B 3 , Ei-Eg, and Hi-H 3 . 

A Smale-Fomenko diagram for small A (0 < A < Ai) is shown in Fig. 54 (with small smooth 
distortions of the picture). Recall that, as is seen on the diagrams of the third critical subsystem (see, 
e.g., Fig. 20), the sets A 3 and h 2 have two common points in the domain t > 0 for small A, namely, 
the points R 3 and R 5 . The curve A 3 passes through the cusp <5 2 of the Smale-Fomenko diagram. 
This is the image of the singular point R 3 ; it has a finite limit as A —>• 0. According to (4.20), the 
L-coordinate of the image of the second point B§ tends to 00 . Therefore, a large distortion is made 
on the curve A 3 for large h in order to show the distant common point of the curves A 3 and h 2 (a 
tangent point) and the boundary between the chambers Ag and Ag. 

In Figs. 55 and 56, changes related to transitions through A = Ai (the chambers B 3 and B disappear) 
and through A = A 3 (the chamber B 2 disappears) are shown. New chambers do not appear here. 

Significant changes in the Smale-Fomenko diagram (as in all separating sets considered above) and 
in chambers occur when transitioning through the value A = A*. At this value, the loop surrounding 
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the chamber Eo, the chamber E3, the swallowtail of the curve 62 along with the chambers A5, Aq, 
A7, and C, and segments of the boundary of the chamber A4 are collapsed and then disappear. For 
A > A*, a swallowtail appears again on the curve 62 , but with other new chambers in its neighborhood, 
namely, the chambers A 10 , An, A 12 , A 13 , Eg, and F. The tansition is shown in Fig. 57. 

The transition through the value A = 1 is quite simple: the common tangent point of A 2 , Ao, and 
A 3 falls on the axis of symmetry t = 0 and then disappears; at the same time, the chambers A 3 and Ei 
also disappear. The transition is shown in Fig. 58. 
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All further surgeries of the Smale-Fomenko diagrams are associated with points on the axis £ = 0. 
The domain containing all these surgeries and the corresponding chambers is shown in Fig. 59. 

The chambers G and Hi appear when transitioning through A*, while the chamber H 2 appears 
when transitioning through A 4 . After the transition through A 5 the chamber A 13 disappears and 
the chamber H 3 appears. After the transition through A 2 , the chamber F disappears, and after 
the transition through the last separating value y/2, the chambers G, Hi, and A 2 disappear. This 
information on the lifetime of chambers is presented in Table 8.1. 

8.2. Fomenko graphs. For each chamber in the space M 3 (-£, h, A) that is cut by the Smale-Fomenko 
diagram, we define a rough topological invariant, the Fomenko graph or the molecule Wi } h{ A), i.e., the 
graph obtained by the retraction of each connected component of the integral manifold, with indication 
of the additional integral of the type I\ for each critical level and, if necessary, the orientation of the 
corresponding atom. Two Fomenko graphs are considered to be matching (identical) if there exists a 
graph homeomorphism that can be prolonged to atoms (for detail, see [ 8 ]). 

All Fomenko graphs that appear in this problem on smooth isoenergetic levels Qf h do not contain 
rank-0 critical points, and degenerate rank-1 critical points are described in [47] and are presented in 
Fig. 60. Here, graphs are collected in groups that are formally identical, but they cannot be converted 
one into another. The groups 1-6 correspond to graphs of the types W\-Wq found in [14]. The 
groups 7-9 are new. In contrast to the group 3 and the type VF 3 , in the groups 7 and 8 atoms B are 
not joined “head to head” (see [14]); in the group 9, in contrast to the graph Wy, the edge from the 
atom C 2 goes to the “foot,” but not to the “head” of the atom B, which generates another Liouville 
foliation (see [14]). 

Naturally, there exist coinciding Fomenko graphs even for different type of isoenergetic surfaces. To 
distinguish these graphs, it is necessary to apply a fine classification (see [ 8 ]). Using information on all 
arising chambers (see Table 8.1), we see that for most graphs, we can obtain marks directly from the 
analogs for the cases A = 0 and £ = 0. We note that in a similar table in [47], there is an inaccuracy: 
the lifetime of the chamber A 2 is set as unlimited, although this chamber and the corresponding graph 
do not exist for A 2 > 2. 

Marks on graphs that have no analogs can be found by using loop molecules of nodal points of 
bifurcation diagrams found above. A possible representation of the result (with gluing matrices) is 
presented in [32]. This problem is discussed in the following section. 
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(a) X<K - (c) X>K 



(a) X<K (b) X=K (c) X>K 


Fig. 57. Transition through A* 
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Table 8.1 


Chamber 

Group of graph 

Lifetime 

Exit on 

Marked 

(number) 

on A 

O 

II 

O 

II 

< 

molecule 

Ar 

1(1) 

0 < A < +00 

Yes/Yes 

A [7, Table 3], [20, Table 8 ] 

A 2 

2 ( 2 ) 

0 < A < y/2 

No/Yes 

B [20, Table 8 ] 

A 3 

3(6) 

0 < A < 1 

Yes/Yes 

C [7, Table 3], [20, Table 8 ] 

A 4 

2(3) 

0 < A < A* 

Yes/No 

B [7, Table 3] 

A 5 

6(9) 

0 < A < A* 

Yes/No 

J [7, Table 3] 

a 6 

7(16) 

0 < A < A* 

No/No 


a 7 

6(9) 

0 < A < A* 

No/No 


a 8 

2(3) 

0 < A < +00 

No/No 


Ag 

8 ( 20 ) 

0 < A < +00 

No/No 


A 10 

2 ( 2 ) 

A > A* 

No/No 


An 

6(10) 

A > A* 

No/No 


A 12 

7(17) 

A* < A < A 2 

No/No 


A 13 

6(11) 

A* < A < A 5 

No/Yes 

F [20, Table. 8 ] 

Bi 

2(4) 

0 < A < +00 

Yes/No 

F [7, Table 3] 


7(18) 

0 < A < A 3 

No/No 


b 3 

6(15) 

0 < A < Ai 

Yes/No 

G [7, Table 3] 

C 

6(12) 

0 < A < A* 

No/No 


D 

6(13) 

0 < A < Ai 

Yes/No 

/ [7, Table 3] 

Ei 

5(8) 

0 < A < 1 

Yes/Yes 

D [7, Table 3], [20, Table 8 ] 

E 2 

4(7) 

0 < A < A* 

Yes/No 

E [7, Table 3] 

e 3 

6(9) 

0 < A < A* 

Yes/No 

H [7, Table 3] 

e 4 

8(21) 

0 < A < +00 

No/No 


E.5 

9(22) 

0 < A < +00 

No/Yes 

E [20, Table 8 ] 

e 6 

6(11) 

A* <C A <! +00 

No/Yes 

G [ 20 , Table 8 ] 

¥ 

6(10) 

A* < A < A 2 

No/No 


G 

6(14) 

A* < A < y/2 

No/Yes 

H [20, Table 8 ] 

Hi 

7(19) 

A* < A < y/2 

No/No 


h 2 

2(5) 

A 4 < A < -|-oo 

No/Yes 

J [20, Table 8 ] 

h 3 

6(14) 

A 5 < A < +00 

No/Yes 

I [20, Table 8 ] 
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Fig. 59. Surgeries at large values of A 


Now we illustrate the operation of the “constructor” of Fomenko graphs by examples. We consider 
the Smale-Fomenko diagrams for average values A (we take A = 0.8). Let h = 2.5 (Fig. 61). If 
l increases from zero, this level crosses five chambers E 5 , E 4 , Ag, As, Ai (the paths A = const on 
Fig. 61 are marked by the numbers 1,... ,5). In the corresponding /i-section of the diagram £(A), 
Fomenko graphs are determined by the bifurcations along the lines £ = const for increasing k (five 
broken arrows in Fig. 62). Another eight chambers can be observed on the same level of h at small 
values of A. The Smale-Fomenko diagram for A = 0.1 and the level h = 1.8 is shown in Fig. 63. The 
corresponding paths t = const on the diagram £(A) are shown in Fig. 64. The crossed arcs correspond 
to domains determined by the diagrams of critical subsystems. The sequences of these intersections 
and the corresponding atoms for paths marked by the numbers 1,..., 13 are shown in Table 8.2. 

The following interesting phenomenon can occur (in the study of the coincidence of Fomenko graphs, 
this phenomenon is usually not discussed): some groups contain graphs that differ only in a pair of 
atoms on the same critical level. For example, in Fig. 60, all levels containing two critical circles 
possess the following property of “stability”: for any sufficiently small perturbation {&, h ), the number 
of circles on such critical level is invariant. However, this is not so if the critical level contains multiple 
points, for example, the levels i = 0 and k = l + (/r—A 2 /2), where h > A 2 /2, i.e., all levels whose images 
lie on the singular parabola (4.4). Such a level is contained in each Fomenko graph of the form Wq^ 
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Fig. 60. Fomenko graphs 
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Fig. 62. Bifurcation diagram in the (£, fc)-plane for A = 0.8 and h = 2.5 
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Fig. 64. Bifurcation diagram in the (£, £:)-plane for A = 0.1 and h= 1.8 

































Table 8.2 


Number of 

path 

Chamber 

Sequence 

of arcs 

Sequence 

of atoms 

Group of graph 

(number) 

1 

E5 

b\ —y C2 —y clq —y 

—y C4 —y a2 —y cti 

A+ — > B + — ► (A + , A_) — >■ 

->• C 2 ->• A_ ->• A- 

9 ( 22 ) 

2 

e 4 

b\ —y C2 —y clq —y 

—y 0-3 —y cl\ 

A+ — > B + — > (A +1 A_) — > 

B- -> A_ 

8(21) 

3 

Ag 

b\ —y C2 —y ci —y 

—y a4. —y CL3 —y cl\ 

A- 1_ —y B- 1_ —y A— —y A-y. —y 

-> B_ -> A- 

8(20) 

4 


b\ —y CL4 —y CL3 —y cl\ 

A-y. —y A-y. —y B _ — y A— 

2 ( 3 ) 

5 


b\ —y ci\ 

A+ -»• A- 

1(1) 

6 


b 2 ->• a 3 -> c 6 -» 

—y C7 —y cl\ 

2A + —y B— —y A-y. —y 

-> B_ -> A. 

7 ( 18 ) 

7 

Ag 

b\ —y a4 —y CL3 —y 

—y cq —y C'j —y &\ 

A-y. — y A-y. — y B— — y 

A+ -> B- -> A- 

7 ( 16 ) 

8 

C 

b\ — y cl4 — y C3 —y 

—y clj —y C'j —y cl\ 

A-y. — y A-y. — y 2A-y. — y 2B— —y 

-> B_ -> A. 

6(12) 

9 

A ,5 

b\ —y cq —y cLg —y 

—y cl'j —y cj —y cl\ 

A-y. — y A-y. — y 2A-y. — y 2B— —y 

-> B_ -> A. 

6 ( 9 ) 

10 

e 3 

b\ — y Cq —y clq —y 

—y eg —y clq —y cl\ 

A-y. — y A-y. — y 2A-y. — y 2B— —y 

-> B_ -> A- 

6 ( 9 ) 

11 

E2 

b\ —y cq —y clq —y 

—y clq —y cl\ 

A+ ->• A + -> 2 A* -> 

-> B- -> A- 

4 ( 7 ) 

12 

Ei 

b\ — y cq — y clq — y 

—y C4 —y ci2 —y cl\ 

A+ ->• A + -> 2 A* 

C 2 A- -> A- 

5 ( 8 ) 

13 


b\ —y cq —y cj —y cl\ 

A.y. — y A-y. — y B — y A— 

2 ( 3 ) 
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Table 8.3 


Chamber 

Sequence 
of arcs 

Sequence 
of atoms 

Group of graph 
(number) 

Mh) 

^ 

A+ A_ 

1(1) 

A 2 (/ 2 ) 

^ CZ0 —^ ^2 —^ &1 

A4 — y B- 1_ — y A— —y A— 

2(2) 

A 3 (/ 3 ) 

&1 —?*■ C0 —^ C7 —^ 

—^ 0-6 —?*■ (22 —^ &1 

-A_|_ —y A4. —y B _ —y 

—)► B+ —>• A- —>• A- 

3 ( 6 ) 

a 4 (/ 4 ) 

^ C0 —^ C7 —^ 

A4 — y A4 — y B— — y A— 

2 ( 3 ) 


^ C0 —^ 0-8 —^ 

—?*■ (27 —^ C7 —^ ( 2 i 

A- 1_ — y A4 —y 2A4 — y 

2B_ -> -» 

6 ( 9 ) 

A 6 (/ 6 ) 

61 —?*■ cz.4 —^ 0-3 — y 
— y cq — y C7 — y &\ 

A4 — y A4 —y B— — y 
-> A + B- ->• A- 

7 ( 16 ) 

At(/ 7 ) 

b\ —y (24 —y cq —y 
—y eg —y clq —y ci\ 

A4 —y A4 —y 2A4 —y 

2 B_ -> B_ -> A- 

6 ( 9 ) 

A 8 (/ 8 ) 

b\ —y 04 —y clq —y cl\ 

A4 —y A4 —y B— —y A— 

2 ( 3 ) 

Ag(/g) 

b\ —y C2 —y ci —y 
—y a4. —y cl 3 —y cl\ 

A- 1_ —y B4 —y A— —y 
-> A4 B- ->• A- 

8(20) 

Aio(/io) 

b\ —y C2 —y C\ —y cl\ 

A4 —y B- 1_ —y A— —y A— 

2(2) 

An(/n) 

b\ —y C2 —y cm —y 

—^ ®10 —^ Ci —y CL\ 

A + -q> B + ->• 2 B + ->• 

-q> 2 , 4 _ ->• -> T_ 

6(10) 

Al2(/l2) 

b\ —y clq —y 0,2 —^ 

—y C2 —y c\ —y ci\ 

A4 —y B4 —y A— —y 
—)► B4 —>• A- —>• A- 

7 ( 17 ) 

Al 3 (/l 3 ) 

b\ —y (20 —y C5 —y 
—y C3 —y 0-2 —y a 1 

A + -> B + -> 2 B + 

-> 2^1_ ->• A- 

6(11) 

®i(/u) 

b 2 —> 03 —> a 1 

2 A + -> B- ->• 

2 ( 4 ) 

B 2 (/is) 

b2 -» a 3 c 6 -» 

—y cj —y ci\ 

2A4 —y B— —y A4 —y 

B- A- 

7 ( 18 ) 

B 3 (/i 6 ) 

l >2 Cg —> Cg —> 

—y (13 —y o,\ 

2 A + ->• 2 A + 2 £?_ 

—> B_ —> A- 

6 ( 15 ) 

C(/it) 

b\ —y cl4 —y c§ —y 
—y cij —y C7 —y cl\ 

A4 —y A4 —y 2A4 — y 
-¥ 2 B_ -> B- -> 

6(12) 

©(/is) 

62 —^ cs —y cij —y 
—y C7 —y ct\ 

2 A + ->• 2 T + -> 2 B_ -> 

6 ( 13 ) 

Ei(/i 9 ) 

b\ —y cq —y clq —y 
—y c\ —y cl2 —y cl\ 

A + A + 2 A* -> 

-> C 2 -> A_ A_ 

5 ( 8 ) 

E 2 (/ 2 o) 

b\ —y cq —y clq —y 
—y clq —y cl\ 

A + -> A + -> 2 A* 

B_ — >■ A- 

4 ( 7 ) 

E 3 (/ 2 l) 

b\ —y cq —y clq —y 
—y eg —y clq —y cl\ 

A4 —y A4 —y 2A4 —y 
-> 2 B_ -> B- -> A- 

6 ( 9 ) 

E 4 (/ 22 ) 

b\ —y C2 —y cLg —y 
—y clq —y cl\ 

A+ —> B + — > (. A + , A- ) —> 
B- -> T_ 

8(21) 

E 5 (/ 23 ) 

b\ —y C2 —y cLg —y 
—y C4 —y cl2 —y &\ 

A+ —t B + — > (T + , A_) —>■ 

->• C 2 ->• A- 

9 ( 22 ) 

e 6 (/ 24 ) 

b\ —y C2 —y cLu —y 

—y Cq —y (22 —y CLi 

A + -> B + 2 B + -> 

2 ->■ A- A^ 

6 (11) 








































































































Table 8.3 (continued) 


Chamber 

Sequence 
of arcs 

Sequence 
of atoms 

Group of graph 
(number) 

F(/ 25 ) 

b\ —y clq —y C5 —y 
—y ciio —^ ci —y ci\ 

A+ B+ 2 B + ->■ 

-s> 2^4 T A 

6(10) 

g(/ 26 ) 

b\ —y clq —y C5 —y 
y ciio ~^ a i2 

A + -> B + -> 2 B + ->• 

->■ 2A 

6 ( 14 ) 

Hi(/ 27 ) 

b\ —y clq —y ci2 —y 
—>■ C2 —>• a±2 

.A_|_ —y B- 1_ —y A— —y 
-y B+ 2 A- 

7 ( 19 ) 

h 2 (/ 28 ) 

bi -» c 2 -» ai2 

A + -> B + -s> 2 

2 ( 5 ) 

H 3 (/ 29 ) 

bi —> c 2 —> an —»■ 

—>■ aio —> 012 

A + -> B + -> 2 B + ->• 

->■ 2T_ -)• A^ 

6 ( 14 ) 


with h > A 2 /2 (recall that the boundary value h = A 2 /2 determines an isoenergetic manifold with a 
degenerate point; we do not consider this case in the present paper). It is easy to see that under a 
small perturbation of the zero level of t within the chamber, the instable level splits into two levels 
and its atoms diverge. This phenomenon occurs in all chambers that have an exit on the axis l = 0 
for h > A 2 / 2 . Unstable (in the above sense) Fomenko graphs are presented in Fig. 65 , where the mark 
of the chamber equipped with the index 0 means the intersection of the chamber with the axis i = 0. 



Fig. 65 . Unstable Fomenko graphs 


8.3. Fomenko Zieschang invariants. To complete the description of topological invariants of 
isoenergetic levels of reduced systems, we must equip the Fomenko graphs found above with marks. 

First, we write all transitions between atoms that can occur in Fomenko graphs into a table (see 
Table 8 . 4 ). Of course, here we mean atoms in which a wall of the bifurcation diagram (or, which is the 
same, a domain of the corresponding critical subsystem) is indicated. For each transition, we indicate 
all graphs in which this transition occurs. If at least one of these graphs is preserved in the limit 
cases A = 0 and i = 0 , we note this fact in the corresponding column of the table and then take r- 
and e-marks from the results of [ 7 ] or [ 20 ], respectively. Missing marks are taken from loop molecules 
listed in the corresponding column of the table. In this case, necessary comments are given below. 
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Table 8.4 


Number 
of edges 

Edge 

(r, £)-labels 

Graphs 

Exit on 

A = 0/£ = 0 

Loop 

molecule 

1 


r = 0, e = 1 

Ai 

Yes/Yes 


2 

(^[i>i] ’ ^[ae]) 

r = 0, e = 1 

A 2 , A12, A13, 

F, G, Hi 

No/Yes 


3 

(-®[a s ] ) A [a 2 ] ) 

r = 0, £ = 1 

A 2 , A3, A12, Hi 

Yes/Yes 


4 

(■®[as]’ ^[01]) 

r = 0, £ = 1 

A 2 , A3 

Yes/Yes 


5 

(A-W, B [C2] ) 

r = 0, £ = 1 

Ag, A10, An, E4, 
E5, Eg , H 2 , H3 

No/Yes 


6 

(B[c 2 ]', A [c 1]) 

r = 0, £ = 1 

Ag, A10, Ai 2 

No/No 

A31 

7 

(73[c 2 ] 1 4[ai] ) 

r = 00, £ = 1 

A10, Ai 2 

No/No 

R" 

°27 

8 

( B [c 7 ]b 4 [ai ]) 

1 

r= 2’£ = 1 

A4, A5, A6, 

©2, C, ID 

Yes/No 


9 

(T [C6 ];B[ C7 ]) 

r = 0, £ = 1 

A3, A4, b 2 

Yes/Yes 


10 

(^[&i]:^[c 7 ]) 

r = 0, £ = 1 

A3, A4 

Yes/Yes 


11 


r = 0, £ = 1 

A7, Ag, Ag, Bi, 
B3, E 2 , E3, E4 

Yes/No 


12 

(^[a 4 ];S[as]) 

r = 0, £ = 1 

Ag, Ag, Ag 

No/No 

Ai 2 

13 

{A[ bl y, B [a3 ]) 

r = 00, £ = 1 

Ag, Ag 

No/No 

^31 

14 

{ 2A [b2]'i B[a 3 ]) 

r = 00, £ = 1 

Bi, B 2 

Yes/No 


15 

(B [c2 y,2A [ai2] ) 

r = 00, £ = 1 

Hi, H 2 

No/Yes 


16 

(B[ C7 ]; B [ae] ) 

r = 0, £ = 1 

A3 

Yes/Yes 


17 

(^[05 ]’ -®m) 

r = 0, £ = 1 

e 2 

Yes/No 


18 

(^[c 6 ]! ^[a 5 ]) 

r = 0, £ = 1 

Ei, E 2 

Yes/Yes 


19 


r = 0, £ = 1 

Ei, E 2 

Yes/Yes 


20 


r = 0, £ = 1 

Ei, E5 

Yes/Yes 


21 

(C2 [c 4 ];^[a 2 ]) 

r = 0, £ = 1 

Ei, E5 

Yes/Yes 


22 

(^m;C 2W ) 

r = 0, £ = —1 

Ei 

Yes/Yes 


23 

(2 B [o 7] ;B [c7 ]) 

r = 0, £ = 1 

A5, C, ID 

Yes/No 


24 


r = 0, £ = 1 

A5 

Yes/No 
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Table 8.4 (continued) 


Number 
of edges 

Edge 

(r, £)-labels 

Graphs 

Exit on 

A = 0 /£ = 0 

Loop 

molecule 

25 

(A[c 6 ] ) -®[a 7 ]) 

r = 00, £ = 1 

A5 

Yes/No 


26 

(d-[f,i], -B[a 7 ] ) 

r = 00, £ = 1 

A5, C 

Yes/No 


27 

(2 B [cg] ;B [as] ) 

r = 0, £ = — 1 

A7, B3, E3 

Yes/No 


28 

( 2-4 [c 8 ]; 2 B[c 9 ]) 

r = 0, £ = 1 

A7, B3 

Yes/No 


29 

(^[04] ) -®[cg]) 

r=i£=i 

A 7 

No/No 

A03 

30 

(A [bl] ; B [cg] ) 

r = 0, £ = 1 

A 7 , Eg 

Yes/No 


31 

( 2 ^ 4 [a 8 ]; 2 -B[ C9 ]) 

r=i£=i 

E3 

Yes/No 


32 

(d-[ C6 ]; 73 [c 9 ]) 

r = 0, £ = 1 

E3 

Yes/No 


33 

(B [ail] , ^ 4 [ai]) 

r = 0, £ = 1 

An, Fg 

No/Yes 


34 

(B [oil] ;A [ci] ) 

r = 0, £ = 1 

An 

No/No 

A04 

35 


r = 0, £ = 1 

An, H 3 

No/Yes 


36 

(B [c2] ;2-B [oil] ) 

r = 0, £ = — 1 

An, Eg, H3 

No/Yes 


37 

(^[Cb]; ^[04]) 

r = 00, £ = 1 

A13, E 

No/Yes 


38 

( 5 [c 5 ];- 4 [ci]) 

r = 0, £ = 1 

F 

No/No 

A31 

39 

( B [c 5 ]; 2 ^[aio]) 

r = 00, £ = 1 

G, F 

No/Yes 


40 

( B [a 6 ];2-B[ C5 ]) 

r = 0, £ = 1 

A13, G, F 

No/Yes 


41 

( B [c s y, A[a 2 ]) 

r = 00, £ = 1 

A13 

No/Yes 


42 

(2B [cb] ;2A [03] ) 

r = 0, £ = 1 

A13 

No/Yes 


43 

(B[a llh A[ a2 ]) 

r = 0, £ = 1 

Eg 

No/Yes 


44 

(B[au] > ^-^[03]) 

r = 0, £ = 1 

Eg 

No/Yes 


45 

(2^[c 8 ];-B [07] ) 

r = 00, £ = 1 

C, D 

Yes/No 


46 

(2A [62 ];B[ 07] ) 

r = 00, £ = 1 

D 

Yes/No 


47 

(B [cb ];2A[ 0 i 2 ]) 

r = 00, £ = 1 

G 

No/Yes 


48 

(-^[an] 5 2-A[ ai2 ]) 

r = 0, £ = 1 

h 3 

No/Yes 
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Table 8.4 (continued) 


Number 
of edges 

Edge 

(r, £)-labels 

Graphs 

Exit on 

A = 0/£ = 0 

Loop 

molecule 

49 

( 2 A [b2] - 2 B [c9] ) 

r = 0, e = 1 

®3 

Yes/No 


50 

(-®[c 6 ]> -®[c 7 ]) 

r = 0, e = 1 

Ag 

No/No 

A33 

51 

( B [a 3 y, B [c 7 ]) 

r = 0, £ = 1 

Ag, ®2 

No/No 

822 

52 

( B [a 6 ]i B [C2] ) 

r = 0, £ = 1 

A12, Hi 

No/No 

<526 

53 

(B[c 2]) B[a 3 ]) 

r = 0, £ = —1 

Ag, E4 

No/No 

<531 + • • • 

54 

(A[a 9 ];P[o 3 ]) 

r = 0, £ = 1 

e 4 

No/No 

(531 + (532 

55 

(B[ C2 y, A[ ag] ) 

r= \' e=1 

E 4 , E5 

No/Yes 


56 

(A[ 09 j; C2[c 4 ]) 

r = 0, £ = 1 

E5 

No/Yes 


57 

(B[ ca ]-,C 2 [ C4 ]) 

r = 0, £ = —1 

E5 

No/Yes 



Most of the marks on edges that are not present in the limit cases need no explanations. Thus, the 
edge No. 6 is explicitly present in the loop molecule of the point A31 (see Fig. 46 ), the edge No. 12 in 
the loop molecule of the point A12 (see Fig. 43 ), the edge No. 50 in the loop molecule of the point A33 
(see Fig. 47 ), the edge No. 7 in the loop molecule of the point 8'^, and the edge No. 13 in the loop 
molecule of the point £31 (the last two molecules are shown in Table 6 . 2 ). The edge No. 51 is explicitly 
present in the loop molecule of the nondegenerate saddle point 822 and the edge No. 52 in the loop 
molecule of the nondegenerate saddle point 826 (both molecules are given in Table ??). 

Marks on some transition edges can be set immediately if we note that, in fact, these transitions 
are reduced to already known ones. So, consider the edge No. 54 . On the diagrams Hhk(£, A), it 
joins the segments <23 and ag. Consider the diagram for the domain 9 (see Fig. 28 ). It contains the 
images of nondegenerate rank 0 critical points of the classes £25 and £32 (of the type “center-center”) 
and (53i (of the type “center-saddle”). If we slightly spread out the “double” edges (whose preimages 
contain the two critical circles), we obtain the picture shown in Fig. 66; then it becomes obvious that 
the transition discussed is actually a transition from <23 to a part of ag marked by a'g in the figure: it 
is simply a prolongation of the segment 04. Therefore, this transition coincides with the edge No. 12 
whose mark has just been found. Note that in the proof of Proposition 15 (see above) this result for 
the transition along the edge No 12 was proved by using the addition rule for marks (Lemma 4 ) along 
two transitions a'g = a 4 —>• b” with r = 0 from the loop molecule of the point £32 and b 2 —>• 03 with 
r = 00 from the loop molecule of the point <531 (they are shown by broken arrows). 

Similarly, we consider the edge No. 53 and, for example, find it on the diagram T,lk{K, A) for the 
domain 33 (see Fig. 52 ). As in the previous case, we spread out “double edges” (in this case, the 
same edges ag and 62); the result is shown in Fig. 67 . Obviously, the transition discussed C2 -Y 03 is 
the composition of two transitions: the transition from C2 to the part of the chamber 62 indicated in 
the figure by b 2 (it is a prolongation of the segment b \; this is the edge No. 5 with r = 0 ) and the 
transition from b\ to a% (the edge No. 1 , which is a part of the loop molecule of the nongenerated 
point (53i of the type “center-saddle” with the mark r = 00). These two transitions are shown in the 
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figure by broken arrows. Thus, we obtain that the r-rnark on the transition discussed is equal to zero 
by Lemma 4 . 

Consider the edge No. 29 . It is convenient to examine it on the diagram El#(^L A) for the domain 31 
(see Fig. 52 ). Spreading out the double edges as and eg, we take into account that these edges start 
at a degenerate point A03, and both edges c 9 and c 9 terminate at the corresponding points A32 and 
contain nondegenerate points of the type “center-center”: 824 (the end of one edge as, for example, a' 8 ) 
and 832 (the end of the other edge ag, which coincides with 04). Thus, the transition 04 —>• eg is one of 
the two transitions as —» eg, which, according to the loop molecule of the degenerated point A03 (see 
Fig. 39 ) have r = 1 / 2 . The same result directly follows from the well-known mark of the “double” 
edge No. 31 in the graph E3, which exists for A = 0 (see [ 7 ]). 

Note that the analysis of the last two cases imply that 04 = a§ and 04 = a 9 . This fact is not casual: 
it can be clearly seen in the diagrams of the critical subsystem M. For example, in Fig. 15 , for 
small A the domain 04 is separated from as by the curve consisting of the points 824 and from ag by 
the curve consisting of the points 825■ Both these curves correspond to points of the type “center- 
center” (so they generate a new critical circle), whereas the critical circle 04 just passes through it 
without bifurcation. 

On the edge No. 34 , the behavior of transitions is determined by the loop molecule of the point A04 
(see Figs. 40 and 41 ). Since the preimage of the point 8 ' 28 contains a nondegenerate critical point of 
the type “center-center,” the arc C3 on the diagram A) for the domain 24 splits into two arcs 

(see Fig. 52 ): the first arc C3 terminates at the point S' 28 , while the second can be prolonged to the arc 
C3 = ci- Hence, the transition an —>• c\ is one of the two transitions an —>• C3 in the first component 
of the loop molecule of the point A04, where r = 0. 
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Fig. 67 . Diagram S^(/i, A) for 33 region with layers indicated 



Fig. 68. Fragment of the diagram A) for the domain 31 with layers indicated 


Finally, another transition that has no direct analogs in the limit cases studied above is the edge 
No. 38 . Consider the diagram Y,LK(h, A) for the domain 7 (see Fig. 52 ) in a neighborhood of neighbor¬ 
ing nodal points A31 and 5 ' 27 . The latter point has the type “center-saddle” and attracts the part c' 5 
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of the double arc C5 that terminates at this point. The transition to the arc c\ is possible only from 
the arc c' 5 ', which, as a result, can be prolonged to C2. Thus, the transition C5 —>• c\ is the transition 
C2 —> c \, which is contained in the loop molecule of the point A31 (see Fig. 46 ). Hence, according to 
what was proved above, we obtain r = 0. 

Here we have described the process of producing marks using either transitions in the special cases 
A = 0 and t = 0 examined earlier, or loop molecules constructed above; for brevity, we have discussed 
only r-marks. In fact, the information used also contains all e-marks. Thus, all data in Table 8.4 are 
completely justified. 

Table 8.5 contains all Fomenko-Zieschang invariants for the Kovalevskaya-Yehia case. Each of them 
is supplied with the notation of the graph and the corresponding isoenergetic manifold. The complete 
list coincides with the list given in [ 32 ]; however, there is no correspondence between marked molecules 
and topological types Q^ h in [ 32 ] (although the notion of invariant requires such a correspondence). 

After comparing invariants found (supplied with (r, e)-marks defined above) with invariants listed 
in [ 32 ], we transferred all n-marks of families of hyperbolic atoms from [ 32 ]. Calculation of these marks 
requires considerable additional technical work. So, in [ 32 ], using one of the main procedures of the 
method of loop molecules, namely, the expression of all cycle bases on families of regular tori through 
A-cycles of atoms of bifurcations, which are defined unambiguously (see [ 7 ]), a method of arrangement 
of glue matrices along edges of marked graphs is proposed. Then all n-marks of families are calculated 
by the corresponding definitions. 


Table 8.5 


Marked isoenergetic 
molecule 


Marked isoenergetic 
molecule 


A, 




Ar. 


r = 0 
8 = 1 


Ar, 


A, 





r = 0 
s = 1 
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Table 8.5 (continued) 
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Table 8.5 (continued) 
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Table 8.5 (continued) 
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Table 8.5 (continued) 
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Table 8.5 (continued) 
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Table 8.5 (continued) 
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Table 8.5 (continued) 



9. Conclusion 

The study of the phase topology of the Kovalevskaya-Yehia case is complete. We have presented 
all results available today on analytical solutions and topological analysis of the Kovalevskaya-Yehia 
case. As was noted above, this work can be extended by a description of limit cases, namely, of 
the classical Kovalevskaya problem and a gyrostat with zero area constant. Here these results are 
assumed to be known and are used in the construction of fine topological invariants. The title of 
this paper contains the words “topological atlas”; we do not give any strict mathematical sense to 
this term. The authors want to emphasize that all objects studied are classified in the spaces of the 
corresponding parameters; separating sets and stability domains are constructed; the results of these 
studies resemble geographical atlases. We did not give a detailed explanation of the general theory 
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of integrable Hamiltonian systems (except for the very minimum necessary for definitions), since the 
volume of this work would substantially increase. 
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